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1  INTRODUCTION 


The  general  goal  of  this  research  is  to  build  an  experimental  system  based  on  Ultrasound 
Stimulated  Acoustic  Emission  (USAE)  for  in  vitro  imaging  of  breast  tissue  specimens  and  test  its 
performance  in  detecting  microcalcification.  The  general  theory  of  USAE  is  described  in  [1,2, 3,4,5] 
and  some  applications  are  explained  in  [6,7,8,9,10].  The  first  task  of  this  research,  which  covers 
the  first  year  and  a  part  of  the  second  year,  focuses  on  system  development  and  optimization  for 
detection  of  beast  microcalcification.  The  second  task,  covering  parts  of  the  second  year,  is  centered 
on  tissue  scanning  and  studying  the  manifestation  of  microcalcifications  in  USAE  images. 

This  report  addresses  the  activities  performed  in  the  second  year  in  both  areas.  In  particular,  the 
report  includes  experimental  results  demonstrating  the  capability  of  USAE  method  in  imaging  and 
delineating  calcifications  in  breast  tissue. 


2  BODY 

In  this  section  we  describe  research  accomplishments  associated  with  parts  of  Tasks  1  and  2. 
First,  we  report  system  development  activities  that  improve  our  measurements.  Next,  we  report  some 
experimental  results  on  human  tissue  specimens.  We  will  show  and  discuss  USAE  images  obtained 
from  tissue  samples  with  microcalcifications  and  compare  the  results  with  the  corresponding  x-ray 
mammography  histology  images. 

2.a  System  Development 

2.a.l)  Quiet-scanning  System 

The  step  motors  used  in  the  present  scanner  produce  loud  audio  noise  that  can  interferes  with 
the  acoustic  emission  produced  by  the  object.  To  reduce  the  vibration  and  acoustic  noise  of  the 
scanning,  a  new  scanning  system  has  been  designed  that  is  based  on  analog  servo  motors  instead  of 
the  present  step  motors.  The  new  scanning  system  has  optically  encoded  axes  for  precise  positioning. 
Construction  of  the  new  system  has  been  successful.  The  system  is  now  under  evaluation  test. 


2.a.2)  Software  development  for  the  new  scanning  system 

The  software  for  controlling  the  new  scanning  system  and  the  data  acquisition  system  has  been 
developed. 

2. a.  3)  Direct  Displacement  Measurement  by  Laser  Vibrometer 

Recently  we  purchased  a  laser  vibrometer  for  direct  measurement  of  velocity  and  displacement 
in  the  order  of  nano  meters  at  frequencies  up  to  250  kHz  range.  This  system  allows  us  to  directly 
measure  the  vibration  that  we  introduce  in  an  object  by  the  radiation  force  of  ultrasound.  Hence,  we 
can  validate  our  acoustic  measurements  by  direct  optical  measurement. 
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2.b  Breast  Tissue  Imaging  and  Detection  of  Microcalcifications 

The  hypothesis  is  that  ultrasound  stimulated  acoustic  emission  can  be  used  to  image  tissue  and 

detect  microcalcifications.  To  test  this  hypothesis,  we  scanned  several  excised  human  breast  tissue 

and  compared  the  results  with  the  x-ray  mammography  and  histology. 

2.b.l)  Procedure 

Experiments  were  conducted  on  human  breast  tissue  samples  in  the  following  manner: 

1 .  Initial  tissue  selection  —  Breast  tissues  were  enlisted  within  a  few  days  of  surgery  from  Surgical 
Pathology  section  of  Rochester  Methodist  Hospital.  Patients’  records  were  reviewed  to  identify 
those  with  high  probability  of  having  microcalcifications.  Tissues  from  this  group  of  patients 
were  collected. 

2.  Initial  tissue  preparation  —  Tissues  were  cut  into  (  approximately)  3x3  cm  pieces  and  arranged 
flat  and  side-by-side  in  a  plastic  bag.  The  bag  was  constrained  in  a  rectangular  bracket  to  keep 
them  from  moving.  Each  bag  could  contain  several  tissue  pieces,  but  all  were  from  the  same 
patient. 

3.  Screening  mammography  —  High-resolution  x-ray  mammography  of  each  tissue  bag  was  ob¬ 
tained  using  a  tissue  mammography  machine. 

4.  Reading  —  Each  mammogram  was  read  by  our  radiologist  (Dr.  Morton)  to  identify  microcal¬ 
cifications. 

5.  Final  tissue  selection  —  Tissue  pieces  identified  in  the  above  step  were  selected.  Remaining 
tissue  pieces  were  put  away. 

6.  Mounting  —  Each  tissue  piece  selected  in  the  above  step  was  mounted  on  a  Scanning  Bracket. 
This  bracket  is  designed  to  hold  the  tissue  piece  in  the  water  tank  for  acoustic  scanning.  Each 
bracket  is  shaped  like  a  ring,  approximately  10  cm  in  diameter,  and  made  from  plexiglass  material. 
The  ring  is  covered  by  a  thin  latex  sheet.  The  tissue  piece  was  mounted  on  this  sheet  with  a 
few  small  drops  of  glue.  Glue  drops  were  carefully  placed  on  locations  that  would  not  affect 
the  region  of  interest  with  microcalcifications.  To  identify  the  position  and  orientation  of  the 
object,  small  pieces  of  suture  were  used  as  identification  marks.  For  this  purpose  small  knots  are 
placed  at  different  position  on  the  tissue.  These  markers  normally  can  be  seen  in  the  acoustic, 
x-ray,  and  photography  images. 

7.  Base  Mammography  —  Tissue  mammography  images  were  obtained  from  each  mounted  sample. 
These  images  were  used  as  a  base  of  comparison  with  the  USAE  image.  These  images  were 
converted  into  digital  form  and  stored  in  the  computer. 

8.  Scanning  —  Mounted  samples  were  scanned  in  the  water  tank.  Acoustic  images  were  obtained 
at  various  difference  frequencies.  These  images  were  displayed/printed  in  the  same  size  and 
orientation  as  the  x-ray  mammography  images. 

9.  Image  evaluation  —  Acoustic  emission  images  were  evaluated  against  the  x-ray  mammography 
to  assess  the  manifestation  of  calcifications  in  the  acoustic  images. 

10.  Histology  —  Sections  of  the  tissue  identified  with  microcalcification  were  cut  for  histology.  The 
results  were  used  to  validate  the  existence  of  microcalcifications  and  determine  their  size  and 
locations. 
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2.b.2)  Experiment  Results 


To  demonstrate  the  procedure  outlined  above,  we  present  results  of  some  tissue  experiments. 

The  first  example  illustrates  detection  of  microcalcifications.  Figure  1  shows  high-resolution  x- 
ray  image  of  a  group  of  tissue  samples  collected  from  a  57  year  old  subject.  This  x-ray  was  done  for 
screening  (step  3  of  the  list).  It  shows  several  bright  spots  which  are  identified  as  microcalcifications 
(<400  microns)  by  the  radiologist.  The  tissue  piece  that  is  identified  with  microcalcifications  was 
chosen  for  further  scanning.  Figure  2  shows  the  photograph  of  the  selected  piece  mounted  on  the 
scanning  bracket.  Sutures  that  are  knotted  on  the  tissue  are  for  identification  purposes.  Figure  3 
shows  the  mammogram  of  the  selected  tissue  pieces  mounted  on  the  scanning  bracket.  This  sample 
was  scanned  in  the  water  tank.  The  acoustic  emission  image  of  the  tissue  at  26  kHz  is  shown  in 
Figure  4.  Microcalcifications  can  be  seen  in  this  image  as  four  bright  spots.  It  is  interesting  to  note 
that  the  number  and  position  of  these  spots  matches  with  those  microcalcifications  identified  in  the 
x-ray  mammography.  Also  noticeable  is  that  tissue  inhomogeneities  appear  dim  and  do  not  interfere 
with  the  calcifications  in  the  acoustic  image. 

The  second  example  illustrates  detection  of  calcification  in  a  vessel.  Figure  5  shows  a  high- 
resolution  x-ray  image  of  a  group  of  tissue  collected  from  a  78-year  old  subject  after  breast  removal. 
This  is  the  screening  x-ray.  It  shows  several  bright  spots  which  are  identified  as  calcifications  by  the 
radiologist.  The  tissue  piece  that  is  identified  with  having  a  calcified  vessel  was  chosen  for  further 
scanning.  We  have  knotted  sutures  on  three  locations  on  the  tissue  for  identification  purpose.  Figure 
6  shows  the  mammogram  of  the  selected  tissue  pieces  mounted  on  the  scanning  bracket.  The  sutures 
are  barely  visible.  This  sample  was  scanned  in  the  water  tank.  The  acoustic  emission  image  of 
the  tissue  is  shown  in  Figure  7.  Calcification  can  be  clearly  seen  in  this  image.  It  is  interesting 
to  note  that  the  shape  of  calcified  region  matches  with  great  detail  with  the  corresponding  region 
in  the  x-ray  mammography. 

To  validate  the  presence  of  calcifications,  the  tissue  was  examined  histologically.  The  portion 
of  the  tissue  sample  used  for  this  purpose  is  shown  in  the  x-ray  image  in  Figure  8.  Figure  9  shows 
the  histology  of  the  mid-portion  of  the  selected  tissue  piece  in  a  plane  tangential  to  the  vessel.  The 
histology  shows  the  cross  section  of  the  vessel  with  some  irregular  calcifications  on  the  wall,  the 
thickness  of  the  calcification  layer  is  less  than  100  microns. 

2.c  Discussion 

Experimental  USAE  images  presented  in  this  report  demonstrates  two  important  facts:  (1)  USAE 
imaging  method  is  capable  of  detecting  small  microcalcifications,  as  well  as  calcifications  in  vessels, 
in  the  breast  tissue;  (2)  Microcalcifications  can  be  delineated  form  tissue  inhomogeneities.  Results 
shown  here  also  indicate  that  such  images  have  high-spatial  resolution,  high  detectability,  good 
contrast,  and  high  signal-to-noise  ratio.  Microcalcifications  shown  in  our  results  are  only  a  few 
hundred  microns  in  diameter,  and  the  calcifications  in  the  vessel  wall  were  not  more  than  100  microns 
thick.  Histology  image  from  the  second  sample  also  revealed  small  (<40  microns)  microcalcifications 
below  the  vessel.  These  small  microcalcifications  could  not  be  detected  by  either  USAE  or  x-ray 
mammography  methods. 

A  comparison  between  the  USAE  image  Figure  7  and  the  corresponding  mammography  Figure  3 
reveals  another  interesting  capability  of  USAE  method.  High  tissue  density  in  the  region  surrounding 
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the  microcalcification  absorbs  a  great  deal  of  x-ray  energy,  hence  the  x-ray  image  shows  very  low 
contrast  within  this  region,  and  as  a  result,  microcalcifications  are  only  barely  visible.  However, 
the  USAE  image  delineates  these  microcalcification  with  high  contrast.  This  is  because  the  acoustic 
emission  is  sensitive  to  the  elastic  properties  of  the  object. 

Acoustic  images  reported  here  were  obtained  at  relatively  high  vibration  frequencies  (about  20 
kHz).  At  these  frequencies  we  are  able  to  acquire  the  data  at  much  higher  speed,  reducing  the 
acquisition  time  to  a  few  minutes  per  image  (for  more  details,  refer  to  High-Speed  Scanning  section 
in  the  last  year’s  report).  This  is  because  at  these  frequencies  the  noise  level  from  the  scanner 
mechanism  is  negligible,  hence  we  can  keep  the  scanning  motor  running  while  collecting  the  data. 
Scanning  speed  is  an  important  factor  in  our  experiments. 


2.d  Key  Research  Accomplishments 

•  Demonstrated  that  small  microcalcification  (also  vessel  calcifications)  can  be  detected  by  USAE 
imaging  method. 

•  Demonstrated  that  breast  microcalcifications  can  be  delineated  from  tissue  inhomogeneities  in 
a  USAE  image. 

•  Demonstrated  that  microcalcification  can  be  imaged  at  high  vibrational  frequencies. 


2.e  Reportable  Outcomes 
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2.f  Conclusions 

Small  microcalcifications  (and  vessel  calcifications)  in  breast  tissue  can  be  detected  by  USAE 
imaging  method.  In  these  images  microcalcifications  appear  with  high  contrast  with  respect  to  the 
soft  tissue,  hence  can  be  delineated  from  tissue  inhomogeneities. 
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2.f.l)  So  what? 


Our  experimental  results  support  our  hypothesis  that  ultrasound-stimulated  acoustic  emission 
method  (or  vibro-acoustography)  is  capable  of  detecting  small  calcifications  in  human  beast  tissue 
samples.  Improvement  of  this  technique  such  that  it  can  be  reliably  used  in  vivo  would  be  a 
significant  contribution  to  the  field,  especially  where  the  conventional  mammography  could  not  be 
used  for  detection  of  breast  microcalcification. 
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APPENDIX  A  —  Figures  and  Legends 


Figure  1  High-resolution  x-ray  image  of  breast  tissue  samples  collected  from  a  patient.  The 
piece  at  the  top  right  comer  has  some  microcalcifications,  indicated  by  the  arrow.  This  piece 
is  selected  for  scanning. 


11 


Figure  2  Photograph  of  the  tissue  sample  selected  from  the  collection  shown  in  the  previous 
figure.  Microcalcifications  are  located  at  the  top-left  corner  of  this  sample  as  indicated  by  the 
arrow.  Two  sutures  are  knotted  at  two  location  on  the  tissue  as  indicators  and  landmarks  that 
will  show  in  the  acoustic  image.  The  tissue  is  glued  (using  three  small  droplets  at  locations  far 
from  the  calcified  region)  to  a  thin  latex  sheet  attached  to  the  scanning  bracket  for  scanning. 
This  sheet  is  acoustically  transparent. 
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Figure  3  Tissue  mammography  of  the  selected  sample  shown  in  Figure  2.  This  sample  includes 
four  microcalcification  (<400  microns)  at  the  top  left  corner.  The  brightness  of  the  image  has 
been  reduced  to  make  the  microcalcifications  stand  out  from  the  high  density  tissue  at  the  top 
left  comer  of  the  sample.  As  a  result  most  of  the  tissue  is  not  visible  in  this  x-ray,  except  for  the 
calcified  region  and  some  bright  irregular  line  regions  (containing  no  calcification)  at  the  center. 


13 


Figure  4  USAE  image  of  the  sample  shown  in  Figures  2  and  3.  This  image  clearly  shows  four 
microcalcifications  (as  bright  spots)  at  the  top  left  comer  of  the  tissue  sample.  The  two  sutures 
are  also  indicated  in  the  image.  The  location  of  microcalcifications  shown  here  matches  well 
with  those  shown  in  the  x-ray  mammography. 
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Figure  5  High-resolution  x-ray  of  breast  tissue  samples  from  a  patient.  One  of  the  samples 
contains  a  calcified  vessel.  This  sample  is  selected  for  scanning. 
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Figure  6  X-ray  mammogram  of  the  selected  tissue  sample  indicated  in  the  previous  figure.  The 
calcified  vessel  (bright  line)  is  about  one  millimeter  in  diameter.  This  sample  is  glued  to  the 
latex  sheet  of  the  scanning  bracket. 
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Figure  7  USAE  image  of  the  tissue  sample  shown  in  the  previous  figure.  The  calcified  vessels 
shown  as  a  bright  line  with  high  contrast  with  respect  to  the  background  tissue.  Also  seen  in 
this  image  is  the  suture  that  is  used  for  identification.  Note  that  the  calcification  is  displayed 
with  much  higher  contrast  than  tissue  inhomogeneities. 
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Figure  8  This  x-ray  image  shows  a  12x12  mm  portion  of  the  tissue  sample  shown  in  Figure  5, 
which  is  cut  out  for  histology  examination.  The  thickness  of  the  tissue  is  about  3  mm.  The 
calcified  vessel  is  seen  as  a  bright  line. 


Figure  9  This  image  shows  a  histology  slide  from  the  center  portion  of  the  tissue  piece.  The 
slice  is  cut  tangential  to  the  vessel.  This  magnified  slide  shows  parts  of  the  vessel  wall  including 
an  irregular  calcified  layer  on  the  wall.  The  diameter  of  the  vessel  is  about  1  mm  and  the 
thickness  of  the  calcified  layer  is  less  than  100  microns.  The  layer  covers  portions  of  the  wall. 
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ABSTRACT  We  describe  theoretical  principles  of  an  imag¬ 
ing  modality  that  uses  the  acoustic  response  of  an  object  to  a 
highly  localized  dynamic  radiation  force  of  an  ultrasound  Field. 
In  this  method,  named  ultrasound-stimulated  vibro-acoustogra¬ 
phy  (USVA),  ultrasound  is  used  to  exert  a  low-frequency  (in  kHz 
range)  force  on  the  object.  In  response,  a  portion  of  the  object 
vibrates  sinusoidally  in  a  pattern  determined  by  its  viscoelastic 
properties.  The  acoustic  emission  field  resulting  from  object 
vibration  is  detected  and  used  to  form  an  image  that  represents 
both  the  ultrasonic  and  low-frequency  (kHz  range)  mechanical 
characteristics  of  the  object.  We  report  the  relation  between  the 
emitted  acoustic  Field  and  the  incident  ultrasonic  pressure  field 
in  terms  of  object  parameters.  Also,  we  present  the  point-spread 
function  of  the  imaging  system.  The  experimental  images  in  this 
report  have  a  resolution  of  about  700  fx.m,  high  contrast,  and  high 
signal-to-noise  ratio.  USVA  is  sensitive  enough  to  detect  object 
motions  on  the  order  of  nanometers.  Possible  applications  in¬ 
clude  medical  imaging  and  material  evaluation. 


The  study  of  objects  in  terms  of  their  mechanical  response  to 
external  forces  is  of  considerable  interest  in  material  science 
and  medical  diagnosis.  Elastic  constants  are  closely  connected 
to  the  thermodynamic  properties  of  materials  and  can  be 
related  to  a  wide  range  of  physical  parameters.  Elastic  con¬ 
stants  can  be  determined  by  measuring  deformation  in  re¬ 
sponse  to  an  applied  force.  Although  a  static  force  can  be  used 
for  this  purpose,  using  a  dynamic  force  is  preferred  if  one  is 
interested  in  measuring  the  dynamic  characteristics  of  the 
material  (1). 

Changes  in  elasticity  of  soft  tissues  are  often  related  to 
pathology.  Palpation  is  a  traditional  example  of  estimating 
mechanical  parameters  for  tissue  characterization,  where  a 
static  force  is  applied  and  a  crude  estimation  of  the  tissue 
elasticity  is  obtained  through  the  sense  of  touch.  In  palpation, 
force  is  exerted  on  the  body  surface,  and  the  result  is  an 
accumulative  response  of  all  the  tissues  below.  Physicians  can 
sense  abnormalities  if  the  response  to  palpation  of  the  suspi¬ 
cious  tissue  is  sufficiently  different  from  that  of  normal  tissue. 
However,  if  the  abnormality  lies  deep  in  the  body,  or  if  it  is  too 
small  to  be  resolved  by  touch,  then  the  palpation  method  fails. 

Elasticity  imaging,  a  subject  extensively  investigated  in  re¬ 
cent  years,  is  a  quantitative  method  that  measures  the  me¬ 
chanical  properties  of  tissue.  The  general  approach  is  to 
measure  tissue  motion  caused  by  an  external  (or,  in  some 
methods,  internal)  force/displacement  and  use  it  to  recon¬ 
struct  the  elastic  parameters  of  the  tissue.  Some  investigators 
have  used  static  force  to  compress  the  tissue  and  measured  the 
resulting  strain  by  ultrasound  (2, 3).  Others  have  used  external 
mechanical  vibrators  to  vibrate  the  tissue  and  detected  the 
resulting  displacement  in  tissue  by  Doppler  ultrasound  (4-7). 
For  a  review  of  elasticity  imaging  methods,  refer  to  ref.  8.  A 
recently  developed  method  uses  an  actuator  to  vibrate  the  body 
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surface  and  then  measures  the  strain  waves  with  phase- 
sensitive  MRI  (9). 

Most  of  the  elasticity  imaging  methods  are  based  on  an 
external  source  of  force  resulting  in  a  spatially  wide  stress-field 
distribution.  This  requires  the  stress  field  to  pass  through  the 
superficial  portion  of  an  object  before  reaching  the  interior 
part.  Analysis  of  the  object  response  can  be  complicated 
because  the  stress-field  pattern  changes,  often  unpredictably, 
at  different  depths  before  it  reaches  the  region  of  interest 
within  the  object.  An  alternative  strategy  is  to  apply  a  localized 
stress  directly  in  the  region  of  interest.  One  way  to  accomplish 
this  is  to  use  the  radiation  pressure  of  an  ultrasound  source(s). 
Based  on  this  strategy,  Sugimoto  etal.  (10)  presented  a  method 
to  measure  tissue  hardness  by  using  the  radiation  force  of  a 
single  focused  ultrasound  beam.  In  this  method,  impulsive 
radiation  force  was  used  to  generate  localized  deformation  of 
the  tissue.  Resulting  transient  deformation  was  measured  as  a 
function  of  time  by  an  ultrasound  Doppler  method.  Radiation 
force  has  also  been  used  to  generate  shear  elastic  waves  in 
tissues  (11). 

In  this  paper,  we  describe  the  principles  of  an  imaging 
technique  that  produces  a  map  of  the  mechanical  response  of 
an  object  to  a  force  applied  at  each  point.  The  method  uses 
ultrasound  radiation  force  to  remotely  exert  a  localized  oscil¬ 
lating  stress  field  at  a  desired  frequency  within  (or  on  the 
surface  of)  an  object.  In  response  to  this  force,  a  part  of  the 
object  vibrates.  The  size  of  this  part  and  the  motion  pattern 
depend  on  object  viscoelastic  characteristics.  The  acoustic 
field  resulting  from  object  vibration,  which  we  refer  to  as 
“acoustic  emissions/’t  is  detected  by  a  sensitive  hydrophone 
and  used  to  form  the  image  of  the  object.  This  method  benefits 
from  the  high  spatial  definition  of  ultrasound  radiation  force 
and  high  motion-detection  sensitivity  offered  by  the  hydro¬ 
phone.  We  call  this  technique  ultrasound-stimulated  vibro- 
acoustography  (USVA).  Some  general  aspects  of  this  method, 
including  some  experimental  results,  have  been  outlined  by  the 
authors  in  ref.  12.  Here  we  present  the  theoretical  foundations 
of  USVA. 

METHODS 

Our  aim  is  to  image  an  object  based  on  its  mechanical 
characteristics.  This  is  achieved  by  vibrating  the  object  by 
applying  a  highly  localized  oscillating  force  to  each  point  of  the 
object.  The  localized  force  is  produced  by  modulating  the 
intensity,  and  thereby  the  radiation  force,  of  the  ultrasound  at 
low  frequencies  (normally  in  kHz  range).  The  resulting  sound 
emitted  by  the  object  is  a  function  of  object  mechanical 


Abbreviations:  USVA,  ultrasound-stimulated  vibro-acoustography; 
PSF,  point  spread  function;  CW,  continuous  wave. 

*To  whom  reprint  requests  should  be  addressed  at:  Ultrasound 
Research,  Mayo  Clinic,  e-mail:  fatemi.mostafa@mayo.edu. 
tThe  term  “acoustic  emission”  is  used  to  describe  the  acoustic  field  in 
response  to  a  cyclic  vibration  of  the  object.  Similar  terminology  is  also 
used  in  the  field  of  nondestructive  testing  of  materials  and  in 
optoacoustic  imaging  to  describe  a  different  phenomenon,  usually  the 
acoustic  field  resulting  from  structural  deformation. 
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characteristics  and  the  location  of  the  excitation  point.  The 
image  is  produced  by  mapping  the  amplitude  or  phase  of  this 
sound,  which  is  detected  by  a  sensitive  hydrophone,  vs.  posi¬ 
tion.  Fig.  1  illustrates  this  method.  In  the  following  section,  we 
describe  the  relationship  between  the  USVA  image  and  the 
properties  of  the  object. 

THEORY 

Generation  of  a  Dynamic  Radiation  Force  on  a  Target.  The 

acoustic  radiation  force  is  the  time-average  force  exerted  by  an 
acoustic  field  on  an  object.  This  force  is  an  example  of  a 
universal  phenomenon  in  any  wave  motion  that  introduces 
some  type  of  unidirectional  force  on  absorbing  or  reflecting 
targets  in  the  wave  path.  Radiation  force  is  produced  by  a 
change  in  the  energy  density  of  an  incident  acoustic  field.  For 
a  review  of  this  phenomenon,  refer  to  ref.  13.  Consider  a  plane 
ultrasound  beam  interacting  with  a  planar  object  of  zero 
thickness  and  arbitrary  shape  and  boundary  impedance  that 
scatters  and  absorbs.  The  radiation  force  vector,  F,  arising 
from  this  interaction  has  a  component  in  the  beam  direction 
and  another  transverse  to  it.  The  magnitude  of  this  force  is 
proportional  to  the  average  energy  density  of  the  incident  wave 
(E)  at  the  object,  where  < )  represents  the  time  average  and  S, 
the  projected  area  of  the  object  (14) 

F  =  d  ,5<£),  [1] 

where  d,  is  the  vector  drag  coefficient  with  a  component  in  the 
incident  beam  direction  and  another  transverse  to  it.  The 
coefficient  dr  is  defined  per  unit  incident  energy  density  and 
unit  projected  area.  For  a  planar  object,  d,  is  numerically  equal 
to  the  force  on  the  object.  Physically,  dr  represents  the  scat¬ 
tering  and  absorbing  properties  of  the  object  and  is  given  by 
(14) 

dr  =  pS"1  (nfl  +  n,  -  ;  y  cos <xsdS)  +  qS~l  f  y  sin  asdS ,  [2] 

where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and 
normal  to  it,  respectively.  The  quantities  Ua  and  ns  are  the 
total  absorbed  and  scattered  powers,  respectively,  and  y  is  the 
scattered  intensity,  all  expressed  per  unit  incident  intensity. 
Also,  a5  is  the  angle  between  the  incident  and  the  scattered 
intensity,  and  dS  is  the  area  element.  The  drag  coefficient  can 
also  be  interpreted  as  the  ratio  of  the  radiation  force  magni¬ 
tude  on  a  given  object  to  the  corresponding  value  if  the  object 
were  replaced  by  a  totally  absorbing  object  of  similar  size.  This 
is  because  |dr|  =  1  for  a  totally  absorbing  object.  This  coeffi¬ 
cient  can  be  determined  for  objects  of  different  shapes  and 
sizes.  For  simplicity,  we  assume  a  planar  object  normal  to  the 
beam  axis.  In  this  case,  the  transverse  component  vanishes, 
thus  the  drag  coefficient  (force)  will  have  only  a  component 
normal  to  the  target  surface,  which  we  denote  by  scalar  dr  (F). 
Values  of  dr  for  spheres,  in  terms  of  the  diameter  and  the 
wavelength,  are  given  in  ref.  14. 

To  produce  a  dynamic  radiation  force,  one  can  use  an 
amplitude-modulated  beam  (15).  Consider  an  amplitude- 
modulated  incident  (ultrasonic)  pressure  field,  p(t),  as 

p(t)  =  PWo  cos(Aw//2)  coscoo t,  [3] 


Image 


m —  Acoustic 
Emission 


Fig.  1.  Principle  of  ultrasound-stimulated  vibro-acoustography. 


where  PWo,  Aw/2,  and  w0  are  the  pressure  amplitude,  modu¬ 
lating  frequency,  and  center  frequency,  respectively.  In  our 
analysis  and  experiments,  we  assume  that  the  condition  Aw 
wo  holds.  In  such  a  case,  the  energy  density  of  the  incident  field 
has  slow  variations  in  time.  To  discriminate  the  slow  time 
variations  of  a  function,  let  us  define  the  short-term  time  average 
of  an  arbitrary  function  £(/)  over  the  interval  of  T  seconds  at 
time  instance  /,  as  (((t))r  =  I/T’  /I-r/2  £(TWT'  which  is  a 
function  of/.  The  long-term  time  average  (or  simply  the  time 
average)  is  obtained  by  setting  T  —>  00 .  To  compute  the 
short-term  time  average  of  the  acoustic-energy  density  rele¬ 
vant  to  field  variations  at  Aw/2,  we  choose  T  longer  than  the 
ultrasound  wave  period  but  much  shorter  than  the  modulation 
period,  that  is  27r/w0  «  T  <§c  4tt/Aw.  Under  this  condition, 
the  short-term  time  average  ofp2(/)  is  (p2(t))T  =  0  + 

cos  Aw/).  The  energy  density  is  given  by p2(t)/ pc 2,  where  p  and 
c  are  the  density  and  propagation  speed  in  the  medium  (16). 
We  are  interested  in  the  time-varying  component  of  the 
short-term  time  average  of  the  energy  density.  Denoting  tins 
component  by  eA<t)(t),  we  can  write:  eAo}(t)  =  (Pm/Apc  ) 
cosAw/.  This  component  of  the  energy  density  produces  a 
time-varying  radiation  force  on  the  target  (Eq.  1)  at  frequency 
Aw.  The  amplitude  of  this  force,  FAw,  is 

Fsto  —  Pit,,  Sdr/4pc2.  [4] 

This  equation  states  that  the  time-varying  force  amplitude  is 
proportional  to  the  square  of  incident  ultrasound  pressure,  or 
equivalently,  to  the  incident  power.  If  the  object  moves  in 
response  to  this  force,  then  the  high-frequency  ultrasound 
energy  would  convert  to  low-frequency  mechanical  energy.^ 

Acoustic  Emission  from  a  Target  Caused  by  a  Dynamic 
Force.  The  radiation  force  F Aoi  vibrates  the  target  object  at 
frequency  Aw.  Object  vibration  results  in  an  acoustic  field  in 
the  medium  (acoustic  emission).  This  field  is  related  to  object 
shape,  size,  and  viscoelastic  properties.  To  present  a  conclusive 
analysis  of  this  relationship,  we  have  to  assume  an  object  with 
specific  characteristics.  Here  we  assume  that  the  vibrating 
object  has  a  circular  cross-section  of  radius  b  and  uniformly 
vibrates  back  and  forth  like  a  piston.  This  choice  allows  us  to 
illustrate  the  concept  in  a  simple  form.  We  also  consider  an 
area  S  <  7 rb2  of  the  piston  surface  to  be  projected  normally  by 
the  beam.  Similar  solutions  can  be  carried  out  for  other 
objects.^ 

The  steady-state  normal  velocity  amplitude  of  a  piston,  UAo)7 
caused  by  a  harmonic  force  FAo>  at  frequency  Aw,  can  be 
described  in  terms  of  the  mechanical  impedance  ZAo} , 

UAu)  ^  FAJZAo),  [5] 

where  ZAa)  -  Z'm  +  Zr  is  comprised  of  the  mechanical 
impedance  of  the  object  in  vacuum  ZJn,  and  the  radiation 
impedance  of  the  object  Zr,  all  defined  at  Aw.  Modeling  the 
object  as  a  mass-spring  system,  Z'm  can  be  written  in  terms  of 
Aw  as  (16,  17) 

Z'm  =  R'm  -/(mAo)  -  K’/ Aw),  [6] 

where  m,  R'm ,  and K!  are  the  mass,  mechanical  resistance,  and 
spring  constants  of  the  object,  respectively.  The  radiation 
impedance  of  the  piston  can  be  written  (17)  as  Zr  —  tt b2(Rr  — 
jXr)  [7] ,  where Rr  =  pc[  1  -  (c/Awb)/i(c/2Aw6)]J  [8],andZr 
=  (4pc  /  7r)  JJ/2  sin  [(2b^/c)  cos  a]  sin2  a  da  [9],  where /i(-) 
is  the  first-order  Bessel  function  of  the  first  kind.  In  many 
applications  of  our  interest,  the  wavelength  is  much  greater 
than  the  object  size,  hence,  (Aw /c)b  — »  0.  In  such  cases  Zr 
assumes  a  simpler  form  as  Zr  =  7r63pAw(6Aw/2c  —  ;8/3ir) 


tThe  theory  can  be  extended  to  include  arbitrary  vibrating-part  shapes 
and  nonuniform  displacement  of  the  object.  Nonuniform  displace¬ 
ment  would  be  an  important  issue  when  the  vibration  wavelength  in 
the  object  material  is  smaller  than  2b. 
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[10].  The  mechanical  impedance  of  the  piston  object  can  now 
be  written  as 

{R'„,  +  t Tb2Rr)  -  /(/?/ A (u  —  /C'/Aoj  +  'nb2Xr) 

=  (/?,'„  4-  7rp/?4Aw:/2f )  -  -  K' /Ao) 

+  8p/r’A<o/3).  /;(Aco/c)  — >  0.  [11] 

Once  we  calculate  £/*<„.  we  can  calculate  the  pressure  field 
it  produces  in  the  medium.  We  assume  that  the  acoustic 
emission  signal  propagates  in  a  free  and  homogenous  medium. 
The  farficld  acoustic  pressure  caused  by  a  piston  source  of 
radius  b  set  in  a  planar  boundary  of  infinite  extent  is  given  by 
(17), 

2/,f 0Ao>/c)  sin  d]  cos 
{b\co/c)  sin  d  cos  d  +  j3/* 

X  (2 [12] 

where  /  is  the  distance  from  the  observation  point  to  the  center 
of  the  piston,  -6  is  the  angle  between  this  line  and  the  piston 
axis,  and  (3/?  is  the  specific  acoustic  admittance  of  the  boundary 
surfaced  The  factor  of  two  comes  from  the  presence  of  the 
boundary  wall,  which  would  be  replaced  by  unity  if  the 
boundary  wall  were  not  present  (16).  The  acoustic  emission 
field  resulting  from  object  vibration  can  be  written  in  terms  of 
the  incident  ultrasound  pressure  by  combining  Eqs.  4,  5,  and 
12,  as 

27,[(M<t>/c)  sin  fl]  cos  ~j] 
(b&w/c)  sin  i3  *  cos  d  4-  /3/J  J 

X  {1  /[(/C  +  77 b2Rr)  -;(/hAw  -  K'  / Aw  +  irh2Xr)^(2id)2)P\fdr. 

[13] 


f  Aa> 


cxp(/AaV/c) 


Aid 


Pi.,  =  J±(op 


cxp(/Aa>//c) 
4  t tI 


For  wavelengths  long  compared  to  the  object  size,  i.e.,  when 
b&oj/c  0,  the  term  in  the  first  brace  approaches  a  constant, 
hence  we  may  consider  the  contents  of  the  first  brace  to  be  an 
object-independent  function  (the  specific  acoustic  admittance 
Pb  relates  to  the  surrounding  boundary  surface).  Under  these 
conditions,  the  first  brace  in  the  above  equation  represents  the 
effect  of  the  medium  on  the  acoustic  emission  field,  which  we 
may  call  the  medium  transfer  function  ,  and  denote  it  by 


Aw 

=  j  ~^2~  X 


exp(/Aw//c) 
477 / 


U^bAw/c)  sin  0] 
(bAu)/c)  sin  0 


X 


cos 

cos  0  +  (3b. 


[14] 


The  second  brace  in  Eq.  13  is  l/ZAa),  or  the  mechanical 
admittance  of  the  object  at  the  frequency  of  the  acoustic 
emission  (Aco),  and  we  denote  it  by  Ya<u*  It  is  convenient  to 
combine  this  term  with  the  next  term  (2irb2)  in  Eq.  13,  as 
Gao)  =  2irb2YtuJ)  -  2irb2/Z^w  which  is  the  total  acoustic 
outflow  by  the  object  per  unit  force  (acoustic  outflow  is  the 
volume  of  the  medium  (e.g.,  the  fluid)  in  front  of  the  object 
surface  that  is  displaced  per  unit  time  because  of  object 
vibration.).  Function  represents  the  object  characteristics 
at  the  acoustic  frequency.  We  may  thus  rewrite  Eq.  13  in  a 
more  compact  form  as 

Pa.  =  H^JDQ^iSd,.  [15] 


Eq.  15  indicates  that  the  acoustic  emission  pressure  is  propor¬ 
tional  to:  (/)  the  square  of  ultrasound  pressure  P Wo;  (ii)  the 
ultrasound  characteristics  of  the  object,  dr ,  in  the  projected 


area  S;  (iii)  the  acoustic  outflow  by  this  object,  repre¬ 
senting  the  object  size  b  and  its  mechanical  admittance  at  the 
acoustic  frequency,  Y±(0;  and  (/V)  the  transfer  function  of  the 
medium  at  the  acoustic  frequency,  H±to(l).  The  above  equation 
illustrates  the  basic  nonlinear  relationship  between  the  ultra¬ 
sound  and  acoustic  emission  pressure  amplitudes.  Note  that 
neither  the  medium  nor  the  object  needs  to  be  nonlinear  for 
this  relationship  to  hold.  It  is  interesting  to  note  that  the 
projection  area  S  and  the  vibrating  area  rrb2  play  different 
roles.  The  projection  area  determines  the  extent  of  the  force 
applied  to  the  object  (Eq.  4).  The  vibrating  area,  however, 
influences  the  total  acoustic  outflow  in  the  medium  caused  by 
object  vibration.  The  mechanism  of  object  vibration  is  some¬ 
what  analogous  to  that  of  a  loudspeaker,  where  the  electro¬ 
motive  force  is  exerted  at  a  small  area  of  the  membrane 
(usually  at  the  center),  causing  the  entire  membrane  surface  to 
vibrate.  In  our  method,  the  size  of  the  vibrating  area  depends 
on  the  object  structure.  For  a  free  suspended  point  object, 
smaller  than  the  beam  cross-section,  the  vibrating  area  would 
be  the  same  as  the  projected  area.  For  a  large  stiff  plate, 
however,  the  vibrating  area  could  be  much  larger  than  the 
projected  area  (similar  to  a  loudspeaker).  In  some  cases,  it  is 
more  convenient  to  write  the  acoustic  emission  field  in  terms 
of  the  applied  force  F Referring  to  Eq.  4,  wc  can  rewrite  Eq. 
15  as 


Pi,  =  4pc2HxM)QiJi«-  t16l 

Again  in  analogy  to  a  loudspeaker,  F^,  Gao>,  and  H\W{1) 
represent  the  electromotive  force,  dynamic  characteristics  of 
the  membrane,  and  propagation  medium  transfer  function. 

Beam  Forming.  To  probe  an  object  with  the  dynamic 
radiation  force  at  high  spatial  resolution,  it  is  ideal  to  confine 
the  dynamic  stress  field  to  a  very  small  region  in  three- 
dimensional  space.  We  may  define  the  resolution  cell  of  the 
system  as  the  volume  within  which  the  amplitude  of  the 
modulated  field  is  high  enough  to  produce  a  stress  field  on  a 
target.  The  purpose  of  beam  forming  is  to  produce  a  resolution 
cell  as  small  as  possible.  An  amplitude  modulated  single- 
focused  beam  can  provide  a  resolution  cell  that  is  small  in 
diameter  but  long  in  depth  direction.  A  superior  strategy  that 
can  achieve  a  small  resolution  cell  in  all  dimensions  is  to  use 
two  unmodulated  focused  beams  at  slightly  different  frequen¬ 
cies  and  allow  them  to  cross  each  other  at  their  focal  regions. 
This  is  accomplished  by  projecting  two  coaxial  confocal  con¬ 
tinuous-wave  (CW)  ultrasound  beams  on  the  object.  An  am¬ 
plitude-modulated  field  is  produced  only  at  the  interference 
region  of  the  two  unmodulated  beams  around  their  focal  areas, 
resulting  in  a  small  resolution  cell.  For  this  purpose,  elements 
of  a  two-element  spherically  focused  annular  array  (consisting 
of  a  central  disc  with  radius  a  \  and  an  outer  ring  with  the  inner 
radius  of  a  2  and  outer  radius  of  af)  are  excited  by  separate  CW 
signals  at  frequencies  =  a)o  —  Acu/2  and  w2  =  m  +  Ao>/2. 
We  assume  that  the  beams  are  propagating  in  a  lossless 
medium,  in  the  +z  direction  of  a  Cartesian  coordinate  system 
( x ,  y ,  z),  with  the  joint  focal  point  at  z  =  0.  The  resultant 
pressure  field  on  the  z  =  0  plane  may  be  written  as 

p(t)  =  P\(r)  cos (o)it  +  i/>i(r))  +  Pz(f)  cos (a)2t  +  [17] 

where  r  =  VF+p  is  the  radial  distance.  The  amplitude 
functions  are  (16,  18) 

Pfr)  =  pcUoi  (uaJ/A^o)  jinc(ra [18] 


and 


P2(f)  =  pct/02(Wo)  \p\ jinc(ra2/A2Zo)  -  a22 jmc(w2/A2z0)],  [19] 


§The  specific  acoustic  admittance  is  (3b  =  pc/Zs ,  where  Zb ,  the 
acoustic  impedance  of  the  boundary,  represents  the  ratio  between  the 
pressure  and  normal  fluid  velocity  at  a  point  on  the  surface. 


where  A/  =  2Tv/o>iy  i  =  1,  2,  is  the  ultrasound  wavelength,  Uoi 
is  the  particle  velocity  amplitude  at  the  i-th  transducer  element 
surface,  and  jinc  ( X)  =  /i(27iAT)/ ttX.  The  phase  functions, 
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^.(r)  =  —  7rr2/A/Zn.  for  /  =  1,  2,  are  conveniently  set  to  be  zero 
at  the  origin. 

Now.  we  define  a  unit  point  target  at  position  (x0,yo)  on  the 
focal  plane  with  a  drag  coefficient  distribution  as 

</,.( .v.  y)  =  8{x  -  a'o,  y  -  y(>),  [20] 

such  that  d,(x,  y)dxdy  is  unity  at  (x0,  yo)  and  zero  elsewhere. 
This  equation  is  merely  used  as  a  mathematical  model  because 
dr  is  physically  finite.  In  this  case,  the  projected  area  can  be 
considered  to  be  S  =  dxdy.  We  replace  d,S  in  Eq.  1  with  dr(x , 
y)dxdy  and  follow  the  steps  similar  to  those  outlined  in  Eqs.  3 
and  4  for  the  pressure  field  expressed  by  Eq.  17,  then  the 
complex  amplitude  of  the  normal  component  of  the  force  on 
the  unit  point  target  can  be  found  as 

F*Jfco.yo)  =  pUoiUv  (™q/4A,z0)  }mc(r0al/KiZfi)[(Tral/K2Zo) jinc(r0tf2/A2Zo) 

-  {jta'}! A2z0)  jinc(/-0fl2/A2Zo)]cxp(-;^Aw/2czo),  [21] 

where  the  arguments xq  and y0  are  added  to  denote  the  position 
of  the  point  target  and  /*o  =  \/xq~+  yn-  Eq.  21  describes  the 
spatial  distribution  of  the  force  (the  stress  field).  This  equation 
shows  that  the  stress  field  is  confined  to  the  regions  near  the 
beam  axis  (r«  =  0)  and  decays  as  the  radial  distance  r0 
increases.  The  lateral  extent  of  the  stress  field,  and  hence  the 
resolution  cell  diameter,  would  be  smaller  at  higher  ultrasound 
frequencies  (smaller  A,  and  A2).  One  can  calculate  the  total 
force  on  an  arbitrary  object  by  integrating  the  force  over  the 
projected  area.  The  axial  extent  of  the  resolution  cell  (depth 
resolution  or  the  depth  of  field)  can  be  determined  by  calcu¬ 
lating  the  force  FAu,  as  a  function  of  the  depth  variable,  in  a 
fashion  as  outlined  in  Eqs.  17  to  21.  For  conciseness,  we  will 
present  only  the  measured  values  for  the  depth  resolution  in 
Results. 

Loss  in  the  propagation  path  would  attenuate  both  ultra¬ 
sound  beams,  thus  less  radiation  force  would  be  generated  by 
the  remaining  ultrasound  energy.  In  the  case  of  soft  tissues,  the 
force  attenuation  factor  is,4(z0)  =  exp[az0(wi  +  w2)],  where 
a  is  the  attenuation  coefficient  of  the  tissue.  Energy  loss  in  the 
medium  would  also  result  in  generation  of  a  separate  radiation 
stress  on  the  medium  along  the  ultrasound  paths.  However, 
because  the  two  beams  propagate  along  separate  paths  in  the 
CW  form,  they  exert  mainly  steady  radiation  stresses  to  the 
medium,  which  does  not  cause  object  or  medium  vibrations. 
Dynamic  radiation  force  is  produced  only  in  the  interference 
region  around  the  focal  area,  which  is  another  advantage  of 
using  two  unmodulated  beams  over  a  modulated  single  beam. 

Image  Formation.  To  produce  an  image,  we  scan  the  object 
in  a  plane  and  record  the  complex  amplitude  of  the  acoustic 
emission,  PAa„  at  different  positions.  In  this  process,  we  keep 
Aw  feed.  For  transverse  view  images,  the  scan  plane  is  the  focal 
plane  (x-y).  Alternatively,  for  the  parallel  view  the  scan  plane 
is  the  x-z  plane.  In  the  conventional  ultrasound  imaging 
context,  these  two  views  are  called  the  C-scan  and  B-scan, 
respectively.  Our  main  focus  here  is  the  transverse  view 
imaging.  In  this  case,  the  acoustic  emission  data  obtained  by 
vibrating  the  object  at  point  (. x ,  y)  are  assigned  to  the  corre¬ 
sponding  point  (x,  y)  in  the  image. 

Before  defining  the  image,  we  need  to  define  the  function 
that  represents  the  object.  Referring  to  Eq.  15,  the  terms  that 
are  object  dependent  are  the  drag  coefficient  dr  and  the 
function  Qao>  (assuming  that  H&JJ)  is  object  independent). 
The  object  function  g(x,  y)  is  thus  defined  as  the  spatial 
distribution  of  these  terms, 

g{x,  y)  =  Gao/x,  y)dr(x ,  y).  [22] 

Variables x  andy  are  added  to  denote  the  dependency  of  dr  and 
<2a<o  on  position.  In  particular,  Q±Jx,  y)  implies  the  total 


acoustic  outflow  by  the  object  when  unit  force  is  applied  at 

point  (x,y).  . 

Commonly,  an  imaging  system  is  studied  through  its  point- 
spread  function  (PSF),  which  is  defined  as  the  image  of  a  point 
object.  To  determine  the  PSF  of  our  system,  we  consider  a  unit 
point  target  at  the  origin  with  unit  mechanical  response, 
QAw(x,y)  =  1.  Hence,  referring  to  Eqs.  20  and  22,  we  can  write 
g(x,y)  -  S(x,  y).  To  obtain  the  PSF,  we  move  this  point  object 
to  every  possible  position  (xo,yo)  on  thez  =  0  plane  and  form 
the  image  using  the  resulting  acoustic  emission  field,  P aw(-'o* 
yo).  Because  xo  andyo  are  now  being  treated  as  variables,  we 
may  replace  them  by  variables  a'  andy,  respectively.  We  define 
the  normalized  PSF  of  the  coherent  imaging  system  as  the 
complex  function 

/7(x,y)=PUx,y)/PAJ  0,0).  [23] 

Division  by  Pa^O,  0)  cancels  the  constant  multipliers.  Refer¬ 
ring  to  Eqs.  16  and  21,  we  can  write 

h(x,y)  =  (a\  -  fl22)-1  jinc(ra1/A12r0)[^l  jinc(rfl2/A2z0) 

—  a’2  jinc(ra2/A2z0)l  exp(— 7>2Aoj/2cz0).  [24] 

This  equation  illustrates  that  the  system  PSF  is  a  circularly 
symmetric  function  with  the  peak  at  the  origin  and  decaying 
amplitude  with  increasing  the  radial  distance  r .  Amplitude 
decays  faster  for  higher  ultrasound  frequency.  This  function 
will  be  discussed  further  in  the  next  section. 

EXPERIMENTS 

The  experimental  setup  is  shown  in  Fig.  2.  The  confocal 
transducer  is  constructed  by  using  a  spherical  piezoelectric  cap. 
The  two  elements  are  constructed  by  dividing  the  back  elec¬ 
trode  of  the  piezoelectric  wafer  into  a  central  disc  and  the  outer 
ring,  such  that  the  elements  have  identical  beam  axes  and  focal 
lengths.  Radii  of  the  elements  are  ax  =  14.8  mm,  a2  =  22.5 
mm,  and  a2  =  16.8  mm,  and  the  focal  distance  is  70  mm. 
Transducer  elements  were  driven  by  two  stable  radio  fre¬ 
quency  synthesizers  [Hewlett-Packard  33 120 A  and  Analogic 
2045  (Peabody,  MA)]  at  frequencies  of/o  -  A// 2  and/0  + 
A//2,  where  fQ  =  3  MHz,  and  the  value  of  A/  =  Aw/2 tt  is 
stated  separately  for  each  experiment.  The  object  was  placed 
at  the  focal  plane  of  the  ultrasound  beams  in  a  water  tank. 
Sound  produced  by  the  object  vibration  was  detected  by  an 
audio  hydrophone  (International  Transducer,  Santa  Barbara, 
CA,  model  680,  sensitivity  -154  dB  re  lV/jLiPa)  placed  within 
the  water  tank.  The  received  signal  was  filtered  and  amplified 
by  a  programmable  filter  (Stanford  Research  Sunnyvale,  CA, 
SR650)  to  reject  noise,  then  digitized  by  a  12  bits/sample 
digitizer  (Hewlett-Packard  E1429A)  at  a  rate  sufficiently 
higher  than  the  Nyquist  rate  for  the  particular  A/  used.  Data 
were  recorded  on  a  computer  disc.  For  coherent  imaging, 
which  requires  the  phase  information,  the  reference  signal 
(i.e.,  cos  A  cot)  was  obtained  by  electronic  downmixing  of  the 
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Fig.  2.  Ultrasound-stimulated  vibro-acoustography  system.  The 
confocal  ultrasound  annular  array  transducer  with  two  elements  is 
shown  on  the  left. 
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two  driving  signals  and  was  recorded  along  with  the  hydro¬ 
phone  signal.  The  relative  phase  of  the  acoustic  emission  data 
was  then  calculated  at  each  point  by  discrete  Hilbert  transform. 
We  conducted  two  experiments.  The  first  experiment  is  de¬ 
signed  to  verify  the  relationship  between  the  acoustic  emission 
pressure  and  the  ultrasound  pressure  (Eq.  15),  and  the  second 
is  designed  to  experimentally  measure  the  PSF  stated  in  Eq.  24. 
These  experiments  are  presented  to  prove  the  principles  of  the 
method.  Further  experiments,  illustrating  applications  of  the 
method  for  evaluation  of  object  mechanical  properties  and 
tissue  imaging,  are  presented  in  ref.  12. 

RESULTS 

Acoustic  Emission  vs.  Ultrasound  Pressure.  Eq.  15  states 
that  the  acoustic  emission  field  amplitude  is  linearly  propor¬ 
tional  to  the  square  of  the  incident  ultrasound  pressure,  or 
equivalently,  to  the  incident  power.  To  test  this  hypothesis,  a 
calibrated  1-mm-diameter  ultrasound  needle  hydrophone, 
with  its  tip  facing  the  ultrasound  beam,  was  used  to  measure 
the  ultrasound  field  at  the  focal  point.  The  tip  of  the  hydro¬ 
phone  also  served  as  an  object  to  generate  the  acoustic 
emission.  Here.  A f  was  set  at  40  kHz.  and  the  radial  distance 
from  the  tip  of  the  needle  hydrophone  to  the  audio  hydro¬ 
phone  was  about  65  mm.  The  result  is  shown  in  Fig.  3.  The 
slope  of  the  acoustic  emission  intensity  vs.  ultrasonic  intensity 
indicates  that  the  intensity  of  the  acoustic  emission  field  is 
proportional  to  the  square  of  the  ultrasound  intensity,  or 
equivalently,  the  acoustic  emission  pressure  amplitude  is  lin¬ 
early  proportional  to  the  ultrasound  power,  as  predicted  by  Eq. 
15.  In  another  experiment,  a  450-p.m-diameter  glass  bead  was 
used  as  a  point  object.  In  this  case,  A/ was  set  at  either  7  or  40 
kHz.  The  radial  distance  from  the  glass  bead  to  the  audio 
hydrophone  was  about  50  mm.  Again,  the  data  indicate  a 
quadratic  relationship  between  the  acoustic  emission  and 
ultrasonic  intensities.  These  glass-bead  data  also  show  that 
increasing  the  frequency  increases  the  acoustic  intensity.  This 
can  be  better  understood  by  investigating  the  theoretical  model 
presented  in  Eq.  13.  If  we  assume  that  the  mechanical  admit¬ 
tance  of  the  object  is  approximately  constant  at  these  frequen¬ 
cies,  then  the  object  behaves  almost  as  a  point  source,  and  the 
acoustic  pressure  amplitude  is  proportional  to  A (o.  Hence  the 
acoustic  emission  intensities  Ja^  and  1^  at  frequencies  Aou 
and  Aa>2,  respectively,  are  related  by  ~  (Agj2/Awi)2. 

Now,  letting  Acoj  and  'Acjj  correspond  to  7  and  40  kHz, 
respectively,  we  can  write  the  intensity  ratio  as:  [/aw/Ja^Jcib  — 
20  log  (Ao>2/Awi)  =  15  dB.  The  mean  value  of  the  intensity 
ratio  calculated  from  the  glass-bead  data  at  7  and  40  kHz  is  16 
dB,  which  is  in  close  agreement  with  the  theoretical  result. 

PSF  Measurement.  To  demonstrate  the  image-formation 
process  and  support  the  theoretical  derivation  of  the  PSF  (Eq. 
24),  we  evaluated  this  function  experimentally.  For  this  pur¬ 
pose,  we  used  a  380-/utm  diameter  glass  bead  as  a  model  for  a 
point  and  placed  it  on  a  thin  latex  sheet.  The  latex  sheet 
produces  only  a  small  change  in  the  incident  energy,  and  hence 
does  not  produce  significant  radiation  force  or  acoustic  emis¬ 
sion.  The  entire  object  was  placed  in  a  water  tank  and  the  latex 
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Fig.  3.  Acoustic-emission  field  intensity  vs.  the  combined  ultra¬ 
sound  intensity. 
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Fig.  4.  USVA  images  of  a  380-/xm  glass  bead:  (A)  in-phasc,  (B) 
quadrature,  (C)  phase,  and  ( D )  magnitude.  The  phase  in  C  ranges  from 
-77  radians  (black  regions)  to  +77  radians  (white  regions),  and  was 
normalized  to  be  zero  at  the  center  of  the  glass  bead.  (Modified  with 
permission  from  ref.  12,  copyright  1998,  American  Association  for  the 
Advancement  of  Science.) 

sheet  surface  was  scanned  in  a  raster  format  at  0.2-mm 
increments  in  either  direction  at  A/  =  7.3  kHz.  The  amplitude 
and  phase  of  the  acoustic  emission  signal  were  calculated  at 
each  point  relative  to  the  reference  signal  data.  The  phase  was 
normalized  to  the  phase  value  at  the  center  of  the  bead.  The 
resulting  in-phase,  quadrature,  phase,  and  magnitude  images 
are  shown  in  Fig.  4A-D.  Transverse  image  resolution,  defined 
as  the  —  6-dB  width  of  the  bead  image,  is  approximately  700 
l± m  for  the  in-phase  image  in  either  dimension  (refer  to  Fig.  5). 
To  compare  the  experimental  results  with  those  of  the  theory, 
we  calculated  the  profile  of  the  PSF  for  the  transducer 
parameters  used  in  this  experiment  according  to  Eq.  24.  Fig. 
5  shows  the  theoretical  PSF  profile  and  the  glass-bead  profile 
obtained  from  the  experiment  (in-phase  image  profile,  shifted 
to  center  at  zero).  This  figure  shows  excellent  agreement 
between  theory  and  experiment  for  amplitudes  above  20%  of 
the  peak.  The  experimental  data  show  some  background  offset 
about  12%  of  the  peak.  We  believe  that  this  background  offset 
is  caused  at  least  by  the  following  sources:  (i)  acoustic  emission 
by  the  latex  sheet;  (ii)  the  background  acoustic  noise  in  the 
experimental  setup  caused  by  equipment  fans  and  some  struc¬ 
tural  and  building  vibrations;  (in)  nonlinearity  of  water  that 
can  produce  a  nonlinear  mix  of  the  two  beams  even  in  the 
absence  of  the  object;  and  (iv)  streaming  as  a  result  of  energy 
absorption  by  water  (14),  which  in  turn  can  vibrate  the  latex 
sheet  blocking  the  stream.  To  evaluate  the  depth  resolution  (or 
the  slice  thickness),  we  placed  the  glass  bead  on  the  beam  axis 
and  scanned  it  in  the  z-direction  about  the  focal  point.  The 
depth  resolution,  defined  as  the  distance  between  the  points 
where  the  amplitude  of  the  acoustic  emission  field  drops  to  -6 
dB  of  its  peak,  was  9  mm. 

DISCUSSION 

System  Properties.  Spectral  characteristics.  In  general, 
USVA  images  represent  object  characteristics  at  two  ends  of 
the  spectrum:  the  drag  coefficient  at  the  ultrasound  frequency 
and  the  mechanical  admittance  at  the  low  acoustic  frequency 
A/.  The  ultrasound  frequency  is  usually  set  at  a  value  suitable 
to  form  the  beam,  whereas  A f  can  vary  in  a  wide  range 
depending  on  the  application.  If  the  two  beams  are  produced 
by  similar  ultrasound  transducer  elements,  then  the  practical 


Radial  distance  (mm) 

Fig.  5.  The  theoretical  PSF  profile  of  the  USVA  system  according 
to  Eq.  24  and  the  glass-bead  in-phase  image  profile  (Fig.  4 A)  obtained 
from  the  experiment. 
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upper  limit  for  A/ is  about  equal  to  the  transducer  bandwidth. 
The  lower  limit  of  A f  is  zero. 

Sensitivity.  System  sensitivity  in  detecting  very  small  displace¬ 
ments  is  an  important  practical  issue,  especially  when  the  allow¬ 
able  ultrasound  power  is  limited  (for  example,  in  medical  imag¬ 
ine).  Motion  measurement  with  ultrasound  pulse  echo  has  been 
used  previously  to  study  “stiffness"'  of  tissues  (10).  However,  the 
sensitivity  of  ultrasound  pulse  echo  to  motion,  at  common 
medical  ultrasound  frequencies,  is  limited  to  several  micrometers. 
An  advantage  of  USVA  is  its  high  displacement  sensitivity.  Cyclic 
displacement  of  100  nm  at  10  kHz  produces  an  acoustic  intensity 
of  about  3.0  X  If)-3  watts/cm2.  Hydrophones  similar  to  the  one 
used  in  these  experiments  arc  sensitive  to  as  little  as  10  15 
watts /cm2  and  therefore  are  capable  of  delecting  very  small  cyclic 
displacements.  For  instance,  the  hydrophone  detected  an  acoustic 
pressure  of  about  15  X  10“3  Pa  at  a  distance  of  5  cm  from  the  glass 
bead  shown  in  Fig.  4.  Under  assumptions  of  isotropic  vibration, 
this  pressure  would  be  produced  by  a  similar-sized  sphere  vibrat¬ 
ing  with  a  displacement  amplitude  of  about  6  nm.  Sensitivity 
increases  at  higher  frequencies  because  the  acoustic  emission 
pressure  is  proportional  to  frequency  for  constant  mechanical 
admittance  (Eq.  13). 

Comparisons  with  pulse-echo  systems.  Some  contrasting 
features  of  USVA  with  respect  to  the  conventional  ultrasound 
pulse-echo  imaging  (B-modc  and  C-mode)  arc:  (/)  A  pulse- 
echo  image  represents  object  microstructurc  by  displaying  its 
ultrasonic  reflectivity  distribution.  The  acoustic  emission  sig¬ 
nal  in  a  USVA  system  is  proportional  to  the  drag  coefficient 
dr(x,  y),  which  is  a  local  ultrasound  parameter,  and  to  the 
function  y)>  which  represents  the  bulk  response  of  the 

object  at  acoustic  frequency  Aw.  Hence,  a  USVA  image,  in 
general,  represents  both  the  microstructurc  and  macrostruc¬ 
ture  of  the  object,  (ii)  The  echo  signal  in  a  pulse-echo  system 
is  a  linear  function  of  the  incident  ultrasound  pressure  ampli¬ 
tude  and  the  amplitude  reflection  coefficient  of  the  object.  In 
USVA,  however,  the  acoustic  emission  signal  is  proportional 
to  the  ultrasound  power  and  the  power  reflection  coefficient  of 
the  object.  (Hi)  Pulse-echo  systems  are  not  directly  sensitive  to 
medium  absorption.  Absorption  is  indicated  as  relative 
changes  in  the  amplitude  of  the  echoes  resulting  from  the 
scatterers  within  the  medium.  In  a  USVA  system,  the  acoustic 
emission  can  be  produced  directly  as  a  result  of  energy 
absorption  by  the  medium,  even  if  the  medium  is  homogeneous 
(refer  to  Eq.  2).  (iv)  Pulse-echo  systems  are  generally  broad¬ 
band.  The  USVA  method  presented  here  is  basically  a  nar¬ 
rowband  technique,  (v)  Pulse-echo  systems  achieve  high 
depth  resolution  by  transmitting  short  wideband  pulses.  A 
USVA  system  gains  its  depth  resolution  by  tailoring  beam 
geometry  to  limit  the  depth  of  the  region  where  the  two  beams 
interfere.  A  USVA  system  does  not  require  a  wide  bandwidth 
signal  to  achieve  a  high  depth  resolution,  (vi)  The  data  of 
USVA  images  are  acquired  one  point  at  a  time,  which  resem¬ 
bles  the  data  acquisition  in  C-mode  pulse-echo  systems. 
B-mode  pulse- echo  systems,  however,  require  much  less  ac¬ 
quisition  time  because  the  data  are  collected  one  line  at  a  time. 

Applications.  USVA  promises  applications  in  two  general 
areas:  medical  imaging  and  material  evaluation. 

Medical  applications.  USVA  can  be  used  to  image  tissues  and 
evaluate  their  mechanical  characteristics.  To  use  USVA  for  in 
vivo  applications,  one  must  take  into  account  limitations  such  as 
safe  (ultrasound)  power  limit,  tissue  attenuation,  body  noise,  and 
phase  aberration.  The  ultrasound  power  required  to  generate  a 
detectable  acoustic  emission  depends  on  the  object,  acoustic 
noise,  and  receiver  sensitivity.  Experimental  results  shown  in  Fig. 
3  demonstrate  that  ultrasound  intensities  as  low  as  30  mW/cm2 
are  sufficient  to  detect  the  acoustic  emission  from  a  1-mm- 
diameter  object  with  our  hydrophone  in  the  water  tank.  This 
power  value  is  much  smaller  than  the  FDA  limit  for  safe 
diagnostic  ultrasound  applications.  Tissue  attenuation  reduces 
the  ultrasound  intensity  at  the  target,  and  hence  the  acoustic 


emission  [by  the  factor  A(z0)\.  Attenuation  limits  the  usable 
ultrasound  frequency,  and  hence,  lowers  the  resolution.  It  also 
limits  the  signal-to-noise-ratio  (SNR)  because  of  the  loss  in  the 
acoustic  emission  energy  as  a  result  of  ultrasound  attenuation  by 
tissue.  Direct  effect  of  tissue  attenuation  on  the  acoustic  emission 
signal  is  probably  negligible  because  attenuation  of  the  compres- 
sional  waves  at  frequencies  in  the  order  of  a  few  kHz  in  soft  tissues 
is  low.  Sources  of  biological  noise  of  the  human  body  include 
cardiovascular  and  respiratory  systems  and  muscle  movements. 
Body  noise  is  usually  concentrated  below  1  KHz  and  can  be 
filtered  out  if  A f  is  above  this  value.  The  SNR  can  be  improved 
by  increasing  the  time  duration  ot  the  signal  recorded  at  each 
point  (to  increase  the  signal  energy)  and  by  using  very  narrow- 
band  filters  (to  reject  the  noise),  while  keeping  the  ultrasound 
power  within  the  safe  level.  Phase  aberration  in  tissue  can  reduce 
the  sensitivity  of  the  system  by  decreasing  the  effective  ultrasound 
energy  density  at  the  beam  interaction  region.  One  may  use 
known  phase  aberration  correcting  methods  to  reduce  such  an 
effect.  The  practical  value  of  these  methods  for  USVA  remains 
to  be  studied. 

Material  evaluation.  Another  field  in  which  USVA  can  be 
potentially  useful  is  material  characterization,  including  me¬ 
chanical  parameter  evaluation,  imaging,  and  nondestructive 
testing  of  materials.  USVA  can  be  used  for  detection  and 
imaging  flaws  in  materials.  Also,  one  may  use  USVA  to 
evaluate  the  mechanical  frequency  response  of  an  object  at  low 
frequencies.  In  such  case,  we  are  interested  in  determining 
Q£uo(x,y)  vs.  frequency.  We  assume  that  the  object  is  uniform 
within  the  projected  area  S.  Then,  the  total  radiation  force  on 
this  object,  FA<1„  can  be  calculated  by  integrating  FAw(x.y)  over 

S.  Referring  to  Eq.  21,  one  can  show  that  for  Aco  wo,  F Aw 

is  virtually  independent  of  Ao>.  If  //Aa,(/)  is  known  and  nonzero, 
then  the  function  Q^(x,y)  can  be  estimated  using  Eq.  16  as 
Qao,(*,  y)  =  FAw( [x,  y)/F^JT^(l).  In  practice,  PAo)(x,  y)  is 
obtained  by  sweeping  Aco  in  the  range  of  interest  and  recording 
the  resulting  acoustic  emission  (12). 
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Abstract:  In  the  past  few  years  there  has  been  an  increasing  interest  in  using  the  radiation 
force  of  ultrasound  for  evaluating,  characterizing,  and  imaging  biological  tissues.  Of  particular 
interest  are  those  methods  that  measure  the  dynamic  properties  of  tissue  at  low  frequencies.  In 
this  paper  we  present  dynamic  radiation  force  methods  for  probing  tissue  as  a  new  field  and 
discuss  the  inter-relationship  of  several  methods  within  this  field  and  comparing  their  features. 
The  techniques  in  this  field  can  be  categorized  in  three  groups:  transient  methods,  shear  wave 
measurement  methods,  and  a  recently  developed  method  called  vibro-acoustography.  The  last 
method  is  the  focus  of  this  paper.  After  briefly  describing  the  key  concepts  of  the  first  two 
methods,  we  will  present  a  detailed  description  of  vibro-acoustography.  Finally,  we  will  compare 
capabilities  and  limitations  of  these  methods. 

1.  Introduction 

It  is  well  known  that  changes  in  elasticity  of  soft  tissues  are  often  related  to  pathology. 
Traditionally,  physicians  use  palpation  as  a  simple  method  for  estimating  mechanical  properties 
of  tissue.  In  this  method,  a  static  force  is  applied  and  a  crude  estimation  of  tissue  elasticity  is 
obtained  through  the  sense  of  touch.  In  palpation,  the  force  is  applied  on  the  body  surface  and  the 
result  is  a  collective  response  of  all  the  tissues  below.  Clinicians  can  sense  abnormalities  if  the 
response  to  palpation  of  the  suspicious  tissue  is  sufficiently  different  from  that  of  normal  tissue. 
However,  if  the  abnormality  lies  deep  in  the  body,  or  if  it  is  too  small  to  be  resolved  by  touch, 
then  the  palpation  method  fails.  The  dynamic  response  of  soft  tissue  to  a  force  is  also  valuable 
in  medical  diagnosis.  For  instance,  rebound  of  tissue  upon  sudden  release  of  pressure  exerted 
by  the  physician’s  finger  on  the  skin  provides  useful  diagnostic  information  about  the  tissue. 

Quantitative  measurement  of  the  mechanical  properties  of  tissues  and  their  display  in  a  raster 
format  is  the  aim  of  a  class  of  techniques  generally  called  elasticity  imaging  or  elastography. 
The  general  approach  is  to  measure  tissue  motion  caused  by  an  external  (or,  in  some  methods, 
internal)  force  or  displacement  and  use  it  to  reconstruct  the  elastic  parameters  of  the  tissue.  The 
excitation  stress  can  be  either  static  or  dynamic  (vibration).  Dynamic  excitation  is  of  particular 


2 


E-mail:  fatemi.mostafa@mayo.edu 
E-mail:  jfg@mayo.edu 


2 


interest  because  it  provides  more  comprehensive  information  about  tissue  properties  in  a  spectrum 
of  frequencies,  or  alternatively,  the  transient  behavior  of  the  tissue  could  be  deduced  from  the 
measurements.  In  most  elasticity  imaging  methods,  ultrasound  is  used  to  detect  the  motion  or 
displacement  resulting  from  the  applied  stress.  Magnetic  resonance  elastography  is  a  recently 
developed  method  (Muthupillai  et  al  1995)  that  employs  a  mechanical  actuator  to  vibrate  the 
body  surface  and  then  measures  the  strain  waves  with  phase  sensitive  Magnetic  Resonance 
Imaging  (MRI)  machine. 

The  majority  of  elasticity  imaging  methods  are  based  on  an  external  source  of  force.  In  these 
methods  the  object  is  pressed  by  a  known  amount  of  force  or  displacement,  and  the  resulting 
internal  deformations  are  measured  by  means  of  pulse-echo  ultrasound.  Elasticity  of  the  region 
of  interest  is  then  calculated  based  on  the  resulting  deformation  in  relation  to  the  magnitude  of 
applied  force  (or  displacement).  Normally,  the  region  of  interest  rests  deep  in  the  body  and  away 
from  the  source  of  the  force.  The  force  that  is  actually  exerted  on  the  region  of  interest  depends 
on  the  elastic  properties  of  the  tissues  located  between  the  source  and  the  region  of  interest. 
Hence,  the  deformation,  and  the  estimated  elasticity  of  the  region  of  interest,  are  subject  to 
variables  of  the  intervening  tissues. 

An  alternative  strategy  is  to  apply  a  localized  stress  directly  in  the  region  of  interest.  One 
way  to  accomplish  this  is  to  use  the  radiation  pressure  of  ultrasound.  Acoustic  radiation  force 
is  the  time  average  force  exerted  by  an  acoustic  field  on  an  object.  This  force  is  produced  by  a 
change  in  the  energy  density  of  an  incident  acoustic  field  (Chu  and  Apfel  1982),  for  example,  due 
to  absorption  or  reflection.  Several  benefits  may  result  from  using  ultrasound  radiation  force  for 
evaluating  tissue  properties,  including:  (a)  acoustic  (ultrasound)  energy  is  a  noninvasive  means 
of  exerting  force,  (b)  existing  ultrasound  technology  and  devices  can  be  readily  modified  for  this 
purpose,  thus  eliminating  the  need  for  developing  a  new  technology,  (c)  radiation  force  can  be 
generated  remotely  inside  tissue  without  disturbing  its  superficial  layers,  (d)  the  radiation  stress 
field  can  be  highly  localized,  thus  allowing  for  precise  positioning  of  the  excitation  point,  and 
(e)  radiation  force  can  be  produced  in  a  wide  range  of  frequencies  or  temporal  shapes.  These 
features  make  radiation  force  methods  more  attractive  than  other,  mostly  mechanical,  excitation 
methods  used  in  elasticity  imaging. 

Tissue  probing  with  the  radiation  force  of  ultrasound  can  be  accomplished  by  a  variety  of 
techniques,  depending  on  the  excitation  and  detection  methods  used.  As  in  elasticity  imaging 
methods  with  mechanical  excitation,  the  radiation  force  methods  can  use  either  a  static  or  dynamic 
stress. 

It  is  the  opinion  of  the  authors  that  using  the  dynamic  radiation  force  to  remotely  probe  tissue 
has  certain  unique  characteristics  and  capabilities  that  can  open  a  new  chapter  in  the  field  of  tissue 
characterization  and  imaging.  We  believe  that  the  techniques  based  on  this  phenomenon  have 
the  potential  to  provide  a  new  set  of  information  about  tissue  that  has  never  been  available  from 
other  methods.  We  also  believe  that  it  is  important  to  recognize  and  treat  this  set  of  techniques 
as  a  new  field.  Especially,  it  is  insightful  to  set  this  new  field  apart  from  the  conventional 
ultrasound  tissue  characterization/imaging  field.  A  major  difference  being  the  fact  that  using 
the  dynamic  radiation  stress  allows  one  to  analyze  the  object  based  on  its  structural  vibration 
properties  as  opposed  to  its  ultrasonic  parameters. 
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Our  intention  in  writing  this  paper  is  to  present  the  dynamic  radiation  force  methods  for 
probing  tissue  as  a  new  field  by  discussing  the  inter-relationship  of  several  methods  within  this 
field  and  compare  their  features.  To  the  best  of  our  knowledge,  such  a  discussion  has  not  been 
presented  before  in  the  literature. 

The  dynamic  radiation  force  methods  may  be  categorized  as:  (a)  transient  methods,  where  an 
impulsive  radiation  force  is  used  and  the  transient  response  of  the  tissue  is  detected  by  Doppler 
ultrasound  (Sugimoto  et  al  1990);  (b)  shear- wave  methods,  where  an  impulsive  or  oscillating 
radiation  is  applied  to  the  tissue  and  the  resulting  shear  wave  is  detected  by  ultrasound  or  other 
methods  (Andreev  et  al  1997;  Sarvazyan  et  al  1998;  Walker  1999);  (c)  vibro-acoustography,  a 
method  recently  developed  by  the  authors,  where  a  localized  oscillating  radiation  force  is  applied 
to  the  tissue  and  the  acoustic  response  of  the  tissue  is  detected  by  a  hydrophone  or  microphone. 

In  this  paper  we  systematically  discuss  the  features  and  capabilities  that  are  offered  by  the 
dynamic  radiation  force  techniques.  Here,  we  describe  the  general  approach  used  in  the  transient, 
shear  wave,  and  vibro-acoustography  methods.  We  pay  particular  attention  to  the  last  method, 
present  its  theory,  and  discuss  its  features.  Then,  we  compare  the  capabilities  and  limitations  of 
all  three  methods  and  discuss  their  possible  applications. 


2.  Radiation  force  methods 

The  acoustic  radiation  force  is  a  time  average  force  exerted  by  an  acoustic  field  on  an  object. 
This  force  is  an  example  of  a  universal  phenomenon  in  any  wave  motion  that  introduces  some 
type  of  unidirectional  force  on  absorbing  or  reflecting  targets  in  the  wave  path.  For  a  review  of 
this  topic  the  reader  may  refer  to  (Chu  and  Apfel  1982). 

Consider  a  plane  ultrasound  beam  interacting  with  a  planar  object  of  zero  thickness  and 
arbitrary  shape  and  boundary  impedance  that  scatters  and  absorbs.  The  radiation  force  vector, 
F,  arising  from  this  interaction  has  a  component  in  the  beam  direction  and  another  transverse 
to  it.  The  magnitude  of  this  force,  is  proportional  to  the  average  energy  density  of  the  incident 
wave  (E)  at  the  object,  where  ()  represents  the  time  average,  and  S,  the  area  of  the  projected 
portion  of  the  object  (Westervelt  1951): 


F  =  d  rS(E),  (1) 

where  dr  is  the  is  the  vector  drag  coefficient  with  a  component  in  the  incident  beam  direction 
and  another  transverse  to  it.  The  coefficient  dr  is  defined  per  unit  incident  energy  density  and 
unit  projected  area.  For  a  planar  object,  d,  is  numerically  equal  to  the  force  on  the  unit  area  of 
the  object  per  unit  energy  density.  Physically,  the  drag  coefficient  represents  the  scattering  and 
absorbing  properties  of  the  object  and  is  given  by  (Westervelt  1951) 


dr 


(na  +  us 


/ 


7  cosasdS 


7  sin  asdS, 


(2) 


where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and  normal  to  it,  respectively.  The 
quantities  IIa  and  Jl8  are  the  total  absorbed  and  scattered  powers,  respectively,  and  7  is  the 
scattered  intensity,  all  expressed  per  unit  incident  intensity.  Also,  as  is  the  angle  between  the 
incident  and  the  scattered  intensity,  and  dS  is  the  scalar  area  element.  The  drag  coefficient 
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can  also  be  interpreted  as  the  ratio  of  the  radiation  force  magnitude  on  a  given  object  to  the 
corresponding  value  if  the  object  were  replaced  by  a  totally  absorbing  object  of  similar  size.  For 
simplicity,  we  assume  a  planar  objected  oriented  perpendicular  to  the  beam  axis.  In  this  case, 
the  transverse  component  vanishes,  thus,  the  drag  coefficient  (force)  will  have  only  a  component 
normal  to  the  target  surface  which  we  denote  by  scalar  dr  ( F ). 

To  produce  a  time  varying  radiation  force,  the  intensity  of  the  incident  beam  can  be  shaped 
in  various  ways.  For  example,  a  short  ultrasound  pulse  can  produce  a  transient  pulsed  radiation 
force,  and  a  sinusoidally  modulated  beam  can  result  in  a  sinusoidally  varying  force. 

2.1.  Transient  methods 

In  the  transient  methods  the  radiation  force  of  ultrasound  is  used  to  make  a  minute  deforma¬ 
tion  in  the  tissue.  The  transient  recoil  of  the  tissue  resulting  from  this  deformation  is  measured 
and  used  for  evaluation  of  tissue  elastic  properties. 

A  method  for  measuring  tissue  hardness,  presented  by  Sugimoto  et  al  (1990),  uses  the 
radiation  force  of  a  single  focused  ultrasound  beam.  Ideally,  hardness  may  be  represented  by  the 
spring  constant  of  the  object,  which  is  the  ratio  of  the  applied  force  to  the  displacement.  Because 
it  is  difficult  to  measure  the  absolute  values  of  force  or  deformation  in  vivo,  this  method  relies 
on  evaluating  the  hardness  based  on  changes  in  relative  values.  In  this  method  an  ultrasound 
pulse  is  used  to  generate  a  short  duration  radiation  force  which  produces  localized  deformation 
of  the  tissue.  Immediately  after  the  force  pulse,  the  resulting  transient  deformation  of  tissue 
is  measured  as  a  function  of  time  with  Doppler  ultrasound  using  a  separate  transducer.  The 
deformation  includes  an  initial  rapid  squeeze  of  the  tissue,  followed  by  a  relaxation  and  possibly 
a  rebound.  This  function  is  a  function  of  tissue  viscoelastic  parameters,  as  well  as  the  applied 
force.  Sugimoto  et  al  (1990)  argue  that  because  measuring  the  internal  radiation  force  in  an 
absolute  sense  is  difficult,  it  is  advantageous  to  derive  a  relative  quantity  that  is  representative 
of  tissue  hardness.  To  derive  a  single  relative  quantity  from  the  deformation  data,  the  relaxation 
part  of  the  function  is  approximated  by  a  sum  of  several  exponential  curves,  and  the  sum  of 
the  first  order  derivatives  of  such  exponential  is  calculated.  It  is  shown  that  this  quantity  is 
correlated  to  the  spring  constant  of  the  tissue,  thus  may  be  used  as  a  measure  of  tissue  hardness. 

2.2.  Shear  wave  methods 

Shear  modulus  is  related  to  the  hardness  or  elasticity  of  the  material.  It  is  known  that  the  shear 
modulus  of  various  soft  tissues  range  over  several  orders  of  magnitude,  while  bulk  modulus,  a 
parameter  that  is  associated  with  the  conventional  pulse-echo  ultrasound  compressional  wave 
speed,  varies  significantly  less  than  an  order  of  magnitude  (Goss  et  al  1978;  Frizzell  and 
Carstensen  1976).  These  features  indicate  that  shear  modulus  may  be  a  better  parameter  for 
tissue  characterization  than  bulk  modulus.  One  way  to  induce  localized  shear  waves  inside 
tissue  is  to  use  the  radiation  force  of  focused  ultrasound  (Andreev  et  al  1997).  In  a  method 
called  shear  wave  elasticity  imaging  (SWEI)  (Sarvazyan  et  al  1998),  an  amplitude  modulated 
single  focused  ultrasound  beam  is  used  to  induce  a  localized  radiation  stress  inside  the  soft  tissue. 
Localization  of  the  stress  field  is  the  key  in  success  of  the  method.  To  achieve  a  high  degree 
of  localization,  the  method  uses  a  focused  ultrasound  beam.  It  is  shqwn  that  the  radiation  stress 
exerted  within  a  dissipative  medium  peaks  near  the  focal  region  of  highly  focused  transducer. 
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It  has  also  been  suggested  that  localization  can  be  improved  by  designing  the  transducer  and 
selecting  the  beam  parameters  such  that  a  nonlinear  shock  wave  is  produced  in  the  focal  region, 
increasing  the  magnitude  of  the  stress  field  in  the  vicinity  of  the  focal  region,  thus  augmenting 
localization  of  the  stress  field  (Sarvazyan  et  al  1998). 

Modulation  of  the  ultrasound  beam  can  be  in  the  form  of  an  oscillating  wave  or  short  pulse. 
The  resulting  radiation  force  elicits  a  shear  wave  propagating  in  the  radial  direction  with  respect 
to  the  beam  axis,  with  particle  motion  parallel  to  the  beam  axis.  Shear  waves  in  soft  tissue 
travel  at  very  low  speed,  typically  around  a  few  meters  per  second  (Frizzell  and  Carstensen 
1976;  Sarvazyan  et  al  1998),  thus  the  corresponding  wavelength  is  much  shorter  than  that  of 
the  compressional  waves  for  the  same  frequency.  Shear  waves  are  also  highly  attenuated  in  soft 
tissue,  with  an  attenuation  coefficient  two  or  three  orders  of  magnitude  higher  than  that  of  the 
compressional  waves  (Frizzell  and  Carstensen  1976;  Sarvazyan  et  al  1998).  Because  of  high 
attenuation  of  shear  waves,  it  is  possible  to  induce  them  in  a  very  limited  region  in  the  vicinity 
of  the  focal  point  of  the  ultrasound  beam,  hence  avoiding  the  influence  of  tissue  boundaries. 
In  SWEI,  a  short  ultrasound  pulse  is  delivered  to  the  tissue  causing  a  radial  shear  oscillation 
around  the  beam  axis.  By  measuring  the  amplitude  and  the  temporal  characteristics  of  this  wave, 
shear  parameters  of  the  tissue,  such  as  shear  modulus  and  shear  viscosity  can  be  calculated.  For 
example,  the  time  required  for  the  wave  front  to  propagate  from  one  point  to  another  can  be 
used  to  calculate  the  shear  wave  speed,  and  consequently,  the  shear  wave  modulus  / 1 i,  as 

/'  =  p4  (3) 

where  p  and  ct  are  the  density  and  the  shear  wave  velocity,  respectively. 

Shear  waves  may  be  detected  optically.  In  this  method  a  laser  source  and  a  photo  detector 
are  used  to  detect  the  displacement  of  particles  due  to  the  shear  wave  in  a  transparent  phantom 
(Sarvazyan  et  al  1998).  Because  this  method  requires  a  transparent  medium,  its  application  in 
vivo  is  difficult.  Phase  sensitive  magnetic  resonance  imaging  (Muthupillai  et  al  1995;  Sarvazyan 
et  al  1998)  is  an  alternative  method  that  can  be  used  to  measure  the  two-dimensional  distribution 
of  particle  displacement  in  a  given  direction  versus  time  in  a  material.  In  an  experiment 
presented  in  (Sarvazyan  et  al  1998),  an  ultrasound  pulse  of  3.6  ms  duration  produced  by  a 
70  mm  diameter  transducer  focused  at  100  mm  was  transmitted  within  a  cylindrical  rubber 
phantom.  The  displacement  was  measured  by  2.0  T  MRI  system  at  two  different  times  after  the 
acoustic  pulse  was  applied.  The  position  of  the  peak  displacement  at  these  time  points  was  used 
to  estimate  the  shear  velocity,  which  was  shown  to  be  consistent  with  the  independently  measured 
value.  Shear  waves  can  also  be  detected  by  Doppler  ultrasound  (Andreev  et  al  1997).  It  has  been 
shown  (Walker  1999)  that  to  achieve  an  appreciable  displacement  needed  for  Doppler  detection, 
most  soft  tissues  require  high  ultrasound  intensities  which  might  be  beyond  the  recommended 
maximum  intensity  level  for  diagnostic  ultrasound  (FDA  report,  see  reference  9). 

2.3.  Vibro-acoustography 

2.3.1.  Method 

This  technique  produces  a  map  of  the  mechanical  response  of  an  object  to  a  dynamic  force 
applied  at  each  point.  The  method  utilizes  ultrasound  radiation  force  to  remotely  exert  a  localized 


6 


oscillating  stress  field  at  a  desired  frequency  within  (or  on  the  surface  of)  an  object  and  records 
the  resultant  acoustic  response  (Fatemi  and  Greenleaf  1998;  Fatemi  and  Greenleaf  1999).  This 
acoustic  response,  which  is  normally  at  low  kHz  range,  is  a  function  of  the  viscoelastic  properties 
of  the  object  and  can  be  used  to  produce  an  image  of  the  object. 

To  produce  an  oscillating  radiation  force  the  intensity  of  the  incident  ultrasound  must  be 
amplitude  modulated  at  the  desired  low  frequency.  Using  a  single  amplitude  modulated  beam 
seems  to  be  the  simplest  means  to  attain  this  purpose.  However,  such  a  beam  will  exert  a 
radiation  force  on  any  object  that  is  present  along  the  beam  path,  producing  undesirable  acoustic 
emission.  To  confine  the  oscillating  radiation  stress  to  the  desired  region,  vibro-acoustography 
uses  two  unmodulated  CW  beams  at  slightly  different  frequency,  propagating  along  separate 
paths.  The  beams  are  arranged  to  cross  each  other  at  their  respective  foci,  and  thus  produce 
a  modulated  field  at  a  confined,  small,  cross-sectional  region.  The  single  and  double  beam 
arrangements  are  illustrated  in  Fig.  1. 

2.3.2.  System  diagram  and  image  formation 

Figure  2  illustrates  a  vibro-acoustography  system.  The  elements  of  the  annular  array  trans¬ 
ducer  are  driven  by  two  CW  signals  at  frequencies  fi  =  wi/27r  and  f\  +  A /  =  (wi  +  Aw)/27t. 
The  object  to  be  imaged  is  placed  at  the  joint  focal  plane  of  the  transducer  elements  (the  scanning 
plane).  The  ultrasound  field  produces  a  localized  radiation  stress  field  at  the  focal  point  on  the 
object  at  frequency  A/.  Depending  on  the  elastic  properties  of  the  object,  the  radiation  force 
may  cause  a  portion  of  the  object,  or  the  entire  object,  to  vibrate  at  this  frequency.  The  acoustic 
emission  resulting  from  object  vibration  is  received  by  a  hydrophone  which  is  located  nearby. 
Normally  the  wavelength  of  vibration  is  large  compare  to  the  object  size,  hence  the  acoustic 
emission  field  is  almost  omnidirectional.  Therefore,  hydrophone  position  is  not  a  critical  param¬ 
eter  in  this  measurement.  The  filter  is  used  to  eliminate  noise  and  other  interfering  signals.  To 
form  an  image  the  focal  point  of  the  transducer  is  moved  across  the  scanning  plane  on  the  object 
in  a  raster  pattern.  The  acoustic  emission  is  received  at  each  position,  and  an  image  is  formed 
by  displaying  the  amplitude  (or  phase)  of  such  signals  at  corresponding  positions  on  the  image 
plane.  The  spatial  resolution  of  this  imaging  method  is  determined  by  the  ultrasound  beamwidth 
at  the  focal  plane,  which  is  normally  in  the  order  of  the  incident  ultrasound  wavelength. 

2.3.3.  Theory 

In  the  beam  forming  method  to  be  described  the  amplitude  modulated  field  is  obtained  by 
the  interference  of  two  unmodulated  ultrasound  beams.  Radially  symmetric  interfering  beams 
are  obtained  when  two  coaxial,  confocal  transducers  are  used  (Fatemi  and  Greenleaf  1999).  For 
this  purpose,  elements  of  a  two-element  spherically  focused  annular  array  (consisting  of  a  central 
disc  with  radius  a\  and  an  outer  ring  with  the  inner  radius  of  n!2  and  outer  radius  of  0,2)  are 
excited  by  separate  CW  signals  at  frequencies  uj\  =  wo  —  Aw/2  and  W2  =  wo  +  Aw/2.  We 
assume  that  the  beams  are  propagating  in  the  +z  direction  with  the  joint  focal  point  at  2  =  0. 
The  resultant  field  on  the  2  =  0  plane  may  be  written  as: 


p(t)  =  Pi(r)  cos  (wi<  +  V’l  O'))  +  P2O*)  cos  (w2 1  +  ^2  O')), 


(4) 
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where  r  =  \f:r?  +  y2  is  the  radial  distance.  The  amplitude  functions  are  (Morse  and  Ingard 
1968;  Kino  1987) 


Pi(r)  =  pcUoi 


7r  a2 
Ai^o 


(5) 


and 


P2(r)  =  pcUQ2 


7T 

A2^0 


a\  jinc 


(  ra 2 

\A2^0 


—  02  jinc 


(6) 


where  Uq,  is  the  velocity  amplitude  at  the  i-th  transducer  element  surface,  and  A*  =  2i r/w,;  for 
i=l,2  are  the  ultrasound  wavelengths.  The  phase  functions, 
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for  i=7,  2,  are  conveniently  set  to  be  zero  at  the  origin.  Also,  jinc(X)  =  Ji(2ttX)/ttX,  where 
,/i(.)  is  the  first  order  Bessel  function  of  the  first  kind. 

The  instantaneous  energy  is  E  =  p2(t)/pc2.  Replacing  p(t)  from  Eq.  4,  this  energy  will 
have  a  time  independent  component,  a  component  at  the  difference  frequency  Au;  =  lo2  —  u>i 
which  results  from  the  cross  product  of  the  two  pressure  fields,  and  high  frequency  components 
at  oj\  and  u>2  and  their  harmonics.  The  energy  component  at  the  difference  frequency  is: 


CAw(*)  = 


f>i(ro)P2(^o) 

pc2 


cos  [A u)t  +  Aip(r<))], 


(8) 


where  AV'(r’o)  =  ^2(^0)  —  V’l^o)- 

Now,  we  define  a  unit  point  target  with  an  area  of  dxdy  at  position  (xo,yo)  on  the  focal 
plane,  and  with  a  drag  coefficient  dr(x 0,  ?/o)  such  that:  dr(x,Q:  yo)dxdy  =  1  on  the  target  and  zero 
elsewhere.  This  equation  is  merely  used  as  a  mathematical  model.  In  this  case,  the  projected  area 
can  be  considered  to  be  S  —  d,xd,y.  Therefore,  if  the  projected  area  is  unity,  then  dr(x,y )  =  1, 
which  according  to  Eq.  2  corresponds  to  a  totally  absorptive  object. 

Referring  to  Eq.  1,  and  replacing  d,rS  with  unity  and  ( E )  by  e&u(t)  of  Eq.  8,  then  we  can 
write  the  low-frequency  component  of  the  radiation  force  on  the  unit  point  target  as: 


/au,(*0,  2/0;  t)  =  — 2Pl(ro)P2(ro)  COS  [A ut  +  A^(ro)], 


pc 


(9) 


where  arguments  xq  and  yo  are  added  to  denote  the  position  of  the  point  target,  and  ro  = 
■sJx.q  +  ijq.  Referring  to  Eqs.  5  and  6,  the  complex  amplitude  of  the  stress  field  can  be  found  as 
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where  A  A  =  2ttc/Auj  is  the  wavelength  associated  with  A u>.  The  above  equation  indicates  that 
the  radiation  stress  is  concentrated  at  the  focal  point  and  decays  down  quickly  with  distance 
as  jinc(.)2. 

As  an  example,  we  consider  a  confocal  transducer  with  dimensions  a\  =  14.8  mm,  a'2  =  16.8 
mm,  a2  =  22.5  mm,  and  the  focal  length  70  mm.  Also,  we  assume  that  the  center  frequency 
is  3  MHz,  and  the  difference  frequency  is  7.3  kHz.  The  radiation  stress  at  the  focal  plane  of 
a  this  transducer  is  plotted  in  Fig.  3. 

In  medical  applications  the  maximum  ultrasonic  intensity  is  regulated  for  safety  reasons.  It 
is,  therefore,  useful  to  write  the  stress  field  in  terms  of  the  peak  ultrasonic  intensity  at  the  focal 
point.  The  long-term  average  of  ultrasonic  intensity  at  the  focal  point  can  be  written  as 


m  =  p2(°)  +  ^1(0) 

^  ’  2  pc 


(11) 


where  -Pi(O)  and  ^2(0)  can  be  found  from  Eqs.  5  and  6  as 


Assuming  Uq\  =  C/02  =  Uq,  we  can  write  the  focal  plane  stress  field  in  terms  of  1(0): 
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The  first  fraction  in  the  above  equation  represents  the  static  radiation  force  produced  by  the  two 
beams  on  a  total  absorber,  the  second  fraction  is  a  constant  factor,  and  the  rest  represents  the 
spatial  distribution  of  the  stress  field  on  the  focal  plane.  If  Au>  <C  u> 2,  then  we  may  replace 
Ai  and  A2  by  with  Ao,  and  simplify  the  expression  for  the  stress  field.  Under  these  conditions, 
the  stress  at  the  focal  point  is 


^Aw(0,0)  = 


2/(0)  aj(fl2  -°42) 


x 


(14) 


oi+(4-a?y 

The  fraction  on  the  right  represents  the  effect  of  transducer  dimension  on  the  stress  field.  For  the 
transducer  used  in  the  previous  example,  the  value  of  this  fraction  is  0.4999.  (This  fraction  can 
be  also  calculated  for  a  single  element  transducer  of  the  same  diameter,  excited  by  an  amplitude 
modulated  signal.  For  this  purpose,  we  may  let  a\  =  0.2  and  a2  =  0.  For  these  values,  the  last 
fraction  in  Eq.  14  is  0.5.)  The  resulting  radiation  stress  (assuming  c  =  1500  m/sec  for  water) 
is  Faw(0,  0)  =  6.67  x  10_4/(0)  N/m2.  Now  letting  1(0)  =  7200  W/m2,  which  is  the  maximum 
spatial  peak,  time  average  intensity  suggested  by  FDA  for  in  vivo  applications  (FDA  report,  see 
reference  9),  the  resulting  stress  at  the  focal  point  is  /  aw(0,  0)  =  4.80  N/m2.  Referring  to  Fig.  3, 
we  note  that  this  stress  field  is  applied  only  in  a  small  region  around  the  focal  point  to  the  object. 


9 


2.3.4.  Acoustic  emission 


To  explain  the  acoustic  emission  we  consider  an  “object”  within  a  homogeneous  infinite 
medium.  This  model  allows  us  to  separate  the  roles  played  by  the  parameters  of  the  object  and 
the  surrounding  medium,  and  can  be  fitted  to  various  applications.  When  an  oscillating  stress 
field  is  applied  to  the  object,  the  object  vibrates  at  the  frequency  of  the  stress  field.  Vibrational 
energy  of  the  object  is  partly  transferred  to  the  surrounding  medium,  which  is  called  the  acoustic 
emission  field.  Here,  we  calculate  the  acoustic  emission  of  an  object  resulting  from  applying 
the  radiation  stress. 

To  explain  the  physics  in  more  generalized  terms,  for  the  moment  we  consider  an  object  in 
the  form  of  a  flat  disk  facing  the  beam.  Here,  we  assume  that  the  vibrating  object  has  a  circular 
cross  section  of  radius  b  and  uniformly  vibrates  back  and  forth  like  a  piston.  This  choice  allows 
us  to  illustrate  the  concept  in  a  simple  form.  We  also  consider  an  area  S  <  7rfr2  of  the  piston 
surface  to  be  projected  normally  by  the  beam.  We  can  always  return  to  our  elementary  point 
object  by  reducing  the  area  of  this  disk  to  dxdy.  Similar  solutions  can  be  carried  out  for  objects 
of  other  forms.  The  theory  can  be  also  extended  to  include  arbitrary  vibrating-part  shapes  and 
nonuniform  displacement  of  the  object.  The  total  radiation  force  on  this  object  F^,  can  be  found 
by  integrating  the  radiation  stress  over  the  area  of  the  object.  This  force  vibrates  the  target  object 
at  frequency  A uj.  The  steady  state  normal  velocity  amplitude  of  a  piston  at  frequency  Aw,  U^, 
due  to  a  harmonic  force  Fau,  can  be  written  as 


Uauj 


Falo 

Zaw 


(15) 


where  Z is  the  mechanical  impedance  of  the  object  at  Aw.  The  mechanical  impedance  of 
the  object  has  two  components,  one  resulting  from  the  inertia,  friction,  and  the  elasticity  of  the 
object  itself,  and  the  other  resulting  from  the  loading  effect  of  the  surrounding  medium  on  the 
vibrating  object.  The  mechanical  impedance  can  be  interpreted  as  a  measure  of  object  rigidity 
and  how  much  it  yields  to  the  applied  force.  For  example,  for  a  rigid  object,  Zau  is  high,  and 
hence  the  object  resists  the  force. 

Knowing  Uau,  we  can  calculate  the  pressure  field  it  produces  in  the  medium.  We  assume 
that  the  acoustic  emission  signal  propagates  in  a  free  and  homogenous  medium.  The  farfield 
acoustic  pressure  due  to  a  piston  source  of  radius  b  set  in  a  planar  boundary  of  infinite  extent 
is  given  by  (Morse  and  Ingard  1968), 


Pacj  =  -J  Awp 


exp  (jAuil/c) 
47 rl 


2  J\  ( sin  i9)  cos  d 
_  -  X  - 


Ac jb 
c 


sin$ 


cos  d  +  (3b 


(2t Tb2UAu),  (16) 


where  l  is  the  distance  from  the  observation  point  to  the  center  of  the  piston,  t?  is  the  angle 
between  this  line  and  the  piston  axis,  and  (3  b  is  the  specific  acoustic  admittance  of  the  boundary 
surface  (The  specific  acoustic  admittance  is  f3B  =  J|,  where  ZB,  the  acoustic  impedance  of  the 
boundary,  represents  the  ratio  between  the  pressure  and  normal  fluid  velocity  at  a  point  on  the 
object).  The  factor  two  in  the  parenthesis  comes  from  the  presence  of  the  boundary  wall.  It 
would  be  replaced  by  unity  if  the  boundary  wall  were  not  present  (Morse  and  Ingard  1968). 
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The  acoustic  emission  field  resulting  from  object  vibration  can  be  written  in  terms  of  object 
mechanical  impedance  by  combining  Eqs.  15  and  16,  as 

_  2  f  •  . .  exp  (jAu>l/c)  [ 2 h  (^r  sin  j)  cos  t  1  \ 

Aw  ^  r  C2  4tt/  sin $  cos d  +  (3B  [  .... 


■Faw  . 


For  wavelengths  long  compared  to  the  object  size,  i.e.,  when  bAco/c  — >  0,  the  term  in  the  bracket 
approaches  a  constant  because  2J(X)/X  — >  1  for  X  — *  0.  Hence,  ignoring  the  first  fraction  in 
the  brackets,  we  may  consider  the  contents  of  the  brace  to  be  an  object  independent  function 
(the  specific  acoustic  admittance  Pb  relates  to  the  surrounding  boundary  surface).  Under  these 
conditions,  the  brace  in  the  above  equation  represents  the  effect  of  the  medium  on  the  acoustic 
emission  field,  which  we  may  call  the  medium  transfer  function,  and  denote  it  by 


tt  ,A^  exp  (j Awl /c)  w 

—  J  2  x  x 


cos  P  +  Pb 


The  last  fraction  in  Eq.  17  includes  ,  which  is  the  mechanical  admittance  of  the  object  at 
the  frequency  of  the  acoustic  emission  (Au),  and  we  denote  this  by  Yalj-  It  is  convenient  to 
combine  this  term  with  the  next  term  (27r&2)  in  Eq.  17,  as  Qau  =  27t62Yao>  =  2Trb2/Z^,  which 
is  the  total  acoustic  outflow  by  the  object  per  unit  force  (acoustic  outflow  is  the  volume  of  the 
medium,  e.g.,  the  fluid,  in  front  of  the  object  surface  that  is  displaced  per  unit  time  due  to  object 
motion).  Function  Qauj  represents  the  object  characteristics  at  the  acoustic  frequency.  We  may 
thus  rewrite  Eq.  17  in  a  more  compact  form  as 


Paoj  —  PC  H AojO'^Q AujF Aco¬ 


rn 


Equation  19  indicates  that  the  acoustic  emission  pressure  is  proportional  to:  (1)  the  radiation 
force,  which  itself  is  proportional  to  the  square  of  ultrasound  pressure  and  the  ultrasound 
characteristics  of  the  object,  dr,  and  the  projected  area  5;  (2)  the  acoustic  outflow  by  this  object, 
Qa-jj,  representing  the  object  size  b  and  its  mechanical  admittance  at  the  acoustic  frequency, 
Yaw',  and  (3)  the  transfer  function  of  the  medium  at  the  acoustic  frequency,  #Aw(0-  Note 
the  difference  between  the  projection  area  S  and  the  vibrating  area  irb2.  The  projection  area 
determines  the  extent  of  the  force  applied  to  the  object  (Eq.  1).  The  vibrating  area,  however, 
influences  the  total  acoustic  outflow  in  the  medium  caused  by  object  vibration. 


2.3.5.  Measuring  the  acoustic  emission 

A  hydrophone  is  used  to  measure  the  acoustic  emission  field.  Equation  18  indicates  that  for 
long  wavelengths  compare  to  the  object  size,  the  acoustic  emission  field  pattern  resembles  that  of 
a  dipole  with  an  angle  dependency  represented  by  the  last  fraction  to  the  right  of  this  equation. 
This  angle  dependency  is  normally  weak  and  hence  the  position  of  the  receiving  hydrophone 
is  not  critical.  Also,  because  the  attenuation  of  low  frequency  acoustic  waves  in  tissue  and 
liquids  is  normally  low,  one  may  place  the  hydrophone  at  a  distance  from  the  object  without 
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losing  too  much  intensity.  The  simplified  diagram  presented  in  Fig.  2  describes  how  to  measure 
the  amplitude  of  the  acoustic  emission  field.  To  measure  the  phase  of  the  acoustic  emission,  a 
phase  reference  needs  to  be  set.  One  way  to  accomplish  this  is  to  generate  a  reference  signal 
by  electronic  down  mixing  of  the  two  CW  excitation  signals  to  produce  a  beat  signal  at  the 
difference  frequency  A/.  Then,  using  this  reference  signal,  the  phase  of  the  acoustic  emission 
signal  (hydrophone  output)  can  be  calculated  by  conventional  phase  detection  methods. 

2.3.6.  Applications 

Evaluation  of  the  characteristics  of  an  object  (or  a  medium)  by  listening  to  its  sound  is  a 
traditional  approach  that  has  been  used  for  many  purposes.  Qualitative  evaluation  of  a  crystal 
glass  by  tapping  on  it  is  a  simple  example.  Vibro-acoustography  implements  the  same  approach 
but  in  micro-scale,  and  in  a  way  that  could  be  applied  to  tissues. 

Equation  19  presents  a  general  relationship  between  the  mechanical  parameters  of  the 
objects,  surrounding  medium,  and  the  acoustic  emission  field  resulting  from  object  vibration. 
By  measuring  the  acoustic  emission  field,  it  would  be  possible,  in  principle,  to  estimate  some 
of  the  object  or  medium  parameters,  either  in  an  absolute  or  in  a  relative  sense.  For  example, 
one  can  use  vibro-acoustography  to  measure  the  resonance  frequency  of  an  object,  and  from  that 
information  it  is  possible  to  estimate  some  viscoelastic  parameters  of  the  object  or  the  medium. 
This  method  has  been  used  to  measure  the  Young’s  modulus  of  a  metallic  rod  (Fatemi  and 
Greenleaf  1999a).  In  another  experiment,  by  measuring  the  resonance  frequency  of  a  known 
resonator  in  a  liquid  medium,  the  viscosity  of  the  fluid  has  been  estimated  with  good  accuracy 
(Fatemi  and  Greenleaf  1999b).  These  applications  are  not  necessarily  medical,  but  the  principle 
could  be  used  for  evaluation  of  soft  tissues,  blood,  bone,  etc. 

In  medical  applications,  one  can  use  vibro-acoustography  to  obtain  an  image  of  the  object 
for  diagnostic  purposes.  In  such  applications,  the  image  may  not  represent  a  physical  quantity, 
rather  it  provides  a  means  to  visualize  object  details. 

2.3.7.  Experimental  results 

In  the  following  experiment  vibro-acoustography  has  been  used  to  image  excised  human 
carotid  arteries  that  include  calcifications.  Characteristics  of  the  ultrasound  transducer  used  here 
are  those  described  in  the  example  presented  in  Theory  section.  The  object  consists  of  two 
excised  human  carotid  arteries  secured  on  a  thin  sheet  of  latex.  This  sheet  is  almost  transparent 
to  the  ultrasound  and  x-ray.  Figure  4a  shows  an  x-ray  image  of  these  arteries.  The  left  artery 
is  from  a  young  person  and  has  shows  no  calcification.  The  left  artery  is  from  an  older  person 
and  includes  some  calcium  deposit  which  can  be  seen  as  a  bright  region  near  the  bifurcation. 
The  arteries  were  immersed  in  a  water  tank  for  scanning.  The  difference  frequency  A /  in  this 
experiment  is  7  kHz.  The  lead  number  “2”  was  placed  on  the  latex  sheet  for  identification 
purposes.  Figure  4b  shows  the  vibro-acoustography  image  (displaying  the  amplitude  of  the 
acoustic  emission  signal)  of  the  object.  In  this  image,  the  calcified  region  is  seen  as  a  bright 
region,  resembling  its  x-ray  image.  Calcification,  which  is  formed  from  stiff  (compare  to  the 
arterial  wall)  material,  is  an  efficient  acoustic  radiator,  producing  strong  acoustic  emission  when 
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exposed  to  the  radiation  force  of  the  incident  ultrasound.  Hence,  it  stands  out  in  the  vibro- 
acoustography  image.  This  image  has  a  high  spatial  resolution  (about  0.7  mm),  low  speckle, 
and  high  signal-to-noise  ratio.  The  number  “2”,  which  is  also  made  from  a  hard  material,  is 
seen  with  high  contrast.  This  experiment  demonstrates  the  utility  of  the  method  in  delineating 
stiff  material  in  a  soft  background. 

3.  Comparisons  of  the  methods 

In  this  paper  we  presented  three  approaches  for  using  localized  ultrasound  radiation  force  to 
probe  biological  materials.  Here  we  compare  the  features  and  capabilities  of  the  methods. 

3.1.  Objectives 

The  transient  method  estimates  a  quantity  that  correlates  to  the  hardness  of  tissue.  While 
it  does  not  estimate  a  specific  physical  parameter,  the  estimated  parameter  is  quantitative  and 
may  be  useful  for  diagnosis. 

Shear  wave  methods  aim  specifically  at  the  shear  wave  characteristics  in  soft  tissue,  in 
particular,  estimating  the  shear  wave  velocity.  Shear  wave  velocity  in  soft  tissue  is  very  low 
(typically  in  the  order  of  a  few  m/sec),  allowing  localized  measurement  of  the  parameter.  Also, 
because  shear  wave  attenuation  in  soft  tissue  is  high,  the  boundaries  of  the  object  have  minimal 
effect  on  the  measurement.  This  method  is  not  intended  for  use  in  hard  materials. 

Vibro-acoustography  seeks  a  two-fold  objective.  In  this  method,  the  acoustic  emission  signal 
is  dependent  on  both  microstructure  (local  drag  coefficient)  and  macrostructure  (the  mechanical 
admittance)  of  the  object.  The  acoustic  emission  signal  can  be  analyzed  to  estimate  different 
characteristics  of  the  object.  In  some  applications,  both  the  local  parameters  and  those  that  are 
related  to  the  large-scale  structure  of  an  object  are  needed.  An  instructive  example  presented 
in  Fatemi  and  Greenleaf  (1999)  demonstrates  that  vibro-acoustography  can  be  used  to  produce 
an  image  of  a  tuning  fork  displaying  its  detail  (microstructures)  while  measuring  its  resonance 
frequency  (a  parameter  related  to  the  macrostructure)  represented  by  color. 

3.2.  Excitation  method 

In  the  first  two  methods  the  radiation  force  is  produced  by  a  single  focused  transducer  which 
is  driven  by  an  amplitude  modulated  signal,  in  particular,  pulse  shaped  (the  shear  wave  method 
can,  in  principle,  operate  using  a  single  sinusoidally  modulated  beam).  Vibro-acoustography, 
however,  uses  two  CW  excitation  signals  (or  long  tone  bursts  of  CW  signals).  The  major 
difference  in  these  excitation  techniques  is  the  distribution  of  the  radiation  stress  in  the  depth 
direction.  A  possible  drawback  of  using  a  single  amplitude  modulated  beam  method  is  that  the 
pulsed  ultrasound  exerts  a  radiation  stress  on  all  tissue  through  which  it  propagates.  Also,  a  shock 
can  be  produced  on  the  transducer  surface,  and  then  transmitted  to  the  tissue.  Such  phenomena 
result  in  extraneous  motion  of  the  object  that  could  introduce  errors  in  the  measurements.  In 
vibro-acoustography,  the  two  CW  beams  pass  through  separate  paths.  Outside  of  their  interaction 
at  the  focal  region,  they  can  only  produce  a  static  radiation  force  on  the  object,  which  produces 
no  sound.  The  beams  mix  near  the  focal  point  to  produce  a  modulated  field  that  can  generate  an 
oscillating  radiation  force  only  in  that  location.  The  resulting  acoustic  field,  therefore,  is  directly 
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related  to  the  excitation  of  the  object  at  the  focal  area.  This  is  an  important  feature,  because 
the  acoustic  emission  from  other  parts  of  the  object  or  the  transducer  should  not  be  allowed  to 
interfere  with  the  acoustic  emission  from  the  focal  area. 

3.3.  Detection  method 

The  transient  method  and  the  shear  wave  method  are  based  on  localized  detection  of  material 
displacement.  This  is  implemented  by  using  either  Doppler  ultrasound,  optical  displacement 
measurement,  or  a  more  sophisticated  method  such  as  MRI.  Vibro-acoustography  uses  an  audio 
hydrophone  or  microphone  for  detection,  which  is  a  much  simpler  detection  method  than  Doppler 
or  MRI. 

3.4.  Sensitivity 

Detection  sensitivity  is  critical  to  the  success  of  any  of  the  above  methods.  More  ultrasound 
intensity  would  be  needed  if  the  detection  sensitivity  is  poor.  Ultrasound  Doppler  methods,  at 
conventional  ultrasound  frequencies,  probably  can  detect  displacements  in  the  order  of  a  few 
micrometers.  The  sensitivity  of  the  MRI  technique  is  in  the  order  of  100  nm  (Muthupillai  et 
al  1995).  Vibro-acoustography  has  been  shown  to  detect  motions  as  small  as  a  few  nanometers 
(Fatemi  and  Greenleaf  1998).  The  high  sensitivity  of  this  method  is  a  result  of  the  fact  that  small 
motions  of  the  object  can  produce  an  acoustic  emission  pressure  field  that  is  easily  detectable 
by  a  sensitive  hydrophone. 

3.5.  Measured  quantities 

The  quantity  measured  in  the  transient  method  is  somewhat  arbitrary,  and  not  a  direct 
measure  of  any  physical  quantity.  The  merit  of  this  method  is  that  tissue  hardness  is  measured 
from  relative  displacement  values,  without  the  need  to  measure  the  applied  force.  The  shear  wave 
methods  directly  measure  a  physical  parameter  of  the  material  in  an  absolute  sense.  The  acoustic 
emission  field  measured  in  vibro-acoustography  is  a  function  of  several  physical  parameters  that 
relate  to  the  object  and  the  surrounding  medium.  It  is  possible  to  measure  a  single  parameter  only 
if  one  has  enough  knowledge  about  other  parameters  involved.  For  example,  one  can  measure 
the  shear  viscosity  of  a  liquid  using  a  known  resonator  as  an  object  in  the  liquid  (Fatemi  and 
Greenleaf  1999b),  or,  for  a  given  the  geometry,  the  Young’s  modulus  of  a  metallic  object  can 
be  measured  (Fatemi  and  Greenleaf  1999a)  by  vibro-acoustography. 

3.6.  Exposure  safety 

The  shear  wave  method  presented  in  (Sarvazyan  et  al  1998)  uses  high  intensity  pulsed 
ultrasound  for  excitation  and  either  the  optical  or  MRI  methods  for  detection.  It  is  argued  that, 
although  the  intensity  is  high,  the  total  exposure  is  within  the  FDA  limits  (FDA  report,  see 
reference  9).  These  limits  are  recommended  based  on  possible  thermal  and  mechanical  effects 
of  ultrasound  on  tissue.  However,  in  vivo  use  of  such  detection  methods  would  be  difficult. 
Using  Doppler  ultrasound  for  detection  is  less  sensitive,  hence,  requires  higher  peak  ultrasound 
intensity  to  produce  enough  displacement  of  the  tissue.  This  intensity  level  may  exceed  the 
recommended  spatial  peak  pulse  average  intensity  limit. 
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Vibro-acoustography  uses  either  CW  or  tone  burst  ultrasound.  Therefore,  the  continuous 
ultrasound  intensity  (spatial  peak  time  average  intensity)  limit  must  be  observed.  This  limit 
according  to  the  FDA  is  720  mW/cm2  (FDA  report,  see  reference  9).  However,  because  of  high 
sensitivity  of  the  hydrophone  detector,  it  is  possible  to  detect  very  small  levels  of  the  acoustic 
emission  while  using  low  transmit  power.  For  example,  it  has  been  shown  that  a  submillimeter 
object  can  be  detected  at  ultrasound  intensities  far  below  the  FDA  limit  (Fatemi  and  Greenleaf 
1999). 

3.7.  In  vivo  applicability 

In  vivo  application  of  the  transient  method  is  not  clear  because  it  requires  high  ultrasound 
power  to  produce  a  detectable  displacement.  Application  of  the  shear  wave  methods  in  human 
body  largely  depends  on  the  detection  method  used.  Optical  methods  are  not  likely  usable  in 
vivo.  MRI  detection  methods  are  generally  complex,  especially  when  the  operation  is  to  be  linked 
to  the  ultrasound  system  in  the  magnetic  field.  Ultrasound  Doppler  is  a  more  practical  detection 
method  for  shear  wave  detection;  however,  its  sensitivity  may  not  be  sufficient  for  ultrasound 
intensities  within  the  recommended  maximum  level  (Walker  1999).  Vibro-acoustography  uses 
a  hydrophone  for  detection  which  is  simple  to  operate  in  clinical  settings.  It,  however,  requires 
an  acoustically  quiet  environment  for  proper  detection.  The  biological  noise  spectrum  seems  to 
be  below  1  kHz,  which  can  be  easily  filtered  out  if  the  operation  frequency  is  above  this  limit 
(Fatemi  and  Greenleaf  1999). 

3.8.  Image  resolution 

Imaging  complex  objects  using  the  transient  or  the  shear  wave  methods  has  not  been  fully 
explored  in  the  literature.  Vibro-acoustography  can  be  used  for  high  quality  imaging.  The  spatial 
resolution  is  proportional  to  the  width  of  the  stress  field  main  lobe.  For  example,  for  the  3  MHz 
ultrasound  transducer  presented  in  Fig.  3,  the  spatial  resolution  is  about  700  micrometers. 

3.9.  Material  properties 

Both  the  transient  and  shear  wave  methods  are  based  on  local  deformation  of  soft  tissue 
in  the  beam  direction.  These  methods  work  best  if  the  material  under  test  supports  the  shear 
wave,  has  enough  attenuation  to  allow  the  build  up  of  enough  radiation  force,  and  is  compliant 
enough  to  allow  appreciable  displacement.  Application  of  these  methods  for  hard  material,  such 
as  bone  and  calcifications,  would  be  difficult  or  impractical  because:  (a)  stiff  materials  have  low 
compliance,  hence,  their  displacement  in  response  to  the  force  is  relatively  small  and  difficult 
to  detect;  (b)  the  shear  wave  travels  at  high  speed  in  such  material.  In  vibro-acoustography 
hard  materials  produce  structural  vibrations  which  are  often  stronger  than  that  of  soft  tissue. 
Such  vibrations  result  in  strong  and  easily  detectable  acoustic  emissions.  The  bright  appearance 
of  the  arterial  calcifications  in  the  experiment  presented  in  Fig.  4  is  a  demonstration  of  this 
phenomenon. 

In  contrast  to  the  other  two  methods,  vibro-acoustography  can  be  used  to  detect  particles  in 
materials  that  do  not  support  share  waves,  for  example,  detecting  gas  bubbles  in  liquids.  The 
acoustic  emission  resulting  from  particle  vibration  includes  a  compressional  wave  component  that 
can  travel  in  the  surrounding  medium  including  the  liquid.  Tissue  attenuation  for  compressional 
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waves  at  low  frequencies  is  small,  hence  such  waves  can  be  detected  by  the  hydrophone  from 
a  distance.  Shear  waves  in  liquid  medium  decay  very  rapidly  and  are  difficult  to  detect  by 
most  methods. 

4.  Summary 

Recent  progress  on  techniques  for  evaluation  of  dynamic  viscoelastic  parameters  of  tissue, 
based  on  the  radiation  force  of  ultrasound,  has  prompted  us  to  view  these  techniques  as  an 
emerging  new  field.  We  recognize  the  capabilities  and  potentials  of  this  new  field  as  unique  and 
distinguishable  from  those  of  conventional  elasticity  or  ultrasound  imaging  techniques. 

Three  methods  for  probing  biological  tissue  using  the  dynamic  radiation  force  of  ultrasound 
were  presented.  The  first  method  measures  a  minute  transient  tissue  deformation  versus  time.  The 
second  method,  also  uses  a  transient  excitation,  but  measures  the  resulting  shear  wave  amplitude 
and  velocity.  The  last  method,  vibro-acoustography,  measures  the  acoustic  field  resulting  from 
object  vibration  at  a  specified  frequency.  Capabilities  and  limitations  of  these  methods  for 
probing  different  materials  and  their  potentials  for  in  vivo  applications  were  discussed. 
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7.  Legends 


Figure  1  Generation  of  a  modulated  ultrasound  field  by  single  and  confocal  beams.  The  single 
beam  produces  a  field  that  is  modulated  throughout  the  beam  path.  This  field  can  produce  an 
oscillating  radiation  force  on  any  object  in  the  beam  path.  The  confocal  annular  array  transducer 
generates  two  CW  beams  propagating  separate  paths.  These  beams  interfere  near  the  focal  zone, 
producing  a  modulated  field  only  in  that  region.  Hence,  the  oscillating  radiation  force  can  only 
be  generated  on  objects  in  the  beam  crossing  region. 


Figure  2  Vibro-acoustography  system.  The  two  elements  of  the  confocal  annular  array  are 
excited  by  separate  CW  signals  at  two  different  frequencies.  The  object  is  placed  at  the  focal 
plane  of  the  transducer.  The  ultrasound  field  produces  a  radiation  force  at  the  difference  frequency 
A /  on  the  object.  This  force  vibrates  the  object,  which  in  turn  produces  a  sound  field  in  the 
medium.  This  sound  field  is  received  by  an  audio  hydrophone  placed  near  the  object.  To  form 
an  image,  the  beams  are  scanned  in  a  raster  motion  across  the  scanning  plane,  and  the  resulting 
acoustic  emission  is  recorded  from  each  point.  Brightness  of  each  corresponding  point  on  the 
image  plane  is  determined  by  the  amplitude  (or  phase)  of  the  recorded  acoustic  emission. 


Figure  3  Three-dimensional  plot  of  the  stress  field  amplitude  of  a  confocal  annular  array 
transducer  at  the  focal  plane  at  a  point  in  time.  Transducer  dimensions  are:  a\  =  14.8  mm, 
o/2  =  16.8  mm,  a 2  =  22.5  mm,  and  the  focal  length  is  70  mm.  The  center  frequency  is  3  MHz, 
and  the  difference  frequency  is  7.3  kHz. 


Figure  4  Vibro-acoustography  of  excised  human  carotid  arteries:  (a)  x-ray  image  showing  a 
normal  (left)  and  a  calcified  (right)  artery  placed  on  a  sheet  of  latex;  (b)  vibro-acoustography  of 
the  arteries  in  (a)  measured  at  a  vibration  frequency  of  7  kHz.  The  calcification  is  seen  as  a  bright 
region  in  both  x-ray  and  vibro-acoustography  images.  Lead  number  “2”  is  used  for  identification. 
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Vibro-acoustography  is  a  new  imaging  method  that  utilizes  the  radiation  force  of 
ultrasound  to  vibrate  objects  at  low  (kHz  range)  frequencies.  The  acoustic  emission 
resulting  from  such  vibration  is  detected  by  a  hydrophone  and  is  used  to  obtain  an  image 
that  is  a  function  of  object  elastic  characteristics  [Proc  Natl  Acad  Sci  USA  96.6603-6608, 
June  1999].  Elastic  parameters,  such  as  object  stiffness,  are  valuable  in  delineating  hard 
objects  such  as  breast  micro  calcifications  from  the  surrounding  soft  tissue.  This  method 
offers  a  high  spatial  resolution  provided  by  two  focused  ultrasound  beams,  and  produces  a 
highly  localized  dynamic  stress  field  to  probe  the  object  point-by-point  with  low-frequency 
elastic  waves.  The  principles  used  for  beam  forming  in  vibro-acoustography  are 
fundamentally  different  from  those  for  the  conventional  ultrasound  imaging  methods.  The 
reason  is  that  in  vibro-acoustography  the  ultrasound  field  is  designed  to  interact 
nonlinearly  with  the  object  whereas  in  the  conventional  ultrasound  they  interact  linearly.  In 
order  to  understand  the  roles  of  different  parameters  of  the  system  and  the  object  in  the 
final  image  it  is  important  to  analyze  the  input-output  relationship  of  the  system. 

The  approach  taken  in  this  paper  is  to  derive  an  analytical  model  for  vibro-acoustography 
system.  This  model  includes  three  parts:  (a)  the  incident  beam  characteristics  and  the 
resulting  radiation  stress  field,  (b)  a  model  for  the  mechanical  response  of  the  object  to  the 
stress  field,  and  (c)  a  model  for  the  propagation  medium.  Using  this  model  we  can 
calculate  the  acoustic  emission  field  from  an  object  such  as  a  sphere  or  a  disc.  Also  this 
model  allows  us  to  calculate  the  complex  point  spread  function  of  the  imaging  system. 

The  results  show  that  the  acoustic  emission  is  inversely  proportional  to  the  impedance  of 
the  object,  which  is  comprised  of  two  components,  the  mechanical  impedance  and  the 
radiation  impedance.  These  impedances  can  be  written  in  terms  of  object  mass,  shape  and 
elastic  parameters.  The  quadratic  relationship  between  the  ultrasound  pressure  and  the 
resulting  acoustic  emission  pressure  is  determined  theoretically  and  verified 
experimentally.  To  verify  the  theoretical  point  spread  function,  the  experimental  point 
spread  function  is  obtained  and  shown  to  be  in  agreement  with  the  theoretical  results. 

We  conclude  that  the  model  presented  here  can  adequately  represent  the  imaging  system. 

The  U.S.  Army  Medical  Research  and  Materiel  Command  under  DAMD1 7-98- 1-8 121 
supported  this  work. 
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1.  Introduction 

It  is  well  known  that  changes  in  elasticity  of  soft  tissues  are  often  related  to  pathology. 
Traditionally,  physicians  use  palpation  as  a  simple  method  for  estimating  mechanical  properties 
of  tissue.  In  palpation,  a  static  force  is  applied  and  a  crude  estimation  of  tissue  elasticity  is 
obtained  through  the  sense  of  touch.  The  force  is  applied  on  the  body  surface  and  the  result  is  a 
collective  response  of  all  the  tissues  below.  Clinicians  can  sense  abnormalities  if  the  response  to 
palpation  of  the  suspicious  tissue  is  sufficiently  different  from  that  of  normal  tissue.  However, 
if  the  abnormality  lies  deep  in  the  body,  or  if  it  is  too  small  to  be  resolved  by  touch,  then 
the  palpation  method  fails.  The  dynamic  response  of  soft  tissue  to  a  force  is  also  valuable  in 
medical  diagnosis.  For  instance,  rebound  of  tissue  upon  sudden  release  of  pressure  exerted  by 
the  physician’s  finger  on  the  skin  provides  useful  diagnostic  information  about  the  tissue. 

Quantitative  measurement  of  the  mechanical  properties  of  tissues  and  their  display  in  a  raster 
format  is  the  aim  of  a  class  of  techniques  generally  called  elasticity  imaging  or  elastography  [1]. 
The  general  approach  is  to  measure  tissue  motion  caused  by  an  external  (or,  in  some  methods, 
internal)  force  or  displacement  and  use  it  to  reconstruct  the  elastic  parameters  of  the  tissue. 
The  excitation  stress  can  be  either  static,  or  dynamic  (vibration).  The  dynamic  excitation  is  of 
particular  interest  because  it  provides  more  comprehensive  information  about  tissue  properties  in 
a  spectrum  of  frequencies,  or  alternatively,  the  transient  behavior  of  the  tissue  could  be  deduced 
from  the  measurements.  In  most  elasticity  imaging  methods  ultrasound  is  used  to  detect  the 
motion  or  displacement  resulting  from  the  applied  stress.  Magnetic  resonance  elastography  is  a 
recently  developed  method  [2]  that  employs  a  mechanical  actuator  to  vibrate  the  body  surface 
and  then  measures  the  strain  waves  with  phase  sensitive  Magnetic  Resonance  Imaging  (MRI). 
The  majority  of  elasticity  imaging  methods  are  based  on  an  external  source  of  force,  resulting 
in  a  spatially  wide  stress  field  distribution.  This  requires  the  stress  field  to  pass  through  the 
superficial  portion  of  the  object  before  reaching  the  region  of  interest  within  the  interior.  This 
requirement  can  complicate  the  estimation  of  stiffness  because  the  stress  field  patterns  change  at 
different  depths.  Also,  because  the  stress  field  is  widely  distributed  the  response  of  the  object  is 
an  integral  sum  of  the  stress  field.  An  alternative  strategy  is  to  apply  a  localized  stress  directly  in 
the  region  of  interest.  One  way  to  accomplish  this  is  to  use  the  radiation  pressure  of  ultrasound. 

Acoustic  radiation  force  is  the  time  average  force  exerted  by  an  acoustic  field  on  an  object. 
This  force  is  produced  by  a  change  in  the  energy  density  of  an  incident  acoustic  field  [3],  for 
example,  due  to  absorption  or  reflection.  Several  benefits  may  result  from  using  ultrasound 
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radiation  force  for  evaluating  tissue  properties,  including:  a)  acoustic  (ultrasound)  energy  is  a 
non  invasive  means  of  exerting  force,  b)  existing  ultrasound  technology  and  devices  can  be 
readily  modified  for  this  purpose,  thus  eliminating  the  need  for  developing  a  new  technology,  c) 
radiation  force  can  be  generated  remotely  inside  tissue  without  disturbing  its  superficial  layers, 
d)  the  radiation  stress  field  can  be  highly  localized,  thus  allowing  for  precise  positioning  of 
the  excitation  point,  and  e)  radiation  force  can  be  produced  in  a  wide  range  of  frequencies  or 
temporal  shapes.  These  features  make  radiation  force  methods  more  attractive  than  other,  mostly 
mechanical,  excitation  methods  used  in  elasticity  imaging. 

Tissue  probing  with  the  radiation  force  produced  by  ultrasound  can  be  accomplished  in 
a  variety  of  techniques,  depending  on  the  excitation  and  detection  methods  used.  We  may 
categorize  these  method  as:  a)  transient  methods,  where  an  impulsive  radiation  force  is  used 
and  the  transient  response  of  the  tissue  is  detected  by  Doppler  ultrasound  [4];  b)  shear-wave 
methods,  where  an  oscillating  radiation  is  applied  to  the  tissue  and  the  resulting  shear  wave  is 
detected  by  ultrasound  or  other  methods  [5,  6,  7];  c)  vibro-acoustography,  recently  developed 
by  the  authors,  where  a  localized  oscillating  radiation  force  is  applied  to  the  tissue  and  the 
acoustic  response  of  the  tissue  is  detected  by  a  hydrophone  [15,  16].  Figure  [1]  illustrates  the 
relationship  between  various  methods  that  have  been  developed  for  evaluating  or  imaging  the 
elastic  properties  of  tissue. 

Dynamic  radiation  force  methods  seem  to  be  evolving  rapidly  as  a  new  field  in  tissue 
characterization  and  imaging.  The  purpose  of  this  chapter  is  to  systematically  discuss  the  features 
and  capabilities  that  are  offered  in  this  group  of  elasticity  measurement  techniques.  Here, 
we  describe  the  general  approach  used  in  the  transient,  shear  wave,  and  vibro-acoustography 
methods.  We  pay  particular  attention  to  the  last  method,  present  its  theory,  and  discuss  its 
features.  Finally  we  discuss  the  capabilities  and  limitations  of  all  three  methods. 

2.  Theory  of  the  Radiation  Force 

Study  of  radiation  force  and  radiation  pressure  dates  back  nearly  one  century,  to  the  time  of 
Rayleigh  [8].  Since  then,  this  subject  has  been  under  continuous  investigation.  A  historical 
review  of  radiation  force  and  radiation  pressure  is  presented  in  [9],  and  a  critical  review  of  the 
subject  can  be  found  in  [3],  Some  recent  analyses  of  radiation  force/pressure  in  attenuating 
medium,  which  may  be  applicable  to  biological  tissues,  are  presented  in  [10]  and  [11]. 

The  acoustic  radiation  force  is  an  example  of  a  universal  phenomenon  in  any  wave  motion 
that  introduces  some  type  of  unidirectional  force  on  absorbing  or  reflecting  targets  in  the  wave 
path.  Radiation  force  in  fluids  is  often  studied  in  the  context  of  radiation  pressure.  Depending 
on  the  boundary  conditions,  radiation  pressure  can  be  defined  differently.  Simple  explanations 
of  these  definitions  can  be  presented  by  considering  a  sound  traveling  inside,  and  along  the 
axis  of,  a  cylindrical  container  toward  the  opposite  wall  [9].  Rayleigh  radiation  pressure  is  the 
excess  pressure  produced  on  the  opposite  wall  when  the  container’s  side  wall  is  confining  the 
fluid  inside.  Langevin  radiation  pressure  is  the  excess  pressure  on  the  opposite  wall  when  we 
remove  the  confining  side  wall,  so  that  the  fluid  is  free  to  move  (fig.  2).  Here  we  will  focus 
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on  Langevin  radiation  pressure  because  the  conditions  for  which  this  pressure  is  defined  apply 
to  our  experimental  situation.  It  can  be  shown  that  the  Langevin  radiation  pressure  of  a  plane 
wave  impinging  normally  on  a  perfectly  absorbing  wall  is  equal  to  the  total  energy  density  ( E ), 
where  ()  represents  the  time  average.  If  the  wall  is  partially  reflecting,  this  pressure  would  be 
equal  to  (1  +  R){E),  where  R  is  the  power  reflection  coefficient  [9].  Thus,  in  general,  we  can 
write  the  radiation  force  of  a  normally  impinging  sound  beam  on  a  wall  as 

F  =  drS(E),  (1) 

where  F  is  the  force  in  the  beam  direction,  and  dr  is  the  “radiation  force  function”  or  the  drag 
coefficient.  This  dimensionless  coefficient  is  defined  per  unit  incident  energy  density  and  unit 
projected  area.  For  a  planar  object,  dr  is  numerically  equal  to  the  force  on  the  object.  Physically, 
the  drag  coefficient  represents  the  scattering  and  absorbing  properties  of  the  object  [12].  For  a 
perfectly  absorbing  object  dr  =  1,  and  for  a  perfectly  reflecting  object  dr  =  2.  In  the  case  of 
oblique  incidence,  the  radiation  force  will  have  a  normal  as  well  as  a  transverse  component.  A 
more  detailed  description  of  d.r  is  presented  in  [12]. 

To  produce  a  time  varying  radiation  force,  the  intensity  of  the  incident  beam  can  be  modulated 
in  various  ways.  For  example,  a  short  ultrasound  pulse  can  produce  a  transient  pulsed  radiation 
force,  and  a  sinusoidally  modulated  beam  can  result  in  a  sinusoidally  varying  force. 


3.  Radiation  Force  Methods 
3.1  Transient  Method 

In  transient  methods  the  radiation  force  of  ultrasound  is  used  to  make  a  minute  deformation  in 
the  tissue.  The  transient  recoil  of  the  tissue  resulting  from  this  deformation  is  measured  and  used 
for  evaluation  of  tissue  elastic  properties.  A  method  for  measuring  tissue  hardness,  presented 
by  Sugimoto,  et  al.  [4],  uses  the  radiation  force  of  a  single  focused  ultrasound  beam.  Ideally, 
hardness  may  be  represented  by  the  spring  constant  of  the  object,  which  is  the  ratio  of  the  applied 
force  to  the  displacement.  Principles  of  the  transient  methods  are  presented  in  fig.  3. 

In  this  method  an  ultrasound  pulse  is  used  to  generate  a  short  duration  radiation  force  which 
produces  localized  deformation  of  the  tissue.  Immediately  after  the  force  pulse,  the  resulting 
transient  deformation  of  tissue  is  measured  as  a  function  of  time  with  Doppler  ultrasound  using  a 
separate  transducer.  The  deformation  includes  an  initial  rapid  squeeze  of  the  tissue,  followed  by  a 
relaxation  and  possibly  a  rebound.  The  deformation  is  a  function  of  tissue  viscoelastic  parameters, 
as  well  as  the  applied  force.  The  authors  argue  that  because  quantitatively  measuring  the  internal 
radiation  force  is  difficult,  it  is  advantageous  to  derive  a  relative  quantity  that  is  representative  of 
tissue  hardness.  To  derive  a  single  relative  quantity  from  the  deformation  data,  the  relaxation  part 
of  the  function  is  approximated  by  a  sum  of  several  exponential  curves,  and  the  sum  of  the  first 
order  derivatives  of  such  exponential  is  calculated.  Sugimoto,  et  al.  [4]  show  that  this  quantity  is 
correlated  to  the  spring  constant  of  the  tissue,  thus  may  be  used  as  a  measure  of  tissue  hardness. 
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3.2  Shear  Wave  Methods 

Shear  modulus  is  related  to  the  hardness  or  elasticity  of  the  material.  It  is  known  that  the  shear 
modulus  of  various  soft  tissues  range  over  several  magnitudes  of  order,  while  bulk  modulus, 
a  parameter  that  is  associated  with  the  conventional  pulse-echo  ultrasound  compressional  wave 
speed,  varies  significantly  less  than  an  order  of  magnitude  [13,  14].  These  features  indicate  that 
shear  modulus  may  be  a  better  parameter  for  tissue  characterization  than  bulk  modulus.  Tissue 
attenuation  of  the  shear  wave  is  very  large  even  at  low  kHz  frequencies.  One  way  to  induce 
localized  shear  waves  inside  tissue  is  to  use  the  radiation  force  of  focused  ultrasound  [5].  In  a 
method  called  shear  wave  elasticity  imaging  (SWEI)  [6],  an  amplitude  modulated  single  focused 
ultrasound  beam  is  used  to  induce  a  localized  radiation  stress  inside  the  soft  tissue.  Localization 
of  the  stress  field  is  critical  to  success  of  the  method.  To  achieve  a  high  degree  of  localization, 
the  method  uses  a  focused  ultrasound  beam.  It  is  shown  that  the  radiation  stress  exerted  within 
a  dissipative  medium  peaks  about  the  focal  region  of  the  highly  focused  transducer.  Also,  it  has 
been  suggested  that  localization  can  be  improved  by  designing  the  transducer  and  selecting  the 
beam  parameters  such  that  a  nonlinear  shock  wave  is  produced  in  the  focal  region,  increasing 
the  magnitude  of  the  stress  field  in  the  vicinity  of  the  focal  region,  thus  augmenting  localization 
of  the  stress  field. 

Modulation  of  the  ultrasound  beam  can  be  in  the  form  of  an  oscillating  wave  or  short  pulse. 
The  resulting  radiation  force  elicits  a  shear  wave  propagating  in  the  radial  direction  with  respect 
to  the  beam  axis,  with  particle  motion  parallel  to  the  beam  axis.  Shear  waves  in  soft  tissue  travel 
at  very  low  speed,  typically  around  a  few  meters  per  second,  thus  the  corresponding  wavelength 
is  much  shorter  than  that  of  the  compressional  waves  for  the  same  frequency.  Shear  waves 
are  also  highly  attenuated  in  soft  tissue,  with  an  attenuation  coefficient  two  or  three  orders  of 
magnitude  higher  than  that  of  the  compressional  waves.  Because  of  high  attenuation  of  shear 
waves,  it  is  possible  to  induce  them  in  a  very  limited  region  in  the  vicinity  of  the  focal  point 
of  the  ultrasound  beam,  hence  avoiding  the  influence  of  tissue  boundaries.  Shear  parameters 
of  the  tissue,  such  as  shear  modulus  and  shear  viscosity,  can  be  calculated  by  measuring  the 
amplitude  and  the  temporal  characteristics  of  this  wave.  For  example,  the  time  required  for  the 
wave  front  to  propagate  from  one  point  to  another  can  be  used  to  calculate  the  shear  wave  speed, 
and  consequently,  the  shear  wave  modulus  //.,  as 

/'■  =  p4  (2) 

Where  p  and  ct  are  the  density  and  the  shear  wave  velocity,  respectively. 

Shear  waves  may  be  detected  optically.  In  this  method  a  laser  source  and  a  photo  detector 
are  used  to  detect  the  displacement  of  particles  due  to  the  shear  wave  in  a  transparent  phantom 
[6].  Because  this  method  requires  a  transparent  medium,  its  application  in  vivo  is  difficult.  Phase 
sensitive  magnetic  resonance  imaging  [2,  6]  is  an  alternative  method  that  can  be  used  to  measure 
the  three-dimensional  distribution  of  particle  displacement  in  a  given  direction  versus  time  in  a 
material.  Figure  4  illustrates  a  simplified  system  for  shear  wave  elasticity  measurement  using 
magnetic  resonance  imaging.  In  an  experiment  presented  in  [6],  an  ultrasound  pulse  of  3.6  ms 
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duration  produced  by  a  70  mm  diameter  transducer  focused  at  100  mm  was  transmitted  within 
a  cylindrical  rubber  phantom.  The  displacement  was  measured  by  2.0  T  MRI  system  at  two 
different  times  after  the  acoustic  pulse  was  applied.  The  position  of  the  peak  displacement  at 
these  time  points  was  used  to  estimate  the  shear  velocity,  which  was  shown  to  be  consistent  with 
the  independently  measured  value.  Shear  waves  can  also  be  detected  by  Doppler  ultrasound  [5]. 
It  has  been  shown  [7]  that  to  achieve  an  appreciable  displacement  needed  for  Doppler  detection, 
most  soft  tissues  require  high  ultrasound  intensities  which  might  be  beyond  the  safe  limit. 

3.3  Vibro-acoustography 

3.3.1  Method 

This  technique  produces  a  map  of  the  mechanical  response  of  an  object  to  a  dynamic  force 
applied  at  each  point.  The  method  utilizes  ultrasound  radiation  force  to  remotely  exert  a  localized 
oscillating  stress  field  at  a  desired  frequency  within  (or  on  the  surface  of)  an  object,  and  records 
the  resultant  acoustic  response  [15,  16].  Figure  5  illustrates  the  principle  of  this  method.  This 
acoustic  response,  which  is  normally  at  low  kHz  range,  is  a  function  of  the  viscoelastic  properties 
of  the  object  and  can  be  used  to  produces  an  image  of  the  object. 

This  method  operates  on  an  oscillating  radiation  force  to  vibrate  the  tissue  and  consequently 
produce  an  acoustic  emission  from  the  object.  To  produce  an  oscillating  radiation  force  the 
intensity  of  the  incident  ultrasound  must  be  amplitude  modulated  at  the  desired  low  frequency. 
Using  a  single  amplitude  modulated  beam  seems  to  be  the  simplest  means  to  attain  this  purpose. 
However,  such  a  beam  could  exert  a  radiation  force  on  any  object  that  is  present  along  the  beam 
path,  producing  undesirable  acoustic  emission.  To  confine  the  radiation  stress  to  the  desired 
region,  we  use  two  unmodulated  CW  beams  at  slightly  different  frequency,  propagating  along 
separate  paths.  The  beams  are  positioned  to  cross  each  other  at  their  respective  foci,  and  thus 
produce  a  modulated  field  at  a  confined,  small,  cross-section  region. 

3.3.2  Theory 

In  the  beam  forming  method  to  be  described  the  amplitude  modulated  field  is  obtained  by  the 
interference  of  two  unmodulated  ultrasound  beams.  The  main  advantage  of  this  approach  is  that 
the  interference  volume  size  can  be  limited  to  a  small  region.  Radially  symmetric  interfering 
beams  are  obtained  when  two  coaxial,  confocal  transducers  are  used  [16].  For  this  purpose, 
elements  of  a  two-element  spherically  focused  annular  array  (consisting  of  a  central  disc  with 
radius  a\  and  an  outer  ring  with  the  inner  radius  of  a'2  and  outer  radius  of  0,2)  are  excited  by 
separate  CW  signals  at  frequencies  wi  =  wo  —  Aw/ 2  and  W2  =  wo  +  Aw/2,  respectively.  We 
assume  that  the  beams  are  propagating  in  the  +z  direction  with  the  joint  focal  point  at  2  =  0, 
as  shown  in  fig.  6.  The  resultant  field  on  the  z  —  0  plane  may  be  written  as: 

p(t)  —  Pi(r)cos(wif  +  V'i(r))  + -P2M  cos(w2*  + '02(7’)),  (3) 

where  r  =  a/.t2  +  y2  is  the  radial  distance.  The  amplitude  functions  are  [17,  18] 

Pl(r)  =  ^°‘X^jinc(S)’  <4) 
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P2(r)  =  pct/02_^jmcv^„y  . . VA22#;j, 

where  U(u  is  the  velocity  amplitude  at  the  i-th  transducer  element  surface,  and  A;  =  2-Tr/o-v  for 
i=7,2  are  the  ultrasound  wavelengths.  The  phase  functions, 


a2jmc 


Mr)  = 


A;V 


for  i=7,  2,  are  conveniently  set  to  be  zero  at  the  origin.  Also,  jinc(X)  =  Ji(2%X)/-kX,  where 
Ji(.)  is  the  first  order  Bessel  function  of  the  first  kind. 

The  instantaneous  energy  is  E  =  p2(t)/pc2.  Replacing  p(t)  from  (3),  this  energy  will  have 
a  time  independent  component,  a  component  at  the  difference  frequency  Aw  =  u>2  -  o/i  which 
results  from  the  cross  product  of  the  two  pressure  fields,  and  high  frequency  components  at  uj\ 
and  u)2  and  their  harmonics.  The  energy  component  at  the  difference  frequency  is: 

eA oj(t)  =  cos  [Ao>f  -  A tf(ro)],  (7) 

where  AV’(ro)  =  ip2(r o)  —  V’i(ro)- 

Now,  we  define  a  unit  point  target  with  an  area  of  dxdy  at  position  (xo,yo)  on  the  focal 
plane,  and  with  a  drag  coefficient  dr(x o,  yo)  such  that:  dr(xo,  yo)dxdy  =  1  on  the  target  and  zero 
elsewhere.  This  equation  is  merely  used  as  a  mathematical  model.  In  this  case,  the  projected  area 
can  be  considered  to  be  S  =  dxdy.  Therefore,  if  the  projected  area  is  unity,  then  dr(x,y)  =  1, 
which  corresponds  to  a  totally  absorptive  object. 

Referring  to  (1),  and  replacing  drS  with  unity  and  (E)  by  e^(t)  of  (7),  then  we  can  write 
the  low-frequency  component  of  the  radiation  force  on  the  unit  point  target  as: 

/aw(^0,  yo;  t)  =  -^2  pi(ro)^2(ro)  cos  [A uit  +  A^(ro)],  (8) 

pc 

where  arguments  jc o  and  yo  are  added  to  denote  the  position  of  the  point  target,  and  tq  = 
yfx  q  +  Pq.  Referring  to  (4)  and  (5),  the  complex  amplitude  of  the  stress  field  can  be  found  as 


FAu(x0,y0)  =pUoiUo2~. — ^jinc 

Ai^o 


2  AA^o  J  ’ 

where  AA  =  2ttc/Auj  is  the  wavelength  associated  with  Aw.  The  above  equation  indicates  that 
the  radiation  stress  is  concentrated  at  the  focal  point  and  decays  quickly  as  jinc(.)2.  It  should  be 
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noted  that  because  /aw(-^0,  2/0 ;  t)  and  FAoJ(x0,  y0)  are  defined  per  unit  point  target,  they  have  the 
dimension  of  force  per  unit  area,  or  equivalently  the  unit  of  stress.  As  an  example,  we  consider 
a  confocal  transducer  with  dimensions  a\  =  14.8  mm,  a2  =  16.8  mm,  «2  =  22.5  mm,  and  the 
focal  length  70  mm.  Also,  we  assume  that  the  center  frequency  is  3  MHz,  and  the  difference 
frequency  is  7.3  kHz.  The  radiation  stress  at  the  focal  plane  of  this  transducer  is  plotted  in  fig.  7. 

In  medical  applications  the  maximum  ultrasonic  intensity  is  regulated  for  safety  reasons.  It 
is  therefore  useful  to  write  the  stress  field  in  terms  of  the  peak  ultrasonic  intensity  at  the  focal 
point.  The  long-term  average  of  ultrasonic  intensity  at  the  focal  point  can  be  written  as 

/(0)  =  Pi2(°)  +  (10) 

2  pc 

where  Pi(0)  and  7^(0)  can  be  found  from  (4)  and  (5)  as 


m  = 


Assuming  Uqi 


U0,  we  can  write  the  focal  plane  stress  field  in  terms  of  1(0): 


-Fa»(*o,  Vo) 


27(0) 


(£j)2+fe)2(»J 


"  V  '2 AAV 

The  first  fraction  in  the  above  equation  represents  the  static  radiation  force  produced  by  the  two 
beams  on  a  total  absorber,  the  second  fraction  is  a  constant  factor,  and  the  rest  represents  the 
spatial  distribution  of  the  stress  field  on  the  focal  plane.  If  Aw  <  W2,  then  we  may  replace 
Ai  and  A2  by  with  Ao,  and  simplify  the  expression  for  the  stress  field.  Under  these  conditions, 
the  stress  at  the  focal  point  is 


P Au;(0,  0)  — 


2/(0)  x  ai(a2  ~a22)  _ 

c  af  +  (a2  —  a, 2  ) 


The  fraction  on  the  right  represents  the  effect  of  transducer  dimension  on  the  stress  field.  For 
the  transducer  used  in  the  previous  example,  the  value  of  this  fraction  is  0.4999.  (This  fraction 
can  be  also  calculated  for  a  single  element  transducer  of  the  same  diameter  excited  by  an 
amplitude  modulated  signal.  For  this  purpose,  we  may  let  a\  =  a 2  and  a'2  =  0.  For  these  values, 
the  last  fraction  in  (13)  is  0.5.)  The  resulting  radiation  stress  (assuming  c.  =  1500  m/sec  for 
water)  is  Fau>(0,0)  =  6.67  x  10-4/(0)  N/m2.  Now  letting  7(0)  =  7200  W/m2,  which  is  the 
intensity  suggested  by  FDA  for  safe  in  vivo  applications,  the  resulting  stress  at  the  focal  point 
is  FAoj(0,  0)  =  4.80  N/m2.  Referring  to  fig.  7,  we  note  that  this  stress  field  is  applied  only  in  a 
small  region  around  the  focal  point  to  the  object. 
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3.3.2. 1  Acoustic  Emission 

To  explain  the  acoustic  emission  we  consider  an  “object”  within  a  homogeneous  infinite 
medium.  This  model  allows  us  to  separate  the  roles  played  by  the  parameters  of  the  object 
and  the  surrounding  medium.  Also,  this  model  can  be  fitted  to  various  applications.  When  an 
oscillating  stress  field  is  applied  to  the  object,  the  object  vibrates  at  the  frequency  of  the  stress 
field.  Vibrational  energy  of  the  object  is  partly  transferred  to  the  surrounding  medium,  resulting 
in  an  acoustic  emission  field.  Here,  we  calculate  the  acoustic  emission  of  an  object  subjected 
to  cyclic  radiation  stress. 

To  explain  the  physics  in  an  analytical  form,  we  consider  a  flat  plate  facing  the  beam. 
Here,  we  assume  that  the  vibrating  object  has  a  circular  cross  section  of  radius  b  and  uniformly 
vibrates  back  and  forth  like  a  piston.  We  also  consider  an  area  S  <  tx b2  of  the  piston  surface 
to  be  projected  normally  to  the  beam.  We  can  always  return  to  our  elementary  point  object 
by  reducing  the  area  of  this  disk  to  dxdy.  Similar  solutions  can  be  carried  out  for  objects  of 
other  forms.  The  theory  can  be  also  extended  to  include  arbitrary  vibrating-part  shapes  and  non 
uniform  displacement  of  the  object.  The  total  radiation  force  on  this  object  Fau>,  can  be  found 
by  integrating  the  radiation  stress  over  the  area  of  the  object.  This  force  vibrates  the  target  object 
at  frequency  Aw.  The  steady  state  normal  velocity  amplitude  of  a  piston  at  frequency  Aw,  Uaw, 
due  to  a  harmonic  force  Faw,  can  be  written  as 


U Ato 


Faoj 
Zau  ’ 


(14) 


where  Zau  is  the  mechanical  impedance  of  the  object  at  Aw.  The  mechanical  impedance  of 
the  object  has  two  components,  one  resulting  from  the  inertia,  friction,  and  the  elasticity  of  the 
object  itself,  and  the  other  resulting  from  the  loading  effect  of  the  surrounding  medium  on  the 
vibrating  object.  The  mechanical  impedance  can  be  interpreted  as  a  measure  of  object  rigidity 
and  how  much  it  yields  to  the  applied  force.  For  example,  for  a  rigid  object,  Zaw  is  high,  and 
hence  resisting  the  force. 

Knowing  Uau,  we  can  calculate  the  pressure  field  it  produces  in  the  medium.  We  assume 
that  the  acoustic  emission  signal  propagates  in  a  free  and  homogenous  medium.  The  farfield 
acoustic  pressure  due  to  a  piston  source  of  radius  b  set  in  a  planar  boundary  of  infinite  extent 
is  given  by  [17], 


Palo  =  -jAu/p 


exp  (jAu)l/c) 
47 tI 


2  J\  (  sin  -i?)  cos  d 


A  wb 


sin  d 


cos  d  4-  0b 


(2'k}?UaJ)-,  (15) 


where  l  is  the  distance  from  the  observation  point  to  the  center  of  the  piston,  d  is  the  angle 
between  this  line  and  the  piston  axis,  and  (3b  is  the  specific  acoustic  admittance  of  the  boundary 
surface  (the  specific  acoustic  admittance  is  (3b  =  where  Zb,  the  acoustic  impedance  of  the 
boundary,  represents  the  ratio  between  the  pressure  and  normal  fluid  velocity  at  a  point  on  the 
object).  The  factor  two  comes  from  the  presence  of  the  boundary  wall.  It  would  be  replaced 
by  unity  if  the  boundary  wall  were  not  present  [17]. 
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The  acoustic  emission  field  resulting  from  object  vibration  can  be  written  in  terms  of  object 
mechanical  impedance  by  combining  (14)  and  (15),  as 


PA„  =  ^^xeXPyA“'/C) 


47 rZ 


2</i(^psin7?)  ^  cost) 


A  tvb 
c 


sin 


x 


cost?  +  Pb 


x 


27 r62 


(16) 


For  wavelengths  long  compared  to  the  object  size,  i.e.,  when  bAto/c  — >  0,  the  term  in  the 
bracket  approaches  a  constant,  hence  we  may  consider  the  contents  of  the  brace  to  be  an  object 
independent  function  (the  specific  acoustic  admittance  Pb  relates  to  the  surrounding  boundary 
surface).  Under  these  conditions,  the  brace  in  the  above  equation  represents  the  effect  of  the 
medium  on  the  acoustic  emission  field,  which  we  may  call  the  medium  transfer  function,  and 
denote  it  by 


tt  n\  ,Aw  w  exp(jAul/c)  cost? 

HU‘)  =J—x - 4 ~t - x  — - 


(17) 


The  last  fraction  in  (16)  includes  which  is  the  mechanical  admittance  of  the  object  at 
the  frequency  of  the  acoustic  emission  (Aw),  and  we  denote  this  by  Yaw  It  is  convenient  to 
combine  this  term  with  the  next  term  (2-Kb2)  in  (16),  as  Qau  =  27r 62Yaw  =  JZ^,  which 
is  the  total  acoustic  outflow  by  the  object  per  unit  force  (acoustic  outflow  is  the  volume  of  the 
medium,  e.g.,  the  fluid,  in  front  of  the  object  surface  that  is  displaced  per  unit  time  due  to  object 
motion).  Function  Qau>  represents  the  object  characteristics  at  the  acoustic  frequency.  We  may 
thus  rewrite  (16)  in  a  more  compact  form  as 


P/Sw  =  Pc2Hau>(1)Q/SujFaw  (18) 

Equation  18  indicates  that  the  acoustic  emission  pressure  is  proportional  to:  a)  the  radiation  force, 
which  itself  is  proportional  to  the  square  of  ultrasound  pressure  and  the  ultrasound  characteristics 
of  the  object,  dr,  in  the  projected  area  5;  b)  the  acoustic  outflow  by  this  object,  Qalo,  representing 
the  object  size  b  and  its  mechanical  admittance  at  the  acoustic  frequency,  Ya ad  and  c)  the  transfer 
function  of  the  medium  at  the  acoustic  frequency,  J/aw(0-  Note  the  difference  between  the 
projection  area  S  and  the  vibrating  area  nb2.  The  projection  area  determines  the  extent  of  the 
force  applied  to  the  object  (1).  The  vibrating  area,  however,  influences  the  total  acoustic  outflow 
in  the  medium  caused  by  object  vibration. 

3.3.3  Applications 

Evaluating  the  characteristics  of  an  object  (or  a  medium)  by  listening  to  its  sound  is  a  traditional 
approach  that  has  been  used  for  many  purposes.  Qualitative  evaluation  of  a  crystal  glass  by 
tapping  on  it  is  a  simple  example.  Vibro-acoustography  implements  the  same  approach  but  in 
micro-scale,  and  in  a  way  that  could  be  applied  to  tissues. 
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Equation  18  presents  a  general  relationship  between  the  mechanical  parameters  of  the 
objects,  surrounding  medium,  and  the  acoustic  emission  field  resulting  from  object  vibration. 
By  measuring  the  acoustic  emission  field,  it  would  be  possible,  in  principle,  to  estimate  some  of 
the  object  or  medium  parameters,  either  in  absolute  or  in  relative  sense.  For  example,  one  can  use 
vibro-acoustography  to  measure  the  resonance  frequency  of  an  object,  and  from  that  information 
it  is  possible  to  estimate  some  viscoelastic  parameters  of  the  object  or  the  medium.  This  method 
has  been  used  to  measure  the  Young’s  modulus  of  a  metallic  rod  [19].  In  another  experiment, 
by  measuring  the  resonance  frequency  of  a  known  resonator  in  a  liquid  medium,  the  viscosity 
of  the  fluid  has  been  estimated  with  good  accuracy  [20].  These  applications  are  not  necessarily 
medical,  but  the  principle  could  be  used  for  evaluation  of  soft  tissues,  blood,  bone,  etc. 

In  medical  applications,  one  can  use  vibro-acoustography  to  obtain  images  of  human  body 
for  diagnostic  purposes.  In  such  applications,  the  image  may  not  represent  a  physical  quantity, 
rather  it  provides  a  means  to  visualize  object  details.  To  demonstrate  the  capability  of  vibro- 
acoustography,  we  used  this  method  to  image  a  specimen  of  human  iliac  artery  that  included 
calcifications.  A  3MHz  confocal  transducer  with  dimensions  a i  =  14.8  mm,  a'2  =  16.8  mm, 
«2  =  22.5  mm,  and  the  focal  length  of  70  mm,  was  used  for  this  experiment.  The  difference 
frequency  was  set  at  44  kHz.  Figure  9  shows  a  photo  of  the  artery  that  is  cut  open.  The  areas 
with  calcification  are  marked  by  letter  “C”.  The  vibro-acoustic  image  of  the  same  artery  is  shown 
on  the  right.  Calcium  deposits,  which  are  formed  in  large  (a  few  mm  to  a  few  cm  long)  and  stiff 
(compare  to  the  arterial  wall)  plates,  are  efficient  acoustic  radiators,  producing  strong  acoustic 
emission  when  they  are  exposed  to  the  radiation  force  of  the  incident  ultrasound.  Hence,  they 
stand  out  in  the  vibro-acoustography  image.  The  vibro-acoustic  image  has  high  resolution,  no 
speckle,  and  high  signal-to-noise  ratio.  These  features  allow  vibro-acoustography  to  delineate 
calcium  deposits  with  high  definition  and  contrast. 

4.  Capabilities  and  Limitations 

Although  the  methods  described  above  tend  to  estimate  the  mechanical  properties  of  the  object, 
the  quantities  they  measure  are  different.  The  quantity  measured  in  the  transient  method  is 
somewhat  arbitrary,  and  not  a  direct  measure  of  any  physical  quantity.  The  merit  of  this  method 
is  that  tissue  hardness  is  measured  from  relative  displacement  values,  without  the  need  to  measure 
the  applied  force.  The  shear  wave  method  directly  measures  a  physical  parameter  of  the  material 
in  an  absolute  sense.  The  acoustic  emission  field  measured  in  vibro-acoustography  is  a  function 
of  several  physical  parameters  that  relate  to  the  object  and  the  surrounding  medium.  It  is  possible 
to  measure  a  single  parameter  only  if  one  has  enough  knowledge  about  other  parameters  involved. 
For  example,  one  can  measure  the  shear  viscosity  of  a  liquid  using  a  known  resonator  as  an 
object  in  the  liquid  [20],  or,  for  a  given  the  geometry,  the  Young’s  modulus  of  a  metallic  object 
can  be  measured  [19]  by  vibro-acoustography. 

Detection  sensitivity  is  critical  to  the  success  of  any  of  the  above  methods.  Greater  ultrasound 
intensity  would  be  needed  if  the  detection  sensitivity  is  poor.  As  a  result,  the  peak  ultrasound 
intensity  at  focal  area  may  be  limited  by  the  safe  in  vivo  limit.  Ultrasound  Doppler  methods, 
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at  conventional  ultrasound  frequencies,  probably  can  detect  displacements  in  the  order  of  a  few 
micrometers.  The  sensitivity  of  the  MRI  technique  is  in  the  order  of  100  nm  [2].  Vibro- 
acoustography  has  been  shown  to  detect  motions  as  small  as  a  few  nanometers  [15].  The  high 
sensitivity  of  this  method  is  a  result  of  the  fact  that  small  motions  of  the  object  can  produce  an 
acoustic  emission  pressure  field  that  is  easily  detectable  by  a  sensitive  hydrophone. 

Image  spatial  resolution  is  an  important  parameter  when  assessing  the  capability  of  an 
imaging  method  in  delineating  fine  structures  in  the  object.  Vibro-acoustography  can  be  used 
for  high  quality  imaging.  The  spatial  resolution  is  proportional  to  the  width  of  the  stress  field 
main  lobe.  For  example,  for  the  3  MHz  ultrasound  transducer  presented  in  fig.  7,  the  spatial 
resolution  is  about  700  micrometers.  Imaging  complex  objects  using  the  transient  or  the  shear 
wave  methods  have  not  been  fully  explored  in  the  literature. 

Performance  of  each  method  is  partly  influenced  by  viscoelastic  properties  of  the  material. 
Both  the  transient  and  shear  wave  methods  work  best  if  the  material  under  test  supports  the  shear 
wave,  has  enough  attenuation  to  allow  the  build  up  of  enough  radiation  force,  and  is  compliant 
enough  to  allow  appreciable  displacement.  Application  of  these  methods  for  hard  material,  such 
as  bone  and  calcifications,  would  be  difficult  or  impractical  because:  a)  stiff  materials  have  low 
compliance,  hence  their  displacement  in  response  to  the  force  is  relatively  small  and  difficult 
to  detect;  b)  the  shear  wave  travels  at  high  speed  in  such  material.  In  vibro-acoustography 
hard  materials  produce  structural  vibrations  which  are  often  stronger  that  of  soft  tissue.  Such 
vibrations  result  in  strong  and  easily  detectable  acoustic  emissions.  The  bright  appearance  of  the 
arterial  calcifications  in  the  experiment  presented  in  [15]  is  a  demonstration  of  this  phenomenon. 
Vibro-acoustography  can  be  used  to  detect  particles  in  materials  that  do  not  support  share  waves, 
for  example,  detecting  gas  bubbles  in  liquids.  The  acoustic  emission  resulting  from  particle 
vibration  includes  a  compressional  wave  component  that  can  travel  in  the  surrounding  medium 
including  the  liquid.  Tissue  attenuation  for  compressional  waves  at  low  frequencies  is  small, 
hence  such  waves  can  be  detected  by  the  hydrophone  from  a  distance.  Shear  waves  in  liquid 
medium  decay  very  rapidly  and  are  difficult  to  detect  by  most  methods. 

The  ultimate  goal  of  viscoelastic  parameter  imaging  methods  is  to  measure  these  parameters 
in  human  body.  In  vivo  application  of  the  transient  method  is  not  clear  because  it  requires 
high  ultrasound  power  to  produce  a  detectable  displacement.  Application  of  the  shear  wave 
method  in  human  body  largely  depends  on  the  detection  method  used.  Optical  methods  are 
not  likely  usable  in  vivo.  MRI  detection  methods  are  generally  complex,  especially  when  the 
technique  is  to  be  linked  to  the  ultrasound  system  in  the  magnetic  field.  Ultrasound  Doppler  is 
a  more  practical  detection  method  for  shear  wave  detection,  however  its  sensitivity  may  not  be 
sufficient  for  ultrasound  intensities  in  the  safe  limit  [7].  Vibro-acoustography  uses  a  hydrophone 
for  detection  which  is  simple  to  operate  in  clinical  settings.  It  however  requires  an  acoustically 
quiet  environment  for  proper  detection.  The  biological  noise  spectrum  seems  to  be  below  1  kHz, 
which  can  be  easily  filtered  out  if  the  operation  frequency  is  above  this  limit  [15]. 

Before  any  of  the  methods  discussed  here  are  used  in  the  in  human  body,  one  must  be  certain 
that  the  ultrasound  exposure  does  not  harm  the  subject.  Vibro-acoustography  uses  either  CW 
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or  tone  burst  ultrasound.  Therefore,  the  continuous  ultrasound  intensity  (spatial  peak  temporal 
average  intensity)  limit  must  be  observed  for  the  safety  reasons.  This  limit  according  to  the 
FDA  is  720  mW/cm2.  However,  because  of  the  high  sensitivity  of  hydrophone  detectors,  it  is 
possible  to  detect  very  small  levels  of  the  acoustic  emission  while  using  low  transmit  power. 
For  example,  it  is  shown  that  a  submillimeter  object  can  be  detected  at  ultrasound  intensities  far 
below  the  FDA  limit  [16].  The  shear  wave  method  presented  in  [6]  uses  high  intensity  pulsed 
ultrasound  for  excitation  and  either  the  optical  or  MRI  methods  for  detection.  It  is  argued  that, 
although  the  intensity  is  high,  the  total  exposure  is  within  the  FDA  limits.  However,  in  vivo 
use  of  such  detection  methods  would  be  difficult.  Using  Doppler  ultrasound  for  detection  is  less 
sensitive,  hence  requires  higher  peak  ultrasound  intensity  to  produce  enough  radiation  force  on 
the  tissue.  This  intensity  level  may  exceed  the  safe  spatial  peak  pulse  average  intensity  limit. 

5.  Summary 

Radiation  force  or  ultrasound  is  a  non  invasive  means  for  introducing  highly  localized  vibrating 
force  inside  tissue.  The  three  methods  described  here  utilize  this  tool  to  evaluate  the  dynamic 
viscoelastic  properties  of  tissue.  The  transient  method  measures  minute  transient  tissue  defor¬ 
mation  versus  time.  The  shear  wave  method  measures  the  amplitude  and  velocity  of  the  shear 
waves  resulting  from  excitation  by  radiation  force.  Finally,  Vibro-acoustography  measures  the 
acoustic  field  resulting  from  object  vibration  at  a  specified  frequency. 
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8.  Legends 


Figure  1  Tissue  elasticity  evaluation  and  imaging  methods.  These  methods  can  be  divided  into 
two  general  groups.  The  static  group  employs  a  steady  force  to  deform  the  object.  The  dynamic 
group  use  either  a  momentary  or  an  oscillating  force  to  vibrate  the  object  and  can  be  divided 
into  two  subgroups.  Mechanical  drive  methods  use  a  mechanical  actuator,  whereas  the  radiation 
force  drive  methods  use  the  radiation  force  of  ultrasound  to  excite  the  object.  The  latter  subgroup 
is  divided  into  three  sub-subgroups  based  on  the  detection  method  used.  Vibro-acoustography 
uses  a  hydrophone  to  detect  object  response. 


Figure  2  Radiation  force  produced  by  projection  of  ultrasound  on  an  object.  Interaction  of  an 
ultrasound  beam  with  an  object  that  scatters  and/or  absorbs  results  in  a  force  on  the  object  in 
the  beam  direction.  The  magnitude  of  this  force  is  proportional  to  the  time  average  energy 
density  of  the  incident  beam  and  a  factor  that  represents  the  scattering  and  absorbing  properties 
of  the  object. 


Figure  3  Principle  of  the  transient  method.  The  excitation  transducer  produces  a  momentary 
radiation  force  which  induces  a  minute  deformation  in  the  object.  Doppler  methods  are  used  to 
measure  this  deformation  using  the  probe  transducer. 


Figure  4  Principle  of  shear  wave  elastography.  The  focused  ultrasound  beam  is  absorbed  by  the 
object,  resulting  a  localized  radiation  stress,  which  squeezes  the  elastic  object.  This  deformation 
propagates  in  the  form  of  shear  waves  through  the  object.  The  phase  sensitive  MRI  machine 
detects  the  spatial  distribution  of  the  resulting  displacement  in  the  object. 


Figure  5  Principle  of  vibro-acoustography.  The  ultrasound  beam  is  used  to  produce  a  localized 
oscillatory  radiation  stress  to  vibrate  the  object.  Vibration  of  the  object  produces  an  acoustic 
emission  field  in  the  medium.  This  field,  which  is  a  function  of  the  viscoelastic  properties  of 
the  object,  is  detected  by  a  sensitive  hydrophone. 


Figure  6  The  coordinate  system  describing  the  position  of  the  ultrasound  beam  and  the 
observation  point  relative  to  the  object.  The  beam  travels  in  the  z-direction. 

Figure  7  Normalized  stress  field  of  the  vibro-acoustography  system  with  a  3  MHz  confocal 
transducer. 
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Figure  8  Vibro-acoustography  system.  The  object  is  vibrated  by  the  radiation  force  of  the 
confocal,  two-beam  transducer.  Resulting  vibrations  are  detected  by  the  microphone  (or 
hydrophone)  and  mapped  into  the  image. 


Figure  9  Vibro-acoustic  image  of  an  ex-vivo  human  iliac  artery  near  the  bifurcation.  On  the 
left  is  a  photo  of  the  cut-open  artery  specimen.  Calcified  areas  are  marked  with  letter  “C”. 
Vibro-acoustic  image  of  the  specimen  is  shown  on  the  right.  Calcified  areas  stand  out  as  bright 
spots.  This  image  was  obtained  with  a  3  MHz  confocal  transducer  as  described  in  the  text.  The 
difference  frequency  was  set  at  44  kHz.  The  gray-level  scale  is  in  arbitrary  units. 
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Ultrasound-stimulated  vibro-acoustic  spectrography  is  a  recently-developed  method  that  employs 
the  radiation  force  of  two  intersecting  continuous  ultrasound  beams  to  remotely  vibrate  an  object  at  an 
arbitrary  low  frequency.  Object  vibration  produces  a  sound  field  (acoustic  emission)  in  the  medium, 
which  is  a  function  of  object  mechanical  properties.  By  measuring  the  acoustic  emission  field,  one  can 
obtain  information  about  the  mechanical  parameters  of  the  object.  In  this  paper,  we  use  this  method  for 
remote  (noncontact)  measurement  of  the  dynamic  Young’s  modulus  of  a  rod  based  on  its  fundamental 
resonance  frequency.  Experimental  results  on  an  aluminum  rod  agree  with  the  published  data. 

Key  Words:  Acoustic  emission;  elasticity;  noncontact  measurement;  radiation  force;  resonance;  ultra¬ 
sound-stimulated  vibro-acoustic  spectrography;  Young’s  modulus. 


INTRODUCTION 

Estimation  of  material  parameters  using  their  mechanical  response  to  externally-applied 
forces  is  of  considerable  interest  in  nondestructive  inspection  of  materials,  material  science, 
and  medical  diagnosis.  Normally,  an  applied  force  is  used  to  produce  displacement  from 
which  elastic  constants,  like  spring  constants,  can  be  determined.  In  medical  applications, 
ultrasound  has  been  used  to  measure  tissue  displacement  associated  with  externally  applied 
compressive  and  cyclic  forces.5  In  resonant  ultrasound  spectroscopy,  an  ultrasound  source 
and  a  detector  are  used  to  measure  the  resonance  frequencies  of  a  sample  with  known  size 
and  mass.  The  resonances  are  related  to  mechanical  parameters,  including  the  elastic  con¬ 
stants  of  the  material.2  In  some  applications,  it  is  important  to  measure  the  elastic  constants 
of  a  material  at  a  particular  frequency,  since  it  is  well  known  that  the  static  and  dynamic 
moduli  of  some  materials,  such  as  plastics,  can  differ  substantially  due  to  the  existence  of  re¬ 
laxation  time  effects  that  can  make  the  modulus  measured  by  the  conventional  quasi-static 
stress-strain  curves  significantly  lower  than  the  dynamic  modulus.  Acoustic  resonators  have 
been  found  to  be  valuable  tools  for  investigating  material  properties.  In  general,  resonance 
frequencies  are  functions  of  elastic  constants  and  density,  and  resonator  strength  (or  Q)  is  a 
function  of  the  viscous  loss  and  friction.  Suitable  measurement  of  the  frequency  response  of 
resonators  can  therefore  yield  a  wide  range  of  information  regarding  the  elastic  properties 
and  structure  of  materials.  To  study  a  material,  normally  a  sample  of  known  geometry,  such 
as  a  rod  or  a  tuning  fork,  is  used.  In  the  case  of  a  rod,  the  Y  oung’s  modulus  can  be  determined 
from  its  resonance  frequencies  in  various  modes,  such  as  flexural  or  longitudinal.  The  use  of 
structures  in  their  resonance  modes  insures  high  signal-to-noise  ratio.3 

Estimation  of  the  elastic  moduli  by  the  resonance  method  requires  some  form  of  actuator 
and  a  motion  detection  device,  such  as  piezoelectric  or  magnetic  transducers.4  It  has  been  re¬ 
al  ized  that  noncontact  methods  are  superior  to  contact  methods  because  contact  with  the  ob¬ 
ject  adds  damping  and  stress,  making  object  modeling  more  complex.3  Another  consideration 
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is  limitation  in  accessing  the  object  when  the  object  under  test  is  very  small  compared  to  the 

m^^t^33umulatedVibn)-Acoustic  Spectrography  (USVAS)  method5  ut.hzes  a 

highly-localized  radiation  force,6'7  of  two  intersecting  ultrasound  beams  tc .excite  the ^ob ec 
at  arbitrary  low  frequences  and  measures  the  frequency  response  of this  object  v.a  a  re^ 
microphone.  This  method  provides  high  spatial  definition  because  it  is  “Pab'eof^enm 
an  oscillatory  force  in  a  very  small  region  of  the  object.  Thrs  feature,  plus  the  that  object 
frequency  response  can  be  measured  remotely,  makes  USVAS  a  poten  ra  can  1 

material  characterization.  Q  .  H  fnr„ 

In  this  paper,  we  introduce  a  noncontact  technique,  based  on  the  USVAS  method,  fo 
motely  measuring  the  Young’s  modulus  of  material.  Here,  we  describe  howultrasoundcan 
be  used  to  measure  the  Young’s  modulus  of  a  metalic  rod  at  frequences  much  lower  hat  Ae 
ultrasound  frequency  employed.  This  method  allows  us  to  avoid  physical  contac  s  .n  bo 
the  excitation  and  detection  processes,  while  providing  the  convenience  of  accessing  smaU 
objects  and  avoiding  the  surrounding  structures  by  focusing  the  impinging  ultrasound  beam 
on  the  object.  Our  purpose  is  to  present  the  methodology.  Although  the  accuracy  oft 
method  can  be  similar  to  that  of  other  resonance  methods,  we  shall  not  focus  on  t  is ;  issue ^  er  ^ 
In  the  following  sections,  we  will  present  the  basic  theory,  method  and  some  experimental 

results  that  show  the  utility  of  the  method. 


METHOD  AND  MATERIALS 

Ultrasound-Stimulated  Vibro-Acoustic  Spectrography  (USVAS) 

When  an  acoustic  wave  is  absorbed  or  reflected  by  an  object,  causing  a  change >  in  the 
acoustic  energy  density,  a  steady  force,  which  is  called  the  radiation  force,  is  exerted  in  the 
direction  of  tte  beam  on  that  object.  For  a  normal  incident  beam,  the  radiation  force  is  given 

i  6,7 

by: 

(i) 

F=drW/c 

where  W  c  and  d  are  the  total  incident  power  on  the  object,  speed  of  sound  and  a  constant 

(the  drag  coefficient)  that  depends  on  the  absorption  and  reflection  properties  of  teojrc, 

respectively  The  force  decreases  for  incident  angle  less  than  90  .  For  a  total  re  ’  '  ,  ' 

For  a  total  absorber,  dr  =  1 .  When  an  amplitude  modulated  ultrasound  beam  impinges  on  die 
object  an  oscillating  radiation  force  is  exerted  on  the  object,  causing  vibration.  This  vibra¬ 
tion  is  produced  by  direct  conversion  of  ultrasound  energy  to  low  frequencyelastic  wavesin 
the  object.  This  conversion  is  a  result  of  momentum  change  by  the  object  Noteworthy,  ne¬ 
ther  the  propagation  medium  nor  the  object  material  need  to  be  acoustically  nonlinear 

^USVASemploITthe  radiation  force  of  two  intersecting  continuous  ultrasound  beamsto 
remotely  vibrate  the  object  at  an  arbitrary  low  frequency.  Object  vibration  produces  a  sou 
field  that  we  named  ‘acoustic  emission’  in  the  medium,  which  is  a  function  of  object  me- 
chanica? properties.  The  term  ‘acoustic  emission’  is  used  here  to  describe  the  acoustic  field 
in  response  to  a  cyclic  vibration  of  the  object,  and  should  not  be  confused  with  similartermi 
nology  used  in  the  field  of  nondestructive  testing  of  materials  or  in  opto-acoustic  imagmg 
context,  where  ‘acoustic  emission’  is  used  to  describe  the  acoustic  field  resulting  from  s  rnc- 
tural  deformation,  cracking  or  thermal  expansion  of  the  object.  By '  e  * 

emission  field,  one  can  estimate  some  of  the  mechanical  parameters  of  the  object,  including 

the  elastic  constants. 
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FIG.  1  USVAS  system. 


Figure  1  illustrates  a  USVAS  system.  The  ultrasound  energy  is  produced  by  means  of  a 
two-element  confocal  ultrasound  annular  array  transducer,  consisting  of  a  center  disc  and  an 
outer  ring.  The  elements  are  driven  by  two  cw  sources,  at  frequencies  equal  to  /  and  f2  -fx  + 
A/,  where  these  frequencies  are  very  close  to  the  central  frequency  of  the  elements  (f0 ),  and  A f 
<  <f0.  The  beams  interact  only  in  a  small  region  near  the  joint  focal  point.  The  object  is 
placed  at  the  joint  focal  point.  The  change  in  the  energy  density  due  to  the  reflection  and  ab¬ 
sorption  of  ultrasound  by  the  object  produces  a  radiation  force  component  in  the  beam  direc¬ 
tion  on  the  object.7  This  force  oscillates  at  frequency  A/  in  response  to  the  temporal  changes 
in  the  energy  density  caused  by  pressure  amplitude  modulation.  The  oscillatory  radiation  force, 
in  turn,  vibrates  the  object  at  A/'.  Object  vibration  elicits  a  sound  field  in  the  medium,  which 
is  received  by  a  nearby  microphone  (or  hydrophone)  and  recorded  after  digitization.  Referring 
to  reference  5,  the  amplitude  of  the  acoustic  emission  signal  <t>(A/)  can  be  written  in  terms  of 
the  frequency  response  of  the  object  Q(Af)  and  the  transfer  function  of  the  medium  //(A/j: 

<X>{Af)  =  CdrlI{Af)Q{Af)  (2) 

where  C  is  a  constant  proportional  to  the  intensity  of  the  impinging  ultrasound  beams. 

1  lence,  with  some  knowledge  of  If  (A /),  0(A/‘)  can  be  used  to  estimate  Q(Af).  This  can  be 
accomplished  by  sweeping  Af  in  the  range  of  interest  while  keeping  the  amplitude  of  the  im¬ 
pinging  ultrasound  fields  constant.  The  resonance  frequency(ies)  of  the  object  can  be  deter¬ 
mined  from  the  peak(s)  of  Q{Af).  The  positions  of  such  peaks  arc  independent  of  the 
constants  C’and  clr.  In  some  applications,  //(A f)  may  not  be  known  or  difficult  to  determine, 
liven,  under  these  conditions  we  may  still  be  able  to  use  our  method.  That  is  because  we  only 
need  to  measure  the  resonance  frequencies  oi'Q(Af)  —  not  the  entire  function  of  Q{Af).  In 
many  applications,  the  object  is  made  to  exhibit  a  sharp  (high  Q,  or  small  relative  bandwidth) 
resonances,  which  are  indicated  by  sharp  peaks  in  Q{Af ),  and  hence,  in  0(A/).  In  such  cases, 
the  positions  of  these  peaks  may  not  be  sensitive  to  H(Af).  Therefore,  the  resonance  frequen¬ 
cies  can  be  found  directly  from  d>(A/’),  without  a  complete  knowledge  of  M{Af).  For  exam¬ 
ple.  assuming  that  we  know  II(Af)  docs  not  have  a  resonance  within  the  frequency  range  of 
interest,  then  we  can  locate  the  resonance  frequencies  of  the  test  object  simply  by  determin¬ 
ing  the  frequency  values  at  sharp  peaks  of  <t>(A/~). 


Estimation  of  Young’s  modulus  using  the  resonance  frequency 

To  compress  a  rod  (length  /,  square  cross  section,  b  x  b)  by  an  amount  dl  requires  a  force 
EJ>\dbl).  where  the  constant  Eu  is  called  the  Young’s  modulus8  and  can  be  expressed  in  units 
of  pressure.  The  Young’s  modulus  of  a  material  can  be  found  by  measuring  the  resonance 
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frequency  of  a  rod  made  from  that  material.  The  resonance  frequency  depends  on  the  bound¬ 
ary  conditions  of  the  rod,  e.g.,  clamped  on  one  end  and  free  at  the  other,  clamped  at  both  ends, 
ox  free  at  both  ends.  A  simple  choice  is  to  resonate  a  rod  that  is  clamped  at  one  end  and  free  at 
the  other.  For  this  case,  the  Y oung’s  modulus  of  the  rod,  E0 ,  can  be  calculated  from  the  fol- 
lowing  relation:8 


(3) 


where  fR  is  the  fundamental  resonance  frequency  of  the  rod  in  vacuum  and  p  is  its  mass  den¬ 
sity.  In  MKS  units,  E0  is  expressed  in  Newtons  per  square  meter  (or  the  pressure  unit  Pascal, 
Pa),  /  and  b  in  meters,  p  in  kilograms  per  cubic  meter  and/*  in  Hz. 

In  practice,  it  is  simpler  to  measure  the  resonance  frequency  of  the  rod  in  a  fluid,  like  water. 
The  resonance  frequency  in  water  is  lower  than  in  vacuum  because  of  the  loading  effect  of 
water  mass.  The  resonance  frequency  in  vacuum  fR  can  be  calculated  in  terms  of  the  reso¬ 
nance  frequency  in  water  (/)  as:9 


fR  =  f„ -Jl  +  U962pw  hn  (4) 

where  m  is  the  mass  per  unit  length  of  the  rod,  in  kilograms  per  meter,  and  pw  is  water  density, 
in  kilograms  per  cubic  meter.  The  above  derivation  is  not  limited  to  square  rods,  similar  rela¬ 
tionships  can  be  derived  for  nonsquare  rods. 


Excitation  of  a  rod  by  the  radiation  force 

Consider  a  rod,  clamped  at  the  bottom  and  free  at  the  other  end,  placed  in  water.  T o  vibrate 
the  rod,  we  focus  a  modulated  ultrasound  beam  on  the  free  tip  of  the  rod,  with  the  beam  axis 
in  the  direction  of  expected  deflection  of  the  rod.  For  a  rod  with  a  square  cross  section,  we 
point  the  beam  normal  to  one  of  the  rod  sides.  The  radiation  force  of  the  beam  is  applied  to 
the  rod  via  water  molecules  in  the  projected  area.  This  force  vibrates  the  rod;  hence,  the  ul¬ 
trasound  energy  is  directly  converted  to  vibrational  energy  in  the  rod.  In  contrast  to  some 
other  resonance  methods  that  use  an  attached  actuator,  this  method  excites  the  rod  without 
any  contact  (except  the  surrounding  liquid).  Hence,  there  is  no  additional  boundary  condi¬ 
tion  imposed  on  the  rod.  The  rod  is  free  to  vibrate  (except  at  the  clamped  end)  at  its  natural 
modes.  Therefore,  the  position  and  the  size  of  the  excitation  point  along  the  length  of  the  rod 
is  immaterial,  as  long  as  it  delivers  the  energy  to  the  rod  and  the  force  is  in  the  correct  orienta¬ 
tion.  We  also  note  that  absorption  of  ultrasound  by  water  can  produce  streaming.7  Streaming 
can  potentially  be  another  source  of  excitation  to  the  rod.  Because  water  absorption  is  low, 
we  neglect  this  effect. 


EXPERIMENT 

The  experimental  setup  is  shown  in  figure  1.  The  free  end  of  an  aluminum  alloy  (Alloy 
7129,  Reynolds  Metals  Company,  Richmond,  VA)  rod  ( /  =  6  cm,  b  =  9  mm)  was  placed  at 
the  focal  plane  of  the  ultrasound  beams  in  a  water  tank.  The  other  end  of  the  rod  was  screw 
clamped  solidly  to  a  heavy  articulated  frame  that  was  rigid  enough  for  the  test.  Water  tem¬ 
perature  was  1 8. 5°C  (not  controlled). 

The  cc.ifocal  transducer  is  constructed  using  a  spherical  piezoelectric  cap.  The  two  ele¬ 
ments  are  constructed  by  dividing  the  back  electrode  of  the  piezoelectric  wafer,  with  the  ra¬ 
dius  of  curvature  to  70  mm,  into  a  central  disc  and  the  outer  ring,  such  that  the  elements  have 
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FIG.  2  Normalized  frequency  response  of  the  rod  in  water. 


identical  beam  axes  and  focal  lengths.  Radii  of  the  elements  are:  a,  =  14.8  mm  (inner  disc),  a2 
~22.5  mm  (outer  radius  of  the  ring),  and  16.8  mm  (inner  radius  of  the  ring).  The  focal 
distance  is  70  mm,  and  the  center  frequency,  /,  is  3  MHz.  The  two  beams  emitted  by  this 
transducer  interact  at  the  focal  region,  producing  an  oscillatory  radiation  stress  on  the  object 
in  a  circular  area  with  a  diameter  of  approximately  700  pm.5  For  a  more  detailed  description 
of  the  stress  field  distribution  on  the  target,  refer  to  reference  10. 

Transducer  elements  were  driven  by  two  stable  rf  synthesizers  (HP  33 120 A  and  Analogic 
2045)  at  frequencies  of/  =  3  MHz  and  f2  =/+  A f  where  A f  was  swept  slowly  (in  four  sec¬ 
onds)  from  1,600  to  1,800  Hz.  The  ultrasound  intensity  at  the  focal  area  was  about  3.8xl05 
W/m2.  Sound  produced  by  the  object  vibration  was  detected  by  a  submerged  audio 
hydrophone  (1TC  model  680,  sensitivity:  -154  dB  re  lV/pPa)  placed  within  the  water  tank. 
The  received  signal  was  filtered  and  amplified  by  a  programmable  filter  (Stanford  Research 
Systems,  SR650)  to  reject  the  noise,  then  digitized  by  a  12  bits/sample  digitizer  (HP  model 
E 1 429 A)  at  20,000  samples  per  second.  A  total  of  80,000  samples  of  acoustic  emission  data 
were  recorded  on  a  computer  disc  (Sun  Sparc  computer),  and  the  spectrum  was  calculated 
using  fast  Fourier  transform. 

To  ensure  that  the  acoustics  of  the  experiment  set-up,  including  the  water  tank,  was  not  af¬ 
fecting  the  estimation  of  the  rod  resonance  frequency,  we  made  an  additional  measurement 
with  the  rod  replaced  by  a  point  object,  consisting  of  a  small  glass  bead  (diameter<  0.5  mm), 
and  obtained  the  resulting  spectrum  in  the  same  frequency  range  (1,600-1,800  Hz).  The 
glass  bead  was  placed  on  a  thin  latex  sheet.  The  latex  sheet  is  almost  transparent  to  the  ultra¬ 
sound,  hence  the  radiation  force  on  this  sheet  is  negligible.  The  bead,  however,  presents  a 
strong  scatterer  and  acts  as  a  point  source  when  this  sheet  is  vibrated  by  the  radiation  force. 
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RESULTS 

Figure  2  shows  the  frequency  response  of  the  rod  (acoustic  emission  amplitude  versus  fre¬ 
quency).  The  sharp  peak  indicating  the  fundamental  resonance  frequency  of  the  rod  in  wa¬ 
ter, /„  =  1 ,7 1 2  Hz.  Referring  to  Eq.  (4),  the  corresponding  resonance  frequency  in  vacuum  is 
/  =  2,051  Hz.  Finally,  using  Eq.  (3),  the  Young’s  modulus  of  the  aluminum  alloy  rod  is 
found  to  be:  E„  =  6.97  xlO'°  N/m‘  (or  Pa).  This  value  is  in  agreement  with  the  value  pub¬ 
lished"  by  the  Aluminum  Association  for  this  aluminum  alloy  (£„=  7  xlO'°Pa),  which  is  also 
close  to  the  value  for  aluminum  (  E„  =  7  xlO10  Pa )  published  in  reference  1 1  The  spectrum 
obtained  from  the  glass  bead  showed  no  peaks  in  the  frequency  range  of  interest,  and  hence, 
validated  our  assumption  that  our  measurement  of  was  not  influenced  by  the  frequency  re¬ 
sponse  of  the  propagation  medium. 


DISCUSSION 


The  aim  of  this  paper  is  to  introduce  a  novel,  noncontact  method  for  evaluation  of  the  dy¬ 
namic  Young’s  modulus  of  materials.  This  method  employs  the  radiation  force  of  ultra¬ 
sound  according  to  the  US  V  AS  technique  to  vibrate  the  object  into  its  resonance.  Hence,  this 
method  can  be  viewed  as  a  resonance  technique.  Here,  we  discuss  some  features,  advan- 
lages,  limitations  and  applications  of  the  method. 


Measurement  accuracy 

Because  the  operation  of  the  method  is  based  on  material  resonance,  its  accuracy  in  mea- 
surement  of  the  Young’s  modulus  is  basically  similar  to  those  of  other  resonance  methods 
used  for  this  purpose.  In  particular,  since  £0  is  proportional  to  the  fourth  power  of  /  (Eq.  (3)), 
accurate  measurement  of  this  parameter  is  of  critical  importance  in  the  accuracy  of  the  fina 
result.  The  accuracy  of  this  method  can  be  improved  by  taking  special  care  to  measure  the 
geometrical  and  frequency  parameters  with  higher  precision.  However,  we  have  not  at¬ 
tempted  to  discuss  these  issues  here.  _ 

In  this  paper,  we  used  only  one  resonance  frequency  to  determine  the  Young  s  modulus  of 
the  object.  One  can  increase  the  accuracy  of  this  method  by  measuring  resonance  frequen¬ 
cies  in  more  than  one  mode.  In  that  way,  we  have  more  data  that  relate  to  the  same  elastic  pa¬ 
rameter.  This  approach,  in  general,  can  lead  to  higher  accuracy  of  the  final  result. 


Advantages 

The  main  advantage  of  this  method  over  other  resonance  methods  is  its  capability  of 
noncontact  and  remote  excitation  and  detection.  This  can  be  desirable  in  a  variety  of  applica¬ 
tions,  especially  when  the  source  of  excitation  and  the  detector  must  be  set  at  some  distance 
away  from  the  resonator  e.g.,  when  measuring  the  elastic  parameters  of  a  material  at  high 

temperatures).  ,  . 

Another  advantage  of  this  method  is  its  high  sensitivity.  Detection  of  object  motion  using 

the  acoustic  emission  field  can  be  very  sensitive.  It  is  shown  that  the  hydrophone  used  in  the 
experiment  is  sensitive  enough  to  detect  the  acoustic  emission  due  to  oscillatory  disp  ace- 
ments  of  the  order  of  a  few  nanometers.5  An  implication  of  high  sensitivity  is  that  it  allows 
the  resonance  to  operate  at  its  linear  regime.  It  is  well  known  that,  at  sufficiently  high  ampli¬ 
tudes,  most  resonators  exhibit  nonlinear  behavior.12-15  Two  common  nonlinear  phenomena 
arc:  a  transition  from  sinusoidal  to  nonsinusoidal  motion,  which  leads  to  harmonic  genera- 
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tion,  and  a  shift  in  the  modal  frequency.  Modal  frequencies  can  shift  upwards  or  downwards 
at  large  amplitudes.  If  the  effective  spring  constant  of  the  resonator  increases  with  amp  1- 
tude,  the  frequency  rises;  if  the  effective  spring  constant  decreases  with  amplitude,  the  fre¬ 
quency  falls.  This  phenomena  can,  in  general,  lead  to  errors  in  estimation  of  the  elastic 
constants  by  resonant  methods.  High  sensitivity  of  motion  detection  in  a  USVAS  system  al¬ 
lows  one  to  test  the  resonator  at  its  linear  regime  by  driving  it  at  low  oscillation  amplitudes, 

and  hence  avoid  the  errors  resulting  from  nonlinear  behavior. 

This  method  also  offershigh  geometrical  definition.  By  that,  we  mean  it  allows  us  to  exert 
a  dynamic  force  at  a  small  projection  region.  The  area  of  projection  region  depends  on  the 
characteristics  of  the  ultrasound  source,  namely  on  the  frequency  and  the  aperture  size,  av 
ing  a  small  projection  area  is  important  when  a  miniaturized  resonator  is  to  be  tested,  t  he 
fact  that  the  excitation  is  applied  via  a  narrow  ultrasound  beam,  permits  system  adaptability 

to  a  variety  of  resonator  shapes.  4 

Another  advantage  of  USVAS  is  its  relative  insensitivity  to  the  constituent  material  ot  the 

resonator.  Since  the  method  does  notrely  on  electrical  or  magnetic  properties  of  the  resona¬ 
tor,  it  can  be  used  to  test  a  wide  range  of  materials. 

Practical  issues 

There  are  some  issues  to  be  considered  when  using  USVAS  to  estimate  the  Young’s 
modulus  as  proposed  here.  The  frequency  spectrum  of  the  acoustic  emission  signal  is  a  prod¬ 
uct  of  the  object  frequency  response  and  the  frequency  response  of  the  propagation  medium 
(including  the  water  tank).5  In  order  to  assume  that  the  latter  would  not  have  a  major  effect  in 
determining  the  resonance  frequency  of  the  object,  one  must  arrange  the  propagation  me¬ 
dium  to  have  a  relatively  flat  frequency  response  in  the  vicinity  of  the  object  resonance  re- 
quenev.  Also  a  high  Q  (small  relative  bandwidth)  resonating  object  would  improve 
measurements  precision.  Temperature  variations  can  change  the  Young’s  modulus  of  the 
object3  and  the  effect  of  water  loading  by  changing  the  viscosity  .  Therefore,  the  temperature 
needs  to  be  controlled  for  more  accurate  estimation. 


Limitations 

The  method  presented  here  relies  on  the  acoustic  coupling  of  the  medium  to  the  resonator, 
both  for  excitation  and  detection.  Hence,  the  acoustic  properties  of  the  medium  surrounding 
the  resonator  are  important  in  the  success  of  the  method.  For  example,  the  method  can  not 
operate  in  a  vacuum.  The  method  also  fails  if  the  acoustic  path  from  the  ultrasound  source  to 
the  resonator  is  not  suitable  for  directing  the  ultrasound  beam.  In  practice,  this  might  happen 
due  to  the  existence  of  obstacles  or  high  attenuation  in  the  beam  path. 

Further  applications 

The  method  can  also  be  used  to  study  the  variations  of  the  elastic  parameters  with  temper¬ 
ature.  For  this  purpose,  one  can  expand  the  method  by  adding  a  phase-Iocked-loop  to  t  e 
system  to  track  the  variations  of  the  resonance  frequency  caused  by  temperature  changes,  in 
a  temperature-controlled  environment.  The  fact  that  the  ultrasound  source  and  the  receiving 
hydrophone  can  be,  in  principle,  thermally  isolated  while  acoustically  coupled  to  the  object, 
provides  the  possibility  of  measuring  the  elastic  moduli  in  a  wide  temperature  range 

In  some  applications,  air  coupling  may  be  preferred  over  liquid  coupling,  for  example,  to 
reduce  the  loading  effects  of  the  coupling  medium  or  to  measure  the  elastic  parameters  at 
high  temperatures.  In  such  cases,  airborne  ultrasound  may  be  used.  However,  one  must  take 
into  account  the  high  attenuation  and  low  sound  speed  of  air  when  designing  the  system. 
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Taking  advantage  of  the  fact  that  soft  tissues  are  good  sound  conductors,  one  may  extend 
the  principle  presented  here  for  biomedical  applications,  such  as  bone  characterization.  A 
modified  version  of  the  system  may  be  used  to  measure  the  frequency  response  of  the  bone  at 
low  frequencies  and  evaluate  its  elastic  parameters. 


CONCLUSIONS 


USVAS  promises  a  new  way  for  noncontact  and  remote  measurement  of  the  Young’s 
modulus  of  materials.  An  airborne  version  of  USVAS  would  provide  further  flexibility  in 
applications  where  water  coupling  is  not  possible. 
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B-5  12:00  noon 

POINT  RESPONSE  OF  VIBRO-ACOUSTOGRAPHY 

SYSTEM 

M.  FATEMI*  and  J.  F.  GREENLEAF,  Mayo  Clinic  and  Foundation,  Rochester, 
MN. 

Vibro-acoustography  is  a  method  for  imaging  objects  in  terms  of  their  response 
to  a  low  frequency  vibrational  force  [Science,  280(4):82-85,  April  1998].  The 
vibrational  force  is  produced  by  projecting  an  intensity  modulated  ultrasound 
beam  on  the  object.  The  resulting  radiation  stress  is  highly  localized  and  suf¬ 
ficiently  strong  to  vibrate  the  object  at  the  modulation  frequency.  Response 
of  the  system  to  a  point  object  is  of  interest  because  it  can  be  used  to  evalu¬ 
ate  image  resolution.  The  point  response  of  a  vibro-acoustography  system  is  a 
function  of  the  ultrasound  beam  shape.  Flere,  we  present  system  point  response 
for  three  different  beamforming  methods:  (1)  single  amplitude  modulated  beam; 
(2)  x-focal  method,  in  which  two  crossing  single  frequency  ultrasound  beams  are 
used  to  generate  a  modulated  field;  and  (3)  confocal  method  which  uses  two 
coaxial  single  frequency  beams.  To  compare  the  point  responses,  we  considered 
the  ultrasound  sources  used  in  these  methods  to  be  of  equal  areas.  Theoretical 
results  show  that  the  confocal  beam  system  has  the  smallest  mainlobe  width. 
The  x-focal  method  offers  the  best  sidelobe  attenuation.  The  magnitude  of  the 
first  sidelobes  for  single  beam,  confocal,  and  the  x-focal  systems  are  17.6  dB, 
19.5  dB,  and  34.7  dB,  below  their  respective  mainlobes.  Experimental  results 
are  in  agreement  with  the  theory. 

This  work  was  supported  by  grant  DAMD17-98-1-8121  from  the  Army  Medical  Research  and 
Material  Command. 
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COMPLEX-VALUED  ABSOLUTE  STIFFNESS 
ESTIMATION  USING  DYNAMIC  DISPLACEMENT 
MEAUREMENTS  AND  LOCAL  INVERSION  OF 
CONSERVATION  OF  MOMENTUM 
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Ultrasound  Research,  Mayo  Foundation,  Rochester,  MN,  USA,  and  2Magnetic 
Resonance  Research,  Mayo  Foundation,  Rochester,  MN,  USA. 

| 

It  has  been  shown  that  elasticity  can  be  estimated  from  direct  displacement 
measurements  measured  using  either  ultrasonic  or  magnetic  resonance  (MRE) 
elastography  methods.  Most  of  the  inversion  techniques  presented  in  the  litera¬ 
ture  are  based  on  inverting  a  displacement  dataset  that  satisfies  a  static  equation 
of  shear  motion  resulting  in  only  relative  estimates  of  shear  modulus.  We  present 
a  robust  method  for  estimating  absolute  shear-wave  speed  and  attenuation  from 
dynamic  displacement  data  based  on  local  inversion  of  the  linearized  differential 
equation  satisfied  by  small-amplitude  displacements  in  a  source-free,  possibly  I 
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INTRODUCTION 

It  is  well  known  that  changes  in  elasticity  of  soft  tissues  are  often  related  to  pathol¬ 
ogy.  Traditionally,  physicians  use  palpation  as  a  simple  method  for  estimating  mechanical 
properties  of  tissue.  The  dynamic  response  of  soft  tissue  to  an  oscillating  or  impulsive  force 
is  also  valuable  in  medical  diagnosis.  For  instance,  rebound  of  tissue  upon  sudden  release  of 
pressure  exerted  by  the  physician’s  finger  on  the  skin  provides  useful  diagnostic  information 
about  the  tissue.  In  vibro-acoustography  ultrasound  is  used  to  produce  an  oscillating  radia¬ 
tion  force  which  then  produces  an  acoustic  response  from  tissues.  Specifically  acoustic 
energy  is  emitted  from  tissues  in  response  to  an  oscillatory  radiation  force  produced  by 
interfering  focused  beams  of  ultrasound.  Frequency  spectra  of  ultrasound  stimulated  acous¬ 
tic  emission  exhibit  object  resonances.  Raster  scanning  the  radiation  force  over  the  object 
and  recording  the  amplitude  and  phase  of  the  emitted  sound  results  in  data  from  which 
speckle  free  images  related  to  elastic  compositions  of  the  acoustically  emitting  objects  can 
be  computed.  We  show  here  that  intact  and  broken  struts  on  prosthetic  disk  heart  valves 
(convexo-concave  Bjork-Shiley)  can  be  distinguished  by  their  resonant  response  to  vibro- 
acoustography. 


BACKGROUND 

Quantitative  measurement  of  the  mechanical  properties  of  tissues  and  their  display 
in  a  raster  format  is  the  aim  of  a  class  of  techniques  generally  called  elasticity  imaging  or 
elastography  (Gao  et  al.,  1996).  The  general  approach  is  to  measure  tissue  motion  caused  by 
an  external  (or,  in  some  methods,  internal)  force  or  displacement  and  use  it  to  reconstruct 
the  elastic  parameters  of  the  tissue.  The  excitation  stress  can  be  either  static,  or  dynamic 
(vibration).  Dynamic  excitation  is  of  particular  interest  because  it  provides  more  compre¬ 
hensive  information  about  tissue  properties  in  a  spectrum  of  frequencies.  Because  the 
motion  is  dynamic,  derivatives  of  displacement  in  both  space  and  time  can  be  measured 
allowing  the  computation  of  elastic  properties  with  low  dependence  on  boundary  condi- 


tions.  Magnetic  resonance  elastography  is  a  recently  developed  method  that  employs  a 
mechanical  actuator  to  vibrate  the  body  surface  and  then  measures  the  strain  waves  with 
phase  sensitive  Magnetic  Resonance  Imaging  resulting  in  a  method  termed  Magnetic  Reso¬ 
nance  Elastography  (Muthupillai  et  al.,  1995).  The  majority  of  elasticity  imaging  methods 
are  based  on  an  external  source  of  force,  resulting  in  a  spatially  wide  stress  field  distribution. 
This  requires  the  stress  field  to  pass  through  the  superficial  portion  of  the  object  before 
reaching  the  region  of  interest  within  the  interior.  This  requirement  can  complicate  estima¬ 
tion  of  stiffness  because  the  stress  field  patterns  vary  with  depth.  An  alternative  strategy  is  to 
apply  a  stress  directly  to  the  localized  region  of  interest.  One  way  to  accomplish  this  is  to 
use  the  radiation  pressure  of  ultrasound.  An  example  is  the  use  of  ultrasound  radiation  pres¬ 
sure  to  vibrate  the  tissue  and  the  use  of  MRI  to  measure  the  displacement  (Sarvazyan  et  al., 
1998  and  Tau  et  al.,  2000).  The  approach  taken  in  this  paper  is  to  estimate  the  mode  of 
vibration  of  an  object  from  the  acoustic  emission  resulting  from  the  cyclic  displacement. 

Acoustic  radiation  force  is  the  time  average  force  exerted  by  an  acoustic  field  on  an 
object.  This  force  is  produced  by  a  change  in  the  energy  density  of  an  incident  acoustic  field 
(Beyer,  1978),  for  example,  due  to  absorption  or  reflection.  Several  advantages  may  result 
from  using  ultrasound  radiation  force  for  evaluating  tissue  stiffness,  including:  a)  acoustic 
(ultrasound)  energy  is  a  noninvasive  means  of  applying  force,  b)  current  ultrasound  technol¬ 
ogy  and  devices  can  be  readily  modified  for  this  purpose,  thus  eliminating  the  need  for 
developing  a  new  technology,  c)  radiation  force  can  be  applied  remotely  inside  tissue  with¬ 
out  disturbing  its  superficial  layers,  d)  the  radiation  stress  field  can  be  localized,  thus  provid¬ 
ing  for  precise  positioning  of  the  excitation  point,  and  e)  radiation  force  can  be  applied  in  a 
wide  range  of  frequencies  or  temporal  shapes.  These  features  make  radiation  force  methods 
highly  attractive  for  creating  remote  stress  fields  within  objects. 


VIBRO-ACOUSTOGRAPHY  THEORY 

Study  of  radiation  force  and  radiation  pressure  dates  back  nearly  one  century,  to  the 
time  of  Rayleigh  (1902  and  1905).  A  historical  review  of  radiation  force  and  radiation  pres¬ 
sure  is  presented  in  Beyer  (1978),  and  a  critical  review  of  the  subject  can  be  found  in  Chu 
and  Apfel  (1982).  Some  recent  developments  of  radiation  force/pressure  in  attenuating 
medium,  which  may  be  applicable  to  biological  tissues,  are  presented  in  Rudenko  et  al. 
(1996)  and  Jiang  and  Greenleaf  (1996).  Acoustic  radiation  force  is  a  universal  phenomenon 
in  any  wave  motion  that  introduces  some  type  of  unidirectional  force  on  absorbing  or 
reflecting  targets  in  the  wave  path.  Radiation  force  in  fluids  is  often  studied  in  the  context  of 
radiation  pressure.  Definitions  and  theoretical  developments  of  radiation  force  depend  on 
the  boundary  conditions.  Simple  explanations  of  these  definitions  can  be  presented  by  con¬ 
sidering  a  sound  traveling  inside,  and  along  the  axis  of,  a  cylindrical  container  toward  a  wall 
(Beyer,  1978).  Rayleigh  radiation  pressure  is  the  excess  pressure  produced  on  the  wall  when 
the  container’s  side  wall  is  confining  the  fluid  inside.  Langevin  radiation  pressure  is  the 
excess  pressure  on  the  wall  when  the  confining  side  wall  is  absent,  so  that  the  fluid  is  free  to 
move.  Here  we  will  focus  on  Langevin  radiation  pressure  because  the  conditions  for  which 
this  pressure  is  defined  apply  to  our  experimental  situation. 

It  can  be  shown  that  the  Langevin  radiation  pressure  of  a  plane  wave  impinging  nor¬ 
mally  on  a  perfectly  absorbing  wall  is  equal  to  the  total  energy  density  <E>,  where  <>  repre¬ 
sents  the  time  average.  If  the  wall  is  partially  reflecting,  this  pressure  would  be  equal 
to(l+R)<E>,  where  R  is  the  power  reflection  coefficient  (Beyer,  1978).  Thus,  in  general,  we 
can  write  the  radiation  force  of  a  normally  impinging  sound  beam  on  a  wall  as  F  =  dr(E), 


where  F  is  the  force  in  the  beam  direction,  and  dr  is  the  “radiation  force  function”  or  the 

drag  coefficient.  This  dimensionless  coefficient  is  defined  per  unit  incident  energy  density 
and  unit  projected  area.  For  a  planar  object,  F  is  numerically  equal  to  the  force  on  the  object. 
Physically,  the  drag  coefficient  represents  the  scattering  and  absorbing  properties  of  the 
object  (Westervelt,  1951).  For  a  perfectly  absorbing  object,  dr  =  1  and  for  a  perfectly 
reflecting  object  dr  =  2 .  In  the  case  of  oblique  incidence,  the  radiation  force  will  have  a 

normal  as  well  as  a  transverse  component.  A  more  detailed  description  is  presented  in  West¬ 
ervelt  (1951).  To  produce  a  time  varying  radiation  force,  the  intensity  of  the  incident  beam 
can  be  modulated  in  various  ways.  For  example,  a  short  ultrasound  pulse  can  produce  a  tran¬ 
sient  pulsed  radiation  force,  and  a  sinusoidally  modulated  beam  can  result  in  a  sinusoidally 
varying  force. 


METHODS 

We  interrogate  the  object  by  arranging  the  intersection  of  two  focused  continuous 
wave  (CW)  ultrasound  beams  of  different  frequencies  to  occur  at  a  selected  point  on  the 
object.  Interference  in  the  intersection  region  of  the  two  beams  produces  sinusoidal  modula¬ 
tion  of  the  ultrasound  energy  density.  Modulation  of  the  energy  density  creates  an  oscilla¬ 
tory  force,  effectively  vibrating  the  object  at  the  selected  region.  The  resulting  vibration  of 
the  object  produces  an  acoustic  field  (Fatemi  and  Greenleaf,  1998)  that  can  be  measured 
some  distance  away. 


Transducer 
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Figure  1.  Ultrasound  is  made  to  interfere  on  the  region  of  interest  and  the  resulting  acoustic  emission  is 
detected  with  a  hydrophone  or  microphone  to  estimate  the  response  of  the  object.  The  signal  is  digitized  for 
analysis. 


Ultrasound  beams  can  be  constructed  in  a  variety  of  ways  for  this  purpose  (Fatemi 
and  Greenleaf,  1999).  We  used  two  coaxial,  confocal  transducer  elements  of  a  spherically 
focused  annular  array  (consisting  of  a  central  disc  and  an  outer  annulus)  driven  by  two  CW 
signals  at  slightly  different  frequencies  /j  and  f2  (Figure  1).  The  energy  density  at  a  point  in 
this  ultrasound  field,  say  at  the  focus,  is  proportional  to  the  square  of  the  sum  of  the  ultra- 


sound  fields  from  the  two  elements.  Squaring  the  sum  of  two  sines  gives  rise  to  sum  and  dif¬ 
ference  frequency  terms.  Thus,  high  frequency  and  low  frequency  variations  in  energy 
density  result  at  the  intersection  of  the  two  beams  produced  by  the  two  elements.  Ultra- 
sound-stimulated  acoustic  emission  results  from  the  energy  term  that  produces  a  low-fre¬ 
quency  vibration.  The  low-frequency  force  on  a  target  at  the  focal  point  can  be  computed 

with  the  following  integration:  Ft  =  drj(E)da,  S  is  the  area  over  which  (E) ,  the  total 

s 

energy  density  in  the  focal  plane,  has  significant  value,  and  ( >  represents  a  short-term  time 
average.  For  focused  beams,  the  intersection  region  can  be  small  enough  that  Ft,  the  low 
frequency  force,  can  be  considered  to  be  an  oscillating  point  force  applied  to  the  object  at 
the  focal  intersection  of  the  beams.  To  produce  an  ultrasound  stimulated  vibro-acoustic 
emission  spectrogram,  we  vibrate  a  small  region  of  the  object  with  an  oscillating  radiation 
force  of  varying  frequency.  The  complex  amplitude  of  the  resulting  acoustic  emission  field, 

<3>(A if) ,  can  be  written:  ^(4 f)  =  drH{Af)Q(Af)  t  where  Q(Af)  is  a  complex  function  rep¬ 
resenting  the  mechanical  frequency  response,  or  admittance,  of  the  object  at  the  selected 
point  and  H(Af)  represents  the  combined  frequency  response,  or  transfer  function,  of  the 
propagation  medium  and  receiver  and  is  assumed  to  be  fixed  and  known  (Fatemi  and  Green- 
leaf,  1999).  Recording  0(A f)  allows  us  to  obtain  <2(4/)  for  each  point,  within  a  constant 
multiplier.  We  raster  scan  the  radiation  force  over  the  object  to  produce  data  that  can  be 
mapped  into  a  pictorial  format.  The  spatial  resolution  of  the  resulting  image  is  determined 
by  the  region  in  which  significant  interference  between  the  ultrasound  beams  occurs  and  is 
of  the  order  of  a  few  wavelengths  at  the  ultrasound  frequency.  In  a  slight  alteration  of  vibro- 
acoustography  called  vibrography  we  attach  an  accelerometer  to  the  object  to  measure  its 
vibration.  The  scan  of  the  vibration  inducing  ultrasound  beam  remains  the  same.  We  show 
both  types  of  images  here. 


RESULTS 

To  test  the  feasibility  of  using  the  technique  to  image  mechanical  properties  of  tis¬ 
sues,  we  measured  the  amplitude  of  acoustic  emission  from  intact  and  broken  struts  on 
Bjork-Shiley  convexo-concave  (BSCC)  prosthetic  heart  valves.  The  valves  were  scanned  in 
a  plane  perpendicular  to  the  ultrasound  beam  axis.  Vibro-acoustography  of  the  valve  at  40 

kHz  resulted  in  an  image  of  Q(40k)  and  is  illustrated  in  Figure  2.  Also  illustrated  in  Figure 
2  are  vibrography  images  representing  <2 (Ik)  of  a  valve  in  which  one  leg  of  the  small  strut 

is  broken.  One  can  see  that  the  strut  vibrates  differently  when  broken.  The  amplitude  images 
are  highly  detailed  and  exhibit  variations  in  acoustic  emission  from  mechanical  objects  hav¬ 
ing  different  vibration  characteristics.  These  differences  are  due  to  variations  in  the  product 
of  the  reflection  properties,  dr ,  and  the  effective  mechanical  vibration  admittance  properties 
of  the  object.  Thus,  both  vibro-acoustography  and  vibrography  are  similar  to  conventional 
pulse  echo  ultrasound  imaging,  which  is  sensitive  to  the  ultrasonic  parameters  of  the  object, 
but  these  two  new  imaging  methods  have  the  advantage  of  also  being  sensitive  to  the 
mechanical  admittance  at  low  frequencies. 


DISCUSSION 

Several  implications  can  be  deduced  from  the  results  of  this  study.  The  results  shown 
here  imply  that  fractured  struts  could  be  found  in  similar  types  of  valves  implanted  in 


patients  assuming  an  appropriate  clinical  instrument  can  be  built.  Vibration  modes  of  com¬ 
plex  objects  could  be  characterized  without  having  to  touch  the  object.  Motion  induced  by 
ultrasound  and  measured  with  ultrasound  pulse  echo  has  been  used  previously  to  study 
“hardness”  (Dutt  et  al.,  2000).  However,  the  sensitivity  of  ultrasound  pulse  echo  to  motion 
at  common  ultrasound  frequencies  is  limited  to  several  micrometers.  The  advantage  of  ultra¬ 
sound  stimulated  vibro-acoustic  emission  is  its  high  displacement  sensitivity.  Cyclic  dis¬ 
placement  of  100  nm  at  10  kHz  produces  an  acoustic  intensity  of  about  3.0  x  103  W/cm2. 
Hydrophones  similar  to  the  one  used  in  these  experiments  are  sensitive  to  as  little  as  10' 

15W/cm2,  and  therefore,  can  detect  very  small  cyclic  displacements.  The  method  will  be 
more  sensitive  for  higher  frequency  sound  because  acoustic  power  is  proportional  to  the 
square  of  frequency  for  constant  displacement  amplitude. 


Broken  strut  2.04  kHz  Intact  strut  7.4  kHz 


Figure  2.  Upper;  vibro-acoustography  image  of  artificial  heart  valve  (BSCC)  showing  high  detail  and  no 
speckle.  Highlights  on  strut  are  due  to  reflection  of  ultrasound  into  the  hydrophone  which  was  at  the  right. 
Lower;  images  obtained  with  vibrography  at  the  resonance  of  the  broken  (left)  and  intact  (right)  struts.  The 
strut  is  broken  at  the  arrow  and  the  free  end  vibrates  strongly  when  driven  at  frequency  of  2.04  kHz,  much 
lower  that  the  7.4  kHz  resonance  of  the  intact  strut  on  the  right. 


Ultrasound  stimulated  vibro-acoustic  spectrography  may  have  application  in  at  least 
two  general  areas.  The  first  is  nondestructive  evaluation  of  materials,  where  material  charac¬ 
teristics  and  structural  flaws  can  be  identified  by  measuring  changes  in  the  mechanical 
response  to  vibration  at  a  point.  The  object  under  test  could  be  remotely  probed  by  beams 
propagating  and  interfering  in  either  water  or  air.  Beams  propagating  within  the  object  could 
also  be  used  to  produce  acoustic  emission  from  flaws.  The  second  general  area  of  applica¬ 
tion  is  medical  imaging  and  detection.  Vibro-acoustography  appears  particularly  suitable  for 
noninvasive  detection  of  hard  tissue  inclusions,  such  as  imaging  arteries  with  calcification, 
detecting  breast  micro-calcifications,  visualizing  hard  tumors,  and  detecting  foreign  objects. 
Because  changes  of  stiffness  alter  the  vibration  frequency  response  or  damping  of  tissue, 
vibro-acoustography  can  potentially  provide  a  noninvasive,  remote,  high-resolution,  speckle 


free  “palpation”  technique  that  can  reach  small  abnormalities  that  are  otherwise  untouchable 
by  conventional  methods. 


SUMMARY 

In  an  in  vitro  experiment  the  acoustic  emission  characteristics  of  broken  and  intact 
struts  on  BSCC  valves  exhibited  different  frequency  responses  and  different  vibro-acoustog- 
raphy  images.  The  images  indicated  that  the  vibration  modes  of  the  object  could  be  deduced 
in  addition  to  the  structural  integrity  of  the  object.  Structural  integrity  of  such  prostheses 
could  be  deduced  using  vibro-acoustography  in  vivo. 
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2.  Project  Summary 

The  general  direction  of  our  research  activities  is  measurement  of  sound  and  vibration  in 
response  to  micro-force.  The  majority  of  our  projects  are  centered  on  the  hypothesis  that 
by  measuring  the  sound  and/or  vibration  resulting  from  such  small  forces  we  will  be  able 
to  obtain  important  information  about  the  object.  In  most  of  our  projects  this  micro-force 
is  produced  by  the  radiation  force  of  ultrasound,  which  is  in  mN  range,  and  the  resulting 
vibration  is  in  nanometer  or  angstrom  range. 

We  have  shown  that  by  measuring  the  acoustic  field  resulting  from  such  vibration  we  can 
obtain  estimates  of  the  mechanical  properties  of  the  object,  image  the  object  at  high 
resolutions,  and  detect  small  particles.  We  also  have  shown  that  the  radiation  force  of 
ultrasound  can  stimulate  the  fetus.  It  appears  that  this  field  of  science  has  great 
potential  in  material  science,  non-destructive  evaluations,  and  medicine. 

The  aim  of  this  proposal  is  to  acquire  instrumentation  that  would  enable  us  to  measure  low 
level  sound  and  vibration.  With  such  capabilities,  we  would  be  able  to  validate  our 
models,  measure  different  modes  of  vibration  and  differentiate  them,  and  improve  our 
measurement  techniques. 

This  proposal  is  for  the  acquisition  of  4  main  components:  (1)  an  anechoic  chamber,  (2)  a 
scanning  laser  vibrometer,  (3)  a  dual  channel  compact  fiber  vibrometer,  and  (4)  a 
de-aeration  system.  The  anechoic  chamber  allows  us  to  reduce  the  background  noise  and 
vibration  in  our  experiments.  Presently,  the  background  noise  is  our  main  limiting  factor 
in  our  newly  developed  imaging  method  called  vibro-acoustography.  The  scanning  laser 
vibrometer  allows  us  to  directly  measure  vibrations  in  angstrom  range.  Also  we  can  use 
this  device  to  obtain  the  spatial  distribution  of  motion  in  response  to  a  small  point 
force.  The  dual-channel  fiber  vibrometer  will  be  used  to  simultaneously  measure  vibration 
at  two  points  of  the  object.  These  instruments  would  enable  us  to  measure  the  transition 
time  and  the  phase  of  vibration.  Also  this  instrument  allows  us  to  measure  vibrations 
directly  under  water,  where  most  of  our  experiments  are  performed.  Finally,  the 
de-aeration  system  is  necessary  to  provide  clean  and  completely  degassed  water  for  our 
water  tank  experiments.  Noting  that  gas  bubbles  interfere  with  our  imaging  technique  and 
produce  a  significant  amount  of  noise  in  the  experimental  water  tank,  de-aeration  system 
would  be  a  valuable  tool  for  reducing  errors  in  most  of  our  experiments. 

Graduate  students,  research  trainees  and  fellows,  and  senior  investigators  will  use  the 
above  instrumentation.  Addition  of  these  equipment  to  our  laboratory  will  help  us  to 
attract  more  qualified  students  in  the  Biomedical  Imaging  track. 


